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The equilibrium constants of the reactions HX0+2D 
=D,.0+2H and H,O+D.0=2HDO are calculated and 
from these an expression is obtained for the electrolytic 
separation coefficient, valid when zero-point energy and 
tunnelling are negligible in the activated state. The change 
in this coefficient with increase in concentration of the 
heavy isotope is calculated. Experimental values have 
been determined for the separation coefficient at different 
metallic cathodes, ranging from 7.6 to 2.8. Local concen- 


tration disturbances round the cathode are not in general 
important, but may be the cause of the poor separation 
found with platinum black. The importance of half-quanta 
in the activated state is discussed. The experimental 
results considered in the light of the calculations indicate 
that tunnelling is not very important in the electrolytic 
separation. The vapor pressures of liquids containing 
isotopes are considered. 





SECTION 1. INTRODUCTION 


HE separation of the hydrogen isotopes by 
electrolysis of solutions of electrolytes in 
water depends essentially upon the removal of 
hydrogen from water molecules. The rapidity of 
the equilibrium reaction 2H,O=H;0++OH- 


_ compels us to consider the water as the funda- 


mental source of the hydrogen evolved at the 
cathode. 

Reactions at the anode come into question only 
insofar as the products of the anodic reactions 
diffuse to the cathode and modify the com- 
position of the solution there. Confining attention 
to solutions which evolve hydrogen and oxygen 
in equivalent quantities, the only changes in the 
nature of the electrolyte during electrolysis are 
the increase in electrolyte concentration and the 
change in the relative number of H and D atoms. 
The effect of local changes in composition round 
the cathode depends upon the experimental con- 
ditions; this point is taken up in Section 3, where 
it is shown that under the usual conditions of 
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experiment they may for practical purposes be 
ignored. 
The “electrolytic separation coefficient’’ @ is 
defined by 
d log H=ad log D. (1) 


When the mol fraction of electrolyte exceeds 
the atomic fraction of the isotope present in 
lesser amount, then the nature and concentration 
of the electrolyte will have some effect upon a. 
Experiment shows, however, that this effect is 
not a very large one. 


SECTION 2. THEORY OF THE ELECTROLYTIC 
SEPARATION 


In this section it is assumed that the solution 
surrounding the cathode is in thermodynamic 
equilibrium with the bulk of the solution (cf. 
Section 3) and that the mol fractions of all ionic 
species are small compared with the atomic 
fractions of both H and D. 

The mechanism of the process at the cathode 
is not always the same in acid and alkaline 
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solutions. For example, some of the hydrogen 
obtainable by electrolyzing a KOH solution with 
a mercury cathode is produced in a secondary 
reaction between the potassium amalgam and 
the water; similarly at a lead cathode it is 
usually supposed that a K-Pb alloy is formed 
initially, and at cathodes which show no tendency 
to form alloys with alkali metals, the metallic 
cation may nevertheless be preferentially dis- 
charged if the hydrogen overvoltage is high 
enough. In solutions of strong acids, direct dis- 
charge of hydrogen ions is the only possibility. 
The various theories of overvoltage proposed 
have all been designed primarily to explain the 
empirical relation between current density and 
overvoltage. The main difference between them 
concerns the decision as to which of the steps 
intervening between the original water molecules 
and the final hydrogen molecules is the slowest 


step. 
We have to consider the consecutive stages: 
(a) 2H.O0 = H,Ot++OH-, 
(b) H,0++ M-=H.0+H-M, 
(c) 2H—- M=H.2+2”M. 


Stage (b) is the neutralization of a hydrated 
hydrogen ion by the negatively charged metal 
M, with formation of a hydrogen atom attached 
in some way to the surface atom M. Stage (c) 
represents the desorption of hydrogen in the 
molecular form from the metal surface. 

One theory of overvoltage assumes that (c) is 
the slow process, (a) and (b) being reversible and 
in equilibrium. This theory (the essentials of 
which were originally due to Tafel') has since 
been modified in various directions, but for the 
discussion in this section, it is only important 
that (a) and (b) are assumed to be in equilibrium. 

The well-known theories of Erdey-Gruz and 
Volmer? and of Gurney* both place the slow 
reaction at stage (b); they differ fundamentally 
in that the former associates the slowness with 
the passage of protons over a potential barrier 
lying midway across the Helmholtz double layer, 


1 Tafel, Zeits. f. physik. Chemie 34, 200 (1900); 50, 641 
(1905). 

2 Erdey-Gruz and Volmer, Zeits. f. physik. Chemie 
A150, 209 (1930). 
3 Gurney, Proc. Roy. Soc. A134, 137 (1932). 
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and the latter attributes it to the slow penetra- 
tion of electrons through a potential barrier of a 
different kind situated between the metal and 
the ions in the solution. 

In the remainder of this paper (Part I) we 
discuss a number of points which are relevant 
to any mechanism in which the slow process 
involved the passage of hydrogen particles over 
or through the upper part of a potential barrier. 
We are not here thinking of a mechanism involv- 
ing the slow leakage of the electrons through the 
potential barrier connected with the work 
function. 

Corresponding to any special theory of the 
mechanism of the cathode process there is some 
definite activated state with which the over- 
voltage is correlated. Consider first the simplest 
case, where the activated state consists of entirely 
free hydrogen atoms. The reaction at the cathode 
cannot really involve such free atoms because the 
activation energy would then be prohibitively 
large, but the calculation forms a basis on which 
the consideration of mechanisms involving dif- 
ferent activated states can be built; this cal- 
culation also applies to any activated state which 
fulfills two conditions: (a) that the total zero 
point energy in the activated state is zero, (b) 
that quantum mechanical leakage of the hy- 
drogen particles through the potential barrier is 
negligible in comparison with the rate at which 
they pass over the barrier. By combining the 
result of this calculation with experimental 
values for the electrolytic separation coefficient 
we have a new possibility of deducing properties 
of the activated state. 

We make the usual assumption for reaction 
rates,‘ that the rate of the slow process is the 
same when the reaction is proceeding irrevers- 
ibly, in the forward direction only, as the rate 
in either one direction at equilibrium. This rate 
is the number of particles in the activated state 
at equilibrium multiplied by the average velocity 
with which they travel in the forward direction. 
Since for our present problem we require only 
the relative reaction rates of the two isotopes, 
we first calculate the equilibrium ratio of H to D 
atoms (and H+ to Dt ions) in free space, if 
equilibrium with water containing H,O, HDO 
and D.O. 


4 Tolman, Statistical Mechanics, p. 269 (1926). 





























































with 
and 


this 
used 
in O1 


W 


wher 
phas 


wher 
the z 
lowe: 
B,, ¢ 
the » 
tion 

are t 


Th 
simil. 
and 4 
Since 
istic 
pract 
even 


K,=. 


the st 
For 
three 


Th 
With 
tion f 








tra- 
of a 
and 


we 
ant 
ess 
ver 
‘ier, 
alv- 
the 
ork 


the 
yme 
ver- 
lest 
rely 
ode 
the 
vely 
rich 
dif- 
cal- 
rich 
zer0 


(b) 


or is 
rich 
the 
ntal 
ient 
‘ties 


tion 

the 
ers- 
rate 
rate 
tate 
city 
jon. 
only 
pes, 
oD 
, in 


[DO 





SEPARATION OF THE 
Consider the equilibria 
H.0+2D = D,.0+2H, 
H.0+D.0 =2HDO, 


(1a) 
(1b) 


with equilibrium constants in the gas phase K, 
and Ky, respectively. 

In the equations which follow, and throughout 
this paper, the subscripts 1, 2, 0, 3, 4 will be 
used to indicate the following molecular species 
in order: H2O, D2.O, HDO, H and D. 

We write 

Ku= F.F3/F\F 2, 
where the partition function for H2O in the gas 
phase is 
+o) (2am,kT)} 821°(87°A 1B,C;) 
F,=exp 
kT h’ 2h 
1 


X(RkT)! II x4, (2) 
i=j,, ky, ty 1—exp (—hv;/kT) 








where —¢,=energy of H.O in its lowest state, 
the zero of energy being the three atoms in their 
lowest states; m;=mass of H2O molecule; A, 
B,, Ci, =moments of inertia of HzO molecule; 
the v’s with subscripts j;, k:, /:, denote the vibra- 
tion frequencies of the H2O molecule; h and k 
are the Planck and Boltzmann constants. 


Fy=[(2am3kT)!/h? JA. (3) 


The partition functions Fp,o and Fp are 
similar, except that the statistical weights 4 
and 4 become 9 and 6 for D,O and D, respectively. 
Since, with the vibration frequencies character- 
istic of D,O the quantities exp (—hv/kt) are 
practically negligible in comparison with unity, 
even at 100°C, we find for K, 


K (= ~) (~) ©) (=) ; (4) 
e=€X — —— eee : 
. kT mJ \mJ \AiBiG 


the subscripts having the meanings already given. 
For the equilibrium in the gas phase of the 
three kinds of water molecules we have 


Ky= Fe?/ FFs. (S) 





The partition function Fy is the same as F; 
with the substitution of the appropriate vibra- 
tion frequencies, and the mass mo, energy — fo, 


HYDROGEN 
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and moments of inertia Ao, Bo, Co, of the HDO 
molecule, except that the symmetry number 2 
is omitted and the statistical weight is 6 instead 
of 4. We then have for K, 








(+2f0-—1—¢2) ( me ) 
K,=exp 
kT mM iM 
A oBoCo 
x x4. (6) 
(A1B,CiA2B2C2)} 


The quantities —{1, —{»9 and —f~2 enter into 
Eqs. (4) and (6) as differences, and since the 
potential function is the same for isotopic 
molecules,* we have only to consider the dif- 
ferences in the zero-point energies involved, i.e., 
the half-quanta for the three normal modes of 
vibration of each molecule. For the molecule 
H:.O we have taken the values obtained by 
Mecke’ and his collaborators in their analysis of 
the infrared absorption of water vapor. 

To obtain the frequencies for the isotopic 
molecules HDO and D,O we have made use of 
the very helpful statement of the problem of the 
three particle system recently published by Van 
Vleck and Cross.6 These authors also carried 


* We are here neglecting the change in the force con- 
stants in the water molecule caused by the mass difference 
between H and D. The effect can be estimated on the 
basis of the Heitler-London theory of binding. For a 
hydrogen atom the potential energy with respect to the 
oxygen atom is E=(1/dm)f(r/adm) where r is the inter- 
nuclear distance and a, is the geometric mean of the 
quantum mechanical analogues of the Bohr orbits of the 
valence electrons on hydrogen and oxygen. The same 
function applies to the binding energy of a heavy hydrogen 
atom, but the new value of a, is slightly smaller, since 
a=h?/4x°ye? where yu is the reduced mass of the electron 
in the H or D atom, and a is the (quantum-mechanical) 
orbit of the atom concerned. Since (a,,£) is the same 
function of r/a,, for both H and D, the potential energy 
is a minimum at the same value of 7/am and Emin is in- 
versely proportional to dm, whilst rmin is proportional to 
dm. Thus from this cause the difference |f2—¢;| is in- 
creased by 30 calories and the nuclear distance of OD 
as compared with OH is decreased by 0.00027 Angstrom. 
The latter decrease makes only a negligible difference to 
the moments of inertia. In the calculation which follows 
we neglect also the effect of the correction to the energy, 
which would amount to about 2.5 percent in the value of 
the separation coefficient. 

5 Mecke, Baumann and Freudenberg, Zeits. f. Physik 
81, 313, 445, 465 (1933). 

6 Van Vleck and Cross, J. Chem. Phys. 1, 353, 357 (1933). 
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through a theoretical quantum-mechanical cal- 
culation of the potential function of the H,O 
molecule, and with the corresponding force 
constants calculated the three frequencies to be 
we = 3520, w, = 3560, ws = 1660 reciprocal cm (we 
are here using the subscripts employed by 
Mecke). The nuclear separation between H and 
O came out to 1.00 Angstrom, and the valence 
angle to 100°. Using the same molecular con- 
stants, we have repeated the calculation of the 
frequencies for all three molecules H20, HDO 
and DO. The results for HDO are w,=3560, 
w,= 2600, w;=1450, and for D.,O w,=2580, 
w,= 2590, w3s=1250. As Van Vleck and Cross 
point out, complete quantitative agreement with 
experiment is not to be expected in view of the 
approximations necessary in their quantum- 
mechanical calculation. We do, however, expect 
that the relative values of corresponding frequencies 
in H,O, HDO and D.O will be fairly accurately 
given by the calculation, in spite of the differ- 
ences of 5 and 3 percent, respectively, between 
the valence angle and the nuclear distance cal- 
culated by Van Vleck and Cross, and those found 
by Mecke, Baumann and Freudenberg. Accord- 
ingly we have taken the calculated ratios of w, 
for H30: w, for HDO:w, for DO (similarly 
for the w,’s and the w,’s) and multiplied the 
frequencies found by Mecke and his collaborators 
by the appropriate ratios. In this way we arrive 
at the consistent set of frequencies set out in 
Table I. These frequencies compare in a very 
reasonable way with the Raman bands in liquid 
heavy water observed by R. W. Wood,’ which 
are included in Table I in square brackets. The 
Raman frequencies are naturally all somewhat 
lower in the liquid than in the vapor. 

Mecke and his collaborators show that the 
moments of inertia of the molecule change 
slightly with the vibrational level, and in the 
ground state are not exactly what they would be 
for a rigid structure, because of the half-quanta 
of vibration. The difference is small, and we con- 
sider that a set of consistent moments of inertia 
of quite sufficient accuracy for our purpose can 
be obtained by taking the valence angle of 
105° 6’ and the nuclear distance of 0.970 Ang- 
strom found by Mecke for the ground state, and 


7R. W. Wood, Science 79, 35 (1934). 
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Moments of inertia 





x 10 
Frequencies in cm A B Cc 
H.0 3795 3670 =1615 1.027 1.972 2.999 


[3445] 


HDO 3838 2680 1411 
[3500] [2623] 


D,O 2782 ° 2670 1216 
[2517] 


1.219 3.142 4.361 


1.846 3.938 5.784 








calculating the three moments of inertia for each 
of the three isotopic water molecules, using 
H=1.00777,  D=2.01357, O=16. For com- 
parison, the values (multiplied by 10) given by 
Mecke for the ground state of HO are: A =0.995, 
B=1.908, C=2.980. The axis of least moment 
of inertia in HDO makes an angle of 20° 38’ with 
the line joining the H and D nuclei. 

With these data Eqs. (4) fand (6)fgive for the 
equilibrium constants K,=126.9 at 25° and 
49.60 at 100°, K,=3.26 at 25° and 3.40 at 100°. 

With regard to the reliability of these equi- 
librium constants, we think that since the 
molecules have practically the same potential 
function, the case is a very favorable one; the 
accuracy of the K’s turns on the accuracy of the 
vibration frequencies; our opinion is that the 
equilibrium constants should be accurate to 
within +10 percent. 

A knowledge of K. and K, enables us to cal- 
culate the ratio N3;°/N,° of free atoms of H and 
D in equilibrium with water vapor of any isotopic 
composition. We require, however, an expression 
for N;°/N,° as a function of the (analytical) 
composition of the liquid water. 

For the equilibrium in the liquid phase we 
have (assuming Henry’s law) 


Ky! = No?/NiN2=([p1°h2°/(po°)* ]po?/Pipe- (8) 


The p°’s denote vapor pressures of the pure 
species, and the ’s partial vapor pressures. The 
vapor pressure of pure D.O is known’ over the 
range 20° to 110°. A direct measurement of po 
is of course impossible because of the rapidity 
with which the equilibrium (1b) is established 
through. ionization. In Section 4 reasons aré 





8 Lewis and Macdonald, J. A. C. S. 55, 3057 (1933). 
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given for the assumption we now make, that the 
vapor pressure for the pure species HDO is the 
geometric mean of the vapor pressure of H,O 
and D.O. This makes the quantity in square 
brackets in (8) equal to unity, and therefore 
makes K,’ = Ky. 

The relationships which we now write down 
are derived simply by combining together the 
expression for K, and K, in (1a) and (1b), the 
assumption of Henry’s law, and the assumption 
that K,’=K>». The use of concentrations and 
pressures instead of activities is evidently allow- 
able with such similar molecules. Writing R, 
for the ratio of the (atomic) amounts of H to D 
in the liquid, and R, for the corresponding ratio 
in the vapor: 


R,=(2G1+Go)/(2G2+Go) 

=(2Nipi°+ Nopo?)/(2N2p2°+Nopo’). (9) 
G;, Go, and Gz are the mol fractions of H,O, 
HDO, and D.O in the vapor phase, and the N’s 
are the mol fractions in the liquid. 
R,=[B(po°— pr°) + pr? } 

/LB(po° — p2°)+p2°/R.]J, (10a) 
where 
— K,’ 1+1/R-) 


4—-K, 2 | 


K,’ 1+1/R,. 2 K,’ 1 } } 

(Sey aoa 
4—-K, 2 4—K,' R.J 

positive root when K,’ <4, negative root when 

K,’>4. When K,’ =4 (as it would be if the two 


kinds of hydrogen atoms were distributed at 
random among the oxygen atoms) then 


B=1/(1+R.). 


For the required ratio of m3°/n,° in the gas 
phase we find 








(10b) 


(10c) 


n3° 2+(K,’)'A 
oetil I caumcmanappnn (11a) 
n° 2A +(K,’)! 


in which R, is explicitly given by Eqs. (10a, b, c) 
and 


A=[(p:°(1—8))/(62°(1/Re—B)) },  (11b) 


where 6 is given by (10b) or (10c). 
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Since the ionization potential of D exceeds 
that of H by 84 calories per gram-atom, the 
equilibrium ratio of H+ to D* in free space is 
equal to 13°/n,° exp (84/RT). 

We recall that the electrolytic separation coef- 
ficient for any process subject to two limitations, 
namely that the activated state consists of 
hydrogen atoms or ions without zero-point 
energy, and that quantum-mechanical leakage 
through the potential barrier defining the energy 
of activation is negligible, is given by: 


a°=d log H/d log D=(1/R,)dH/dD 


= (1.414/R,.)n3°/n4°. (12) 
The symbols H and .D denote the content (in 
g atoms) of the isotopes in the liquid undergoing 
electrolysis, and the superscripts ° indicate the 
limitations just stated under which a® is given 
by this expression. The factor 1.414 enters 
because the relative rates at which particles of 
H and D pass out of the activated state are 
inversely as the square root of their masses. The 
ratio n;°/n,4° now refers to free atoms or free 
ions according to which is postulated for the 
activated state. 

Examination of Eq. (11a) shows that in 
general a° varies with the ratio R, of H to D. 
The limiting values of a° in very dilute heavy 
water on the one hand, and nearly pure heavy 
water on the other, are 


1.414K.3(p1°/pe°)'(Ky’/2) 
1.414K ,}(p1°/p2°)3(2/K,’) 


Table II exhibits the dependence of a° upon the 
isotopic composition of the liquid. The figures 
taken in order along each row show in outline 
how a° changes as the proportion of the heavy 
isotope increases from 0 to 100 percent, for the 
values of K,’ indicated on the left. The calcula- 
tion has been made for the activated state con- 
sisting of free ions; for free atoms the figures are 
14 percent less at 25° and 8.4 percent less at 100°. 

It will be noticed that when K,’=4 the effi- 
ciency of the separation is the same over the 
whole range of composition. Also for electrolysis 
of water in the neighborhood of 50 percent D, 
the efficiency is independent of the value of K,’, 
and has the same numerical value as it has for 
all compositions when K,’ = 4. 


(As R.>~), 
(As R,-0). 




















B. 
TABLE II. 
K/ %D= 0 10 50 90 100 
Temperature 25° 
2 a°= 13.77 16.23 19.46 25.20 27.52 
3.3 ao = 17.68 18.66 19.46 20.99 21.42 
+ a= 19.46 19.46 19.46 19.46 19.46 
Eb a°= 21.74 20.65 19.46 17.51 17.41 
Temperature 100° 
2 aos 8.03 9.81 11.07 14.37 15.65 
ao ao = 10.06 10.61 11.07 12.06 12.17 
4 a? = 11.07 11.07 11.07 11.07 11.07 
5 ao= 12.37 11.74 11.07 10.01 9.90 








The figures set out in Table II have a more 
general significance than appears at first sight. 
As they stand, they refer to any separation 
process in which the activated state consists 
essentially of free ions (or, subject to the factor 
exp (84/RT), to free atoms also). In reality, the 
activated state is never simply this: the hydrogen 
particles in the activated state are those just at 
the top of some potential barrier characteristic 
of the particular slow process which is controlling 
the reaction rate, and depending for its form and 
magnitude upon the nature of the cathode 
surface, and possibly upon the applied cathodic 
overvoltage. We therefore consider a direct 
equilibrium between the particles in free space 
and the particles in the actual activated state, 
in which they are present at some other relative 
concentration 3/n, differing from 3°/n4°. We 
write 

n3/ng=K .n3°/ng, (13) 


where K,. is an equilibrium constant analogous 
to a distribution coefficient between two phases. 
To calculate the actual electrolytic separation 
coefficient a, we have only to multiply a° by 
the corresponding value of K.. 

The factors determining K, are discussed in 
Section 5. The energy of a given particle in the 
activated state at the cathode surface or in the 
Helmholtz double layer will not be appreciably 
affected by the substitution of D for H in the 
solution. Thus we can look upon a as being the 
product of two independent factors, one a 
property of the activated state only, the other 
a property of the ground state only (i.e., the 
interior of the solution). 

It follows that for all these mechanisms the 
electrolytic separation coefficient will vary with 
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the composition of the solution quantitatively in 
the same way. 

Our calculation of K,’ gave a value 3.3 at 25°. 
Reference to Table II shows that this requires 
that a shall decrease slightly as the composition 
of the water varies from 0 to 100 percent D; 
corresponding to the estimated +10 percent 
uncertainty in K,’, we conclude that a should 
decrease by a fractional amount lying between 
10 and 30 percent of its value in very dilute 
heavy water. This is a theoretical conclusion 
which has not yet been properly tested by 
experiment. The available evidence shows that 
(in the case of electrolysis of dilute alkali hy- 
droxide solutions with nickel cathodes at about 
10°) there is not much change in a@ as the per- 
centage of heavy isotope increases; the data are 
consistent with the theoretical prediction. 


Atomic 
percent D a 
Bell and Wolfenden® 0.2 4-6 
This paper 7 5.2-5.8 
Selwood”” 97 5.5-6.5 


The theory of the variation of a with the 
nature and concentration of electrolyte can be 
treated in a manner similar to the discussion in 
this section, but it is much more complex. It is 
an interesting point that measurements of the 
electrolyte separation coefficient can be utilized 
to investigate the thermodynamics properties of 
certain solutions, at least in principle. 


SECTION 3. EXPERIMENTAL 


The measurements reported in this section 
were done with the practical aim of finding the 
most efficient conditions for concentrating the 
heavy isotope. An accuracy of about 5 percent 
in the electrolytic separation coefficient was 
obtained. 

The form of apparatus found most satisfactory 
is illustrated in Fig. 1. Samples of electrolyte 
(about 75 cc) are weighed into the electrolysis 
vessel, and occupy the wider part A; this is 
cooled (or heated) by surrounding it up to the 
level B with a suitable bath. The condenser ¢ 
is needed only when the electrolysis is done 
above room temperature. 


9 Bell and Wolfenden, Nature 133, 25 (1934). 
10 Selwood, personal communication. 
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Acetone — 


CO, bath 
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Fic. 1. 


About 10 cc of water is electrolyzed and the 
electrolytic gas burns in a minute flame at the 
jet F. The water formed is condensed in the 
small condenser D and collected in a weighed 
tube R, which is closed with a tin-foil cover. 






ISOTOPES 223 


The electrolytic separation coefficient is cal- 
culated from the relationship: 


a= (log Wy—log Wy’) /(log Wp—log Wy’), 


where Wy, Wp, and Wy’, Wp’ are the weights 
of H and D in the solution before and after a 
portion of the water has been electrolyzed away. 

The ground joint J is not greased but is sealed 
externally with a thin ring of picein wax. The 
stream of electrolytic gas passes out of the elec- 
trolysis vessel at the temperature of the coldest 
part of the vessel (about 5° in the majority of the 
experiments) and loses the water vapor it carries 
in passing through the vertical trap surrounded 
by acetone and solid CO,. The dried gas passes 
through a spray trap T packed with absorbent 
cotton, and then through a 10 cc bulb filled with 
fine sand, to the jet F inside the condenser D. 
This amount of sand is completely effective in 
extinguishing the detonation wave which fires 
back from the jet whenever the gas flow speed 
is reduced below a critical limit. This simple 
device has proved completely satisfactory in 
eliminating all risk of explosion originating at 
the jet. The jets are made simply by partially 
fusing the end of a piece of thick-walled Pyrex 



































TABLE III. 
Apparent Apparent 
current Atomic % current Atomic % 
Cathode density Electrolyte of D a Cathode density Electrolyte of D a 
Pb 1 (a) KOH iN 7 7.4 Fe 1 (s) KOH iN 7 7.6 
Pb 1 (a) . = 8 7.2 Fe 1 (s) _"s 8 69 
Pb 1.3 (s) H.SO, “ 8 6.2 
Pb 1.3 (s)  @ 9 6.6 Ni 1 (s) " «4 7) 655 
Pb 1.5 (b) 2 «® 12 6.3 
Cu 1 (s) KOH 3N 7 668 
Pt 1 (s) KOH “ 7 17.6 Cu 1 (s) H.SO, 4N 7 5.5 
Pt 1 (s) _ = 8 6.5 Cu 1 (s) “ — 43N 8 5.8 
Pt 1 (s) H.SO, “ 8 5.7 
Pt 1 (s) _— 9 5.7 Ag 1 (s) KOH 3N 7 5.3 
Pt 14 (b) Y - 8 5.7 Ag 1 (s) =. 3 8 5.8 
Pt 14 (b) “  3N 12 64 Ag (c) 1 (s) H.SO, “ 8 6.0 
Pt 14 (b) “ —1N 18 }§6 6.3 Ag (c) 1 (s) - « 10 5.6 
Pt (black) 2 “ —3N 12 4.6 
Pt (black) 2 “« 3N 15 4.7 Hg 0.7 _ * $ 2 
Pt (black) 2 “ 1N 21 4.7 Hg 0.7 “ 2 9 29 
Pt (activated) 4 “ 3N 7 4.7 Hg 0.7 “ 3N 10 2.7 
t (activated) 4 “« 3N 9 47 Hg 0.7 “ —1N i5 28 
Pt (activated) 32 “« 4N 8 3.4 
Pt (activated) 32 “ —1N 11 3.6 Ga (supercooled) 1.5 “ —4N 8 4.1 
: Ga (supercooled) 15 <5 ” 10 4.6 
W 6 (b) “ 43N 12 53 








(a) metal disintegrates; (b) smooth surface; (s) surface roughened with sandpaper; (c) metal blackened during use. 





























224 


capillary so that the diameter of the orifice is 
about 1/15 mm. 

The condensed water is always slightly acid 
with nitrogen oxides formed at the outer surface 
of the flame. The water produced was weighed, 
then made alkaline with a few mg of anhydrous 
sodium carbonate, and a few crystals of potas- 
sium permanganate added. The solution is boiled 
for a short time to free it from CO, and then 
distilled, avoiding bubble-formation as far as 
possible. The specific gravity was determined in 
a 6 cc pycnometer to an accuracy of 2 parts ina 
hundred thousand. The composition of the 
water was calculated on the basis of the figure 
1.1056 for pure DO (ds) assuming a linear 
relationship between composition and _ specific 
gravity. 

In making up the electrolyte solutions the 
heavy water mixture was twice distilled from 
slightly alkaline permanganate before use. Or- 
dinary C.P. potassium hydroxide and sulphuric 
acid were used, the acid having first been sub- 
jected to. electrolysis for some time to remove 
metallic impurities. 

In some of the experiments the water vapor 
was not frozen out of the gas stream; in these 
cases the whole of the vapor carried by the gas 
is of course condensed out along with the water 
formed in the flame, since the system operates 
at constant volume, being in fact normally closed 
by a drop of water hanging at the exit of the 
tube from the condenser D. A correction is 
therefore very easily applied to take account of 
this source of slight error. 

The experimental results are recorded in a 
condensed form in Table III. The temperature 
of the electrolyte was between 10° and 26° in all 
cases, and it is sufficiently close to compare the 
results on the supposition that the temperature 
was in the neighborhood of 20°. 

The second column gives the apparent current 
density in amperes per sq. cm when no allowance 
is made for the actual condition of the surface. 
Some idea of the relative current densities is 
obtained by taking into account the notes 
referred to by letters in brackets following the 
figures. The composition of the solution with 
respect to electrolyte is given in the third 
column, and with respect to the isotopic content 
in the fourth. These are average values during 
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the electrolysis of the 10 cc samples. The last 
column gives the values of the electrolytic 
separation coefficient a. 

The main conclusion is that there are very 
considerable differences between different metals, 
the order of efficiency taking a general view of 
the figures both for acid and alkaline solutions 
being: Pt and Pb; Fe; Cu and Ag; Ni; W; liquid 
Ga; Hg. Platinum black and “‘activated”’ plati- 
num are omitted from this series for the reason 
discussed below. 

The results for bright platinum show that the 
current density can be varied considerably 
without much affecting the efficiency. This has 
also been commented on in the case of nickel 
cathodes and alkali hydroxide solutions by 
Urey," and by Bell and Wolfenden.® 

The results for mercury and for bright plat- 
inum show that the electrolyte concentration and 
the proportion of heavy isotope can be doubled 
without noticeable effect on the efficiency. In 
general the value of a@ is about 15 to 20 percent 
lower for the same metal in the acid than in the 
alkaline solution; silver is only an apparent ex- 
ception, since during the electrolysis in sulphuric 
acid the metal is blackened, presumably because 
of slight cathodic reduction to sulphide, so that 
the surface of the cathode is no longer simply 
silver. 

The “platinum black” electrodes were pre- 
pared in the usual way by deposition from a 
platinic chloride solution, without addition of 
any lead acetate. The electrodes referred to as 
“activated” were also platinum electrodes cov- 
ered by a grey-black deposit of platinum, but the 
activation of the surface was the result of elec- 
trolyzing for some time with alternating current 
of 60 cycles. 

Electrolysis of water by 60 cycle alternating 
current has been studied by Shipley and 
Goodeve,'? who found that a large proportion of 
the electrical energy could be used to effect gas 
evolution. Our experience was, that although 
this is so for a short time, the electrodes soon 
became so strongly activated that the yield of 
gas fell to practically nothing. This applies both 


"H.C. Urey, Science 78, 566 (1933). 
12 Shipley and Goodeve, Trans. Am. Electrochem. So. 
52, 375 (1927). 
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SEPARATION OF THE HYDROGEN 


to platinum and to silver electrodes; the latter, 
as Shipley and Goodeve found, are not seriously 
corroded; nevertheless silver becomes in time 
extremely active in causing the recombination of 
the electrolysis products. 


The local concentration of heavy water at the 
cathode 


During electrolysis the liquid round the 
cathode tends to increase in content of D relative 
to the main bulk of the solution, since the mobil- 
ities of the different kinds of hydrogen ion are not 
very different,’* but the light isotope is preferen- 
tially discharged. This tendency is offset in 
several ways: by stirring of the solution by the 
gas evolved, by convection caused by temper- 
ature inequality, and by diffusion. With current 
densities such as are used in practice there is 
certainly vigorous stirring of the bulk of the 
liquid by the stream of bubbles leaving the 
electrodes. It may be, however, that whilst 
discharge of ions takes place all over the surface, 
the actual growth and detachment of bubbles 
tends to occur mainly at certain favored places, 
the hydrogen atoms or molecules moving over 
the metal surface to these places. In that case 
there would be a very thin layer of liquid covering 
parts of the cathode surface relatively remote 
from the places of bubble formation, which would 
not be effectively stirred. In such a stagnant 
layer mixing would take place mainly by con- 
vection and diffusion. As regards convection, 
there is a liberation of about one-tenth of a 
calorie per second per volt overvoltage for a 
current density of 1 ampere per sq. cm, but the 
greater part of this will be taken up by the 
metal in the first place. 

As regards diffusion, an estimate of the effect 
to be expected can be made in the following way: 
we assume that covering the greater part of the 
cathode surface there is a thin layer of liquid 
which does not suffer mixing into the main bulk 
of the solution either by mechanical means or by 
thermal convection. Outside this layer there is 
complete mixing by the stream of bubbles. 

A stationary state will then be set up in which 
a definite composition is maintained at the 
boundary of the layer, and a definite diffusion 





"Lewis and Doody, J. A. C. S. 55, 3504 (1933). 
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gradient is established across the layer. The con- 
ditions for such a stationary state in a solution 
of an acid are that the whole of the current across 
the boundary must be carried by hydrogen ions, and 
that the rate of neutralization of ions of each 
isotope is equal to the rate of migration of the 
isotope into the layer minus the rate of diffusion 
of the isotope out of the layer. 

It is sufficient to consider the special case of a 
dilute solution of a strong acid, in water con- 
taining a small proportion of D. 

Let Cy, Cp, and Cy’, Cp’ define the isotopic 
composition of the liquid, expressed in g atoms 
of H and D per cc in the bulk of the solution, and 
at the surface of the metal, respectively. Prac- 
tically the whole of the D isotope is in the form 
of HDO molecules, and we assume for the 
purpose of obtaining a rough estimate, that the 
relative numbers of H and D ions at any place 
are in the ratio of the total H and D concen- 
trations at that place. We neglect the transport 
of water arising from hydration of the ions. For 
the stationary state we have 


K I Cy’ Cp 
= (Co Co’) ~—( _ ). (14) 
Xo F Cy’ +aCyp’ Cy+Cu 





where K =diffusion coefficient of HDO in H.O, in 
cm? sec.—'; J=current density in amperes per sq. 
cm; F=the Faraday, a=the electrolytic separa- 
tion coefficient; x9= thickness of stagnant layer, 
in cm. It is reasonable to ignore the effect of the 
difference between Cp’ and Cp provided this 
difference does not amount to more than say 2 
percent. We use Eq. (14) to estimate in a typical 
case what thickness x» of the unstirred layer 
would be required to produce a 2 percent excess 
of Cp’ over Cp. We take: 


Cp/(Cp+Cy) =0.05; Cy’/Cp= 1.02; I=1; 


K=2.5X10-. 


a=7; 


This gives x»~4X10-* cm. It is not obvious 
than an unstirred layer of this thickness may not 
cover part of the surface, and it is therefore 
possible that the experimental values of a are a 
little lower through this cause than the correct 
value, which would be found at much lower 
current densities. 
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According to Eq. (14) the value of xo re- 
mains of the same order of magnitude as a 
increases from about 2 to infinity. 

The experimental results, indicate that the 
effect of local concentration round the cathode is 
small: thus when the current density at a smooth 
platinum electrode is increased 14-fold, there is 
no decrease in a whereas Eq. (14) indicates a 
decrease depending upon the value of xo. We 
conclude that with smooth metal cathodes the 
effect of any local increase in concentration in 
heavy water can be neglected in determinations 
of a which on other grounds do not claim an 
accuracy greater than about 5 percent. 

When the cathode surface possesses a very 
irregular structure, as for example when a metal 
is coated with platinum black, the situation is 
different. The tendency will be for ions to be 
discharged at the outer points, but insofar as any 
electrolysis takes place within the irregular 
shaped cavities, stirring, convection and diffusion 
will all be much less effective. This would account 
for the low efficiency of the separation on plati- 
nized platinum cathodes, and the explanation is 
supported by the fact that increase of current 
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density in this case did lower the efficiency still 
further. 


Section 4. THE VAPOR PRESSURE OF LIQUIDs 
CONTAINING IsotoPIC MOLECULES 


In Section 2 the assumption was made that 
the vapor pressure of pure HDO (0°) would be 
the geometric mean of the vapor pressures of 
H.0 (p:°) and D.O (p2°). It happens that in this 
particular case examination of the numerical 
magnitudes involved shows that so far as the 
requirements of the argument in Section 2 are 
concerned, the situation is much the same if one 
merely assumes that fo° lies somewhere near the 
mean of p:° and p2°, because in the range 20° to 
110° p2°/p.° is of the order 0.9. 

It is, however, desirable to discuss the relative 
vapor pressures of liquids that differ only in 
having different isotopic atoms in their molecules. 
We formulate the discussion in terms of the par- 
ticular example water, applying the conclusions 
reached to the question of the vapor pressure of 
the species HDO. 

The partition function for an isolated H,0 
molecule in the gas phase is 


1 82°(87°A 1B,C1)*(RkT)! 4 





fi=ex 
; - kv h 


where x;=the energy in its lowest state of an 


accessible volume per molecule in the gas phase. 
The other symbols have the meanings defined 


nk, ki, ly J —exp (—hv;/kT) 2h* 


H.0 molecule, the zero of energy being the same 
molecule in its lowest state in the liquid. Vi=the 


for Eq. (2). Ideally we should solve the mechan- 





1 


, (15 





ical problem for the whole mass of water and 
use the normal modes of quasi-oscillation and 
with these the appropriate expressions 1! 
—exp(—hv;/kT) to replace the translational part 
of (15). Instead of the complete treatment, which 
is not possible, we write for the partition function 
of a molecule in liquid water: 


1 872°(87°A ;'B,’C;’)} 








oscillation of the molecule as a whole. 


H.0 molecule in the liquid. 


where the v’s with the subscripts m1, 1, 01, 
denote the internal vibration frequencies of the 
H.0 molecule in liquid water, and the subscripts 
11, 51, t;, denote the three frequencies of quasi- 


A’, By’, C;’ are the moments of inertia of the 


(RT)3-4, (16 


=m, m1» 1 1 —exp (—hy;/RT) 71 14 1 —exp (—hv;/kT) 2h’ 





In g; the three frequencies with the subscript: 
ri, Si, ty are identified as the quasi-oscillatory 
motion of a simple molecule in the potential field 
of the surrounding molecules, but they can b 
thought of as averaged values of the quas- 
frequencies in a large aggregate of liquid mole 
cules. The condition that the vapor shall 
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saturated is f;=g,. Assuming that the saturated 
vapor has the molar volume of a perfect gas, we 
can write V;=kT/p,° and obtain at once an 
expression for the vapor pressure of the liquid. 
Similar expressions can be written for the vapor 
pressures po° and 2° of liquid HDO and D,O, 
replacing the subscripts 1 by 0 or by 2 wherever 
it occurs in the partition functions. 

Water vapor is a very imperfect gas, so that 
we ought in the vapor pressure equation to take 
account of interactions between gas molecules as 
well as between molecules in the liquid. The use 
of a harmonic oscillator and Planck functions for 
intermolecular motion in the partition function 
of a liquid is also a very approximate procedure, 
but our conclusion which depends on the dif- 
ferences between molecules containing isotopes 
would not be invalidated by modifications in 
these details. 

A liquid is a system intermediate between a 
gas and a solid, and a range of possible simplified 
mechanical models can be imagined which 
approach more or less closely one of two extremes. 
In the case of water we utilize the measurements 
recently made by Lewis and Macdonald" of the 
vapor pressure of pure D.O to obtain at least 
a clear indication of the proper choice of idealized 
mechanical model from which to calculate the 
vapor pressure of pure HDO. 

The experimental data are reproduced by the 
equation 


po°/p:° = 1.35 exp (—259/RT) (17) 


in the range 20° to 90°. The fact that the latent 
heat of evaporation is larger for the heavier 
molecule, in spite of the intermolecular forces 
being the same, is in part to be attributed to 
the lower zero-point energy of the quasi-oscilla- 
tions of the heavier molecule. The generalization 
that as between chemically similar substances 
the one of higher molecular weight usually 
vaporizes at higher temperatures is a manifesta- 
tion of the part played by the zero-point energy 
in the latent heat of vaporization." 

We consider first the interpretation of the 
factor 1.35 in Eq. (17), and for this purpose we 


Gime 


* Lewis and Macdonald, J. A. C. S. 55, 3057 (1933). 
* Cf. Urey, Brickwedde and Murphy, Phys. Rev. 40, 1 
(1932), 
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can compare the figure 1.35 with the different 
numerical factors obtained by making different 
extreme assumptions as to the molecular me- 
chanics of the liquid. 


Case I 


The molecules rotate freely in the liquid as in 
the gas, and have the same moments of inertia. 
The frequencies of the quasi-oscillation (as well 
as the internal frequencies) are such that hy >kT. 
Then from the partition functions (15) and (16) 
we obtain for the ratio of the vapor pressures of 


DO and H.O 
p2°/pr° = (m2/m,)) exp ((— x2+x1)/RT) 
= 1.17 exp ((—x2+x1)/kT). (18) 


Case II 


This differs from Case I only in the assumption 
that for the quasi-oscillation frequencies we have 
the relationship hy<kT; whence 


1 kT 2xkTm' 


1—exp(—hv/kT) hv hK} 





where m is the mass and K the force constant of 
the quasi-oscillator considered. The assumption 
of so low a frequency is of course physically 
equivalent to assuming practically free trans- 
lation. This leads to 


p2°/pi° =exp (— x2+ x1) /kT. (19) 


Case III 


Here we assume that the rotational degrees of 
freedom in the gas are in the liquid replaced by 
oscillation for which hy >kT, and that there are 
no lateral quasi-oscillations of type assumed in 
Case I. This gives us 


p° (=) (—xe+x1) 
—={ — exp ——--—— 
pi° A,B,C, kT 








(—xetx1) 
kT 


(20) 


= 2.63 exp 


Case IV 


Finally we have the extreme case which is the 
opposite of Case II: both the rotation and the 
translation of the gas molecules are replaced in 
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the liquid by oscillation for which hy >kT. Then 


p2° (“) (==) (—xi tx) 
pi? my A,B,C, kT 


= 1.17 2.63 exp (—x2et+x1). (21) 





In the present instance of water, we attribute 
the excess of 259 calories in the latent heat of 
vaporization of D,O as compared with H:0 to 
the net result of the operation in opposite direc- 
tions of two effects: the zero-point energy dif- 
ferences of the quasi-oscillation of the H2O and 
D.O molecules as a whole, and the difference 
between these molecules in respect to the de- 
crease in the internal vibrational zero-point 
energies associated with the fact that the Raman 
frequencies in the liquid water’ are lower than 
in the vapor; the effect of the former outweighs 
the latter, so that the zero-point energies of the 
quasi-oscillators must be considerably greater 
than 259 calories. Case II is therefore ruled out. 

Comparing the empirical constant 1.35 with 
Case I, III and IV we conclude that the actual 
condition of liquid water approximates more 
closely to Case I, i.e., there is nearly free rotation 
of the water molecules. 

Returning now to the vapor pressure of the 
species HDO we remark that Cases I, III and IV 
all predict for HDO a numerical factor multi- 
plying the exponential term in the expression for 
the vapor pressure ratio po°/p1° which is less than 
4 percent different from the geometric mean of 
the corresponding factors for D2O and for H,O 
(unity for the latter). 

Further, if we consider the fact that the quasi- 
oscillation frequencies depend on numbers as 
close as the reciprocal square roots of 18, 19 and 
20, and if we examine the differences between 
the Raman frequencies and the infrared absorp- 
tion frequencies given in Table II, we see that 
the latent heat of evaporation of HDO would not 
be far from the arithmetic mean of the other two. 
The arithmetic mean appears as an exponent and 
corresponds again to the geometric mean rela- 
tionship for the vapor pressure of HDO. We 
conclude that for the different kinds of water, the 
geometric mean relationship cannot be a bad 
assumption. 

From the data of Taylor and Jungers'® it 


6 Taylor and Jungers, J. A. C. S. 55, 5057 (1933). 
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appears that the vapor pressure of liquid ND; 
is related to NH; by the expression: pyp,/Pnu, 
=1.315e1%/"", The moments of inertia!’ of 
NH; are 2.79, 2.79, 4.37 and of ND; 5.31, 5.31, 
8.74. Reasoning as above, we have for com- 
parison with the experimental factor of 1.315, in 
Case I 1.275, in Case III 2.69. Thus the molecules 
rotate nearly but not quite as completely in the 
liquid as in the vapor. 

Just after the completion of this paper, Lewis 
and Schutz'® report that the vapor pressure of 
deutacetic acid CH;COOD is not lower but 
higher than that of ordinary acetic acid. It seems 
to us that the mechanical condition of acetic 
acid may well approach more closely to Case IV 
than to Case I, since the rotation of the double 
and single molecules in the vapor phase may be 
considerably hindered in the liquid phase by 
association to larger aggregates, and we note 
that if water showed the behavior corresponding 
to Case IV the vapor pressure of D2O would be 
2.2 atmospheres at the boiling point of H,0. 
The mass difference between heavy and light 
acetic acid is relatively much less than with 
water or ammonia, so that the half-quanta of the 
quasi-oscillations are less important in the liquid 
acetic acids. Of course the moments of inertia 
also differ relatively less with the acetic acids 
but the half-quanta unlike the moments of 
inertia enter into the partition functions expon- 
entially. The problem offered by a largely 
associated vapor is naturally more complex than 
when the association can be neglected. 


SecTION 5. THE ACTIVATED STATE FOR 
ELECTROLYSIS 


In Part II of this paper we discuss the details 
of three mechanisms which have been proposed 
for the process occurring at the cathode during 
the electrolysis of water, in the light of the avail- 
able data on the isotopic separation, and from 
the standpoint of the theory of hydrogen over- 
voltage. In this section we confine our discussion 
to the reason why the values of the separation 
coefficient on different metals are as low as 
experiment shows them to be. 

The highest electrolytic separation coefficient 
so far found is only a little over one-third of the 


17 W. S. Benedict, personal communication. 
18 Lewis and Schutz, J. A. C. S. 56, 493 (1934). 
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maximum value of approximately 19 calculated 
on the assumption that penetration through the 
barrier is negligible. This is to be expected. The 
activated state whatever its nature is not at a 
very much higher level than the initial state. 
One sees this from the fact that the temperature 
coefficient of the rate of electrolysis at different 
cathode potentials corresponds to an activation 
energy of about 18,000 calories for zero over- 
voltage at a mercury cathode, and decreases as 
the overvoltage is raised. Thus the activated 
state must be one in which the hydrogen par- 
ticles are still strongly bound to other atoms. 
For example, if the activation is required for the 
passage of an H or D ion from its hydrated state 
to a state in which it is bound to the metal, then 
at the summit of the barrier it is to some extent 
bound to both oxygen and metal. Or, if the slow 
process were the desorption of hydrogen mole- 
cules the activated state would be one in which 
the hydrogen atom is bound both to another 
hydrogen atom and to the metal. 

In any such activated state the particle at the 
summit of the barrier is momentarily not chang- 
ing its potential energy as it moves forward over 
the barrier, and for this motion along the path 
no half-quantum is to be considered (instead we 
have tunnelling, which in practice is treated as a 
separate problem). Motion in a plane normal to 
the classical path of the particle will be accom- 
panied by increased potential energy, and may 
be thought of as giving rise to two half-quanta. 
If the particle at the summit of the barrier is 
partially bound to several metal atoms in the 
surface of the cathode, additional degrees of 
freedom with associated half-quanta will enter. 

Thus the constant K, introduced at the end of 
Section 2 depends on the relative sizes of these 
half-quanta in the activated state for the two 
isotopes. The larger the half-quanta the smaller 
will be the ratio 13/n4 in the activated state 
compared with m3/n, in free space, which means 
a low value for both K, and the electrolytic 
separation coefficient. In general terms one 
would naturally assume that in respect of these 
half-quanta the activated state is intermediate 
between the initial state of a hydrogen particle 
in water and the state of a free hydrogen particle 
with no half-quanta; so that the value of a must 
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lie between 1.4 and 19. The experimentally 
found range of values seems entirely reasonable 
from this point of view. There will naturally be 
differences from one metal to another, but the 
order in the series of separation coefficients is 
not related in any obvious way to any one 
property of the metals. 

We remark here that the importance of taking 
account of the existence of half-quanta in the 
activated state is quite general for the relative 
rates of reactions involving isotopes: consider a 
bimolecular reaction such as H2+I,=2HI where 
in the activated state there are six half-quanta 
to consider. It is conceivable that for certain 
reactions of this type, just as has been pointed 
out by Polanyi" for the particular case when the 
initial state involves free atoms, that the heavier 
isotope will react more rapidly. The general 
condition for this to happen is simply that the 
sum of all the half-quanta of vibration in the 
activated state exceed this sum for the initial 
state by a sufficient amount to outweigh the 
velocity factor 1.414 and any tunnelling which 
may take place. It is to be expected that examples 
of all intermediate cases will be found. 

As long ago as 1919 Lindemann” drew atten- 
tion to the fact that in principle isotopes could 
be separated by means of chemical equilibria. 
Exact calculations of several gaseous equilibria 
involving heavy hydrogen have been given by 
Urey and Rittenberg*! who show also that the 
differences in the (equilibrium) electrode poten- 
tials of the heavy and light isotopes are too small 
to produce any appreciable separation by frac- 
tional electrolysis. The advantage (for the 
separation) of using reaction rates rather than 
an equilibrium process arises entirely from the 
possibility of selecting a reaction for which the 
activated state has a lower zero-point energy than 
the final state (or of finding a reaction in which 
there is a potential barrier which allows some 
tunnelling). For a favorable reaction the lighter 
isotope is most plentiful in the first material 
formed, no matter in which direction the reaction 
proceeds. This is in marked contrast to equi- 
librium processes. 


19 Polanyi, Nature 133, 26 (1934). 
20 F, A. Lindemann, Phil. Mag. [6] 33, 173 (1919). 
21 Urey and Rittenberg, J. Chem. Phys. 1, 137 (1933). 
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In a preliminary note” on the experimental 
results, we expressed the view that if the activa- 
tion is required to loosen the binding of an Ht 
or D* ion to water, and that the potential barrier 
representing this activation lies symmetrically 
in the Helmholtz double layer, then the well- 
known results of Bowden” indicate that this 
barrier is 1.5 Angstroms in total width; the 
height of about 18,000 calories may be in effect 
diminished by the fall of potential across the 
double layer to’ something like 5000 calories for 
the high current densities employed in practice. 
On metals of low overvoltage the barrier is still 
lower. If one assumes that the symmetrical 
barrier is roughly parabolic in form, then the 
penetration of the barrier occurs near the top 
and can be calculated by a formula due to 
Wigner™ which expresses the total amount of 
non-classical penetration of the barrier in terms 
of the curvature at the top of the barrier. For a 
barrier which is approximately an inverted 
parabola with a curvature at the top equal to 
that of a parabola of height 5000 calories and 


22 Topley and Eyring, J. A. C. S. 55, 5058 (1933). 
23 Bowden, Proc. Roy. Soc. A126, 107 (1929). 
4 E. Wigner, Zeits. f. physik. Chemie B19, 203 (1932). 
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half-width at the base of 0.75 Angstrom it 
appears that the effect of the leakage would be 
to increase the rate of discharge of H by about 
20 percent and of D by about 10 percent, which 
is only of secondary importance as a contribution 
to the separation. 

In a personal communication Messrs. Ogden 
and Bawn of the Manchester laboratory inform 
us of the results they obtain by using Eckart’s 
exact solution®® of the Schrédinger equation. 
They find a very much greater leakage. This is 
directly traceable to the form of Eckart’s poten- 
tial, which has a very much greater curvature 
near the top than the surface we have considered. 
This strong separation of the isotopes is not 
found in practice, and it seems to us that the 
Eckart barrier has altogether too sharp a cur- 
vature at the top to correspond to the actual 
conditions. But in any case the interpretation of 
so small a width for the Helmholtz double layer 
is very difficult; if the barrier is in reality wider 
the curvature and hence also the contribution to 
isotope separation made by quantum-mechanical 
leakage may be even less than in our original 


assumption. 


25 C, Eckart, Phys. Rev. 35, 1303 (1930). 
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An x-ray study of the series of n-aliphatic alcohols 
CioHe1OH to CisH37OH has been made. Spacings along 
identity periods perpendicular to the 001 planes (c spacings) 
have been measured at room temperature and at —50°C. 
By plotting these c spacings against the number of carbon 
atoms in the chain, it is shown that there are two distinct 
crystalline modifications, the second appearing only in 
those members having even numbers of carbon atoms. 
The common (a) modification is shown to crystallize with 


the c axis approximately at right angles with the 001 plane 
while the other (8) has the c axis at an angle of approxi- 
mately 51° 55’ with the 001 plane. The packing at the 
ends of the molecules is different in these two modifications. 
A spectrograph to make such x-ray studies at low temper- 
atures has been designed and constructed. These studies 
are the foundation of a calculation of the intensities of 
such ¢ spacing reflections presented in a following paper. 





INTRODUCTION 


ERIES of organic compounds, particularly 

long chain aliphatic compounds, have always 
yielded results of great interest when studied by 
x-ray methods.! Such series as have been studied 
in this manner are the straight-chain hydro- 
carbons, the n-aliphatic acids and the saturated 
aliphatic ketones. Other long chain aliphatic 
compounds have been studied, for example, 
unsaturated acids, dicarboxylic acids, esters and 
soaps, but in these cases, the series have never 
included more than two or three individual 
members. The accuracy of this previous work, 
as regards the spacing along the ¢ axis, is open 
to serious criticism. From the experimental 
results of Slagle and Ott,” it is known that the 
observed c spacings do not necessarily give the 
projection of the chain length or double the 
chain length, as the case may be, of that particu- 
lar chain present in greatest quantity in the 
sample but that the observed spacing corresponds 
to an “‘average”’ length of all the various chain 





* A portion of the thesis submitted by Donald A. Wilson 
to the Board of University Studies of The Johns Hopkins 
University in partial fulfillment of the requirements for 
the degree of doctor of philosophy. 

+ Pennsylvania state fellow under the National Fellow- 
ship Plan; holder of the John Wiley and Sons’ fellowship. 

‘Compare Strukturbericht, Leipsig (1931) and S. B. 
Hendricks, Chem. Rev. 7, 431 (1930). 

* Slagle and Ott, J. A. C. S. 55, 4396, 4404 (1933). 
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lengths present. Consequently, the presence of 
neighboring members in a sample of a supposedly 
pure compound would tend to make the observed 
value inaccurate as regards the true ¢ spacing 
for that compound. This is undoubtedly the 
case with some of the values of c spacings 
reported in previous work and accounts for some 
of the anomalies which have been observed. 

In this laboratory, the series of m-aliphatic 
alcohols, CypH2,;0H to C,sH37OH, of high purity, 
were placed at our disposal through the courtesy 
of Professor Reid. 

The n-aliphatic alcohols have been studied to 
some extent by Malkin.’ With a Muller spectro- 
graph and iron Ka radiation, the c spacings were 
observed for the n-alcohols Cio, Cis, Crs, Cis, 
Cio, Co1, Coo, Cos and Ceo; by measuring the 
reflections from melted and pressed samples on 
glass plates. It is believed that the present study 
of the alcohols is justified, because of the high 
purity of the materials as well as the fact that 
the available series, although shorter than that 
at the disposal of Malkin is complete in every 
member while the other shows many absences of 
individual compounds. 

In connection with the results of Slagle and 
Ott, mentioned above, concerning known simple 
and complex mixtures of fatty acids, a similar 
study of complex mixtures of these alcohols 
should prove to be of some interest. 


3 Malkin, J. A. C. S. 52, 3739 (1930). 



































































-. &. 


EXPERIMENTAL 


The alcohols used in this work were prepared 
by J. D. Meyer,‘ working under the direction of 
Professor E. E. Reid and are of a high degree of 
purity. The melting points were determined with 
good precision from a time-temperature cooling 
curve, according to the method introduced by 
Andrews, Kohman and Johnston.5 These “setting 
point’”’ values are indicative of the purity and 
are to be found in the work mentioned above. 

The source of x-rays used throughout the 
course of this work was a metal hot filament 
type x-ray tube with aluminum windows made 
in the Seemann Laboratories in Freiburg i. Br. 
A copper anticathode was used and when oper- 
ated at ca. thirty kilovolts and twenty to thirty 
milliamperes the radiation consisted essentially 
of copper Ka radiation. A nickel filter removed 
the copper KB radiation as well as some of the 
continuous region. At no time during the course 
of the spectrographic work was any radiation 
other than copper Ka observed; in particular, 
the L series of tungsten was-absent. 

The two spectrographs used were both of the 
Bragg oscillating crystal type. The material to 
be studied was placed on a glass plate either by 
melting, crystallizing from a solvent or pressing 
the solid crystals. This prepared plate was then 
mounted upon the crystal table of the spectro- 
graph and the x-ray reflections were recorded 
photographically. Reflection occurs according to 
the Bragg formula, -”\=2d sin 6. According to 
the usual method, when the crystal plate has 
been oscillated so that sufficient length of ex- 
posure has been afforded the reflections on one 
side of the primary beam, the crystal table is 
rotated 180° and the procedure repeated on the 
other side. Thus, by measuring the distance 
between corresponding lines on either side of the 
primary beam, the glancing angle @ can be 
ascertained, knowing the essential dimensions of 
the spectrograph. Knowing @, sin @ can be found 
and since the order of reflection (obtained 
definitely by observation of several orders) and 
\ the wave-length of copper Ka radiation are 
known, the Bragg formula can be solved for d, 


* Meyer and Reid, J. A. C. S. 55, 1574 (1933). 
5 Andrews, Kohman and Johnston, J. Phys. Chem. 29, 
914 (1925). 


WILSON AND E. OTT 





the interplanar distance. The value used for the 
wave-length of copper Ka radiation is 1.538674A 
as obtained from Siegbahn® and is the value of 
the center of gravity of copper Ka; and Kaz. 
The opening of the slit was 0.1 mm throughout 
the entire work in both spectrographs. 

The value of the interplanar distance measured 
in this manner, when applied to crystal analysis, 
is not to be confused with the spacing along the 
c axis of the crystal. It is rather the projection 
of the spacing along the c axis on a line perpen- 
dicular to the glass plate. For this reason, a 
graphical comparison of these ‘‘projection”’ spac- 
ings can be used to show the presence of various 
modifications when these modifications make 
characteristic angles with the glass plate. For 
the purposes of brevity and to follow the usual 
custom, the projection of the c axis will be 
referred to merely as the c axis and the observed 
spacings will be listed as c spacings. 

Since the materials used have a melting point 
range of from 57°C to 5.9°C, it was necessary’ 
to design and construct a spectrograph which 
would keep the sample cooled sufficiently for 
the material to remain solid. However, the 
majority of the x-ray diagrams of the alcohols 
which melt above room temperatures was taken 
on a spectrograph which was made in the 
Seemann Laboratories at Frieburg i. Br. The 
calibration of this spectrograph is described in 
the work of Slagle and Ott.’ 

The design and construction of our low temper- 
ature vacuum spectrograph follows. It is es- 
sentially a Bragg oscillating crystal spectrograph 
with the crystal table located at the center of a 
circle of approximately 18 cm radius and with 
the first slit located on the circumference and 
the photographic plate located as the tangent to 
a point 180° from the slit. The exposed plate 
area is 1X11 cm which allows reflections up to 
glancing angle @=9° to be observed on either 


6 Manne Siegbahn, Spektroskopie der Réntgenstrahlen, 
Berlin, Julius Springer, 1931. 

7 The construction of the low temperature spectrograph 
was not absolutely necessary for the study of the alcohol 
series alone, but was essential for studies of other series 
with lower ranges of melting points such as the bromides 
and acetates, which are being made in this laboratory. 
The present work shows the applicability of this spectro- 
graph to such studies. 
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Fic. 1. Scale diagram of crystal table and cooling rod of the 
low temperature vacuum spectrograph. 


side of the primary beam, which is a sufficient 
range for the purposes of this and following 
investigations. The arrangement of the recording 
photographic plate as a tangent does not impair 
the focussing, since the length of the plate is 
small as compared with the circumference of the 
circle. As a collimator, a system of two slits is 
used, the first being adjustable and made of a 
hard lead alloy while the second is fixed just in 
front of the crystal table and is made of steel. 
The oscillation of the crystal is accomplished by 
an eccentric cam arrangement driven by an 
electric motor. 

The essential features of the cooling arrange- 
ment and crystal table are shown diagrammati- 
cally in Fig. 1. @ is a solid copper rod which 
extends up through the floor b of the spectro- 
graph and which is surrounded by a cooling 
mixture contained in a Dewar flask. The cooled 
portion is insulated as much as is possible with 
Bakelite parts. A lubricated steel on steel ground 
joint c affords the vacuum tight mechanism by 
means of which the oscillating motion is trans- 
mitted to the crystal table. The lubricating 
medium which is best suited for the carbon 
dioxide temperatures at which this work was 
being done is a mixture of a light aviation grade 
motor oil and a good grade gasoline in the ratio 
of 9 to 4. It is contained in the Bakelite reservoir 
d which acts in addition as a further protection 
against leakage. A brass shield e has the two-fold 





purpose of protecting the body of the spectro- 
graph against oil seepage and to keep any air 
which might leak through the ground steel joint 
in as close contact as possible to the outlet f 
which leads to an ordinary oil vacuum pump. 
The crystal table k is so constructed that a solid 
block of copper m is forced against the ledge o 
by means of a set-screw holding the crystal plate 
tightly between the block and the ledge. Here, 
the plate is cooled by the copper block m which 
is in contact with the cooling rod a. h is the 
rocking arm which leads to the cam mechanism 
not shown here. 

A mixture of CO. snow and acetone was found 
to produce a temperature of approximately 
— 50°C at the surface of the sample. On account 
of this low temperature, it was necessary to keep 
moisture from collecting on the surface of the 
crystal in some manner and this, in addition to 
reduced conduction of heat, was the reason for 
building the spectrograph with its vacuum 
feature. Under ordinary conditions of running, 
the amount of moisture which would collect on 
the crystal surface during a 16 hour period was 
barely noticeable. This, of course, is not sufficient 
to interfere seriously with ordinary x-ray re- 
flection. 

In order to know the temperature obtained in 
this manner, as well as to observe the constancy 
of this temperature with time, several cooling 
curves were plotted. The temperature was meas- 
ured with a copper-constantin thermocouple, one 
junction of which was cemented to the surface 
of the crystal plate by means of universal wax. 
The cooling curves were made under as nearly 
as possible the same conditions as were met with 
in ordinary running. These curves proved that 
by refilling the Dewar cooling vessel at 2 hour 
intervals, a temperature constant to +1° was 
kept for any desired length of time. 

It was necessary of course to calibrate the 
spectrograph, but it was impossible to do this in 
the same manner as was used to calibrate the 
Seemann spectrograph previously described, be- 
cause the range of the glancing angle @ (0°-9°) 
observable on the plate is small and too close 
to the primary beam. This means that even with 
silver radiation and the mica 001 face, an 
insufficient number of reflections would be ob- 
served for accurate calibration. This difficulty 
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was surmounted by using as standard several 
long chain organic compounds whose spacings 
had been determined with known accuracy on 
the Seemann spectrograph in this laboratory. 
Three compounds were used, the C,; ”-aliphatic 
acid, determined by Slagle and Ott? and the Ci; 
and Cis n-aliphatic alcohols determined during 
the course of this work. The combination of 
these three substances gave a series of 19 points 
on the calibration curve, plotting distance in 
millimeters between corresponding lines on the 
film against the glancing angle @. All fall accu- 
rately on a very smooth calibration curve which 
is nearly a straight line. The calibration curve is 
curved somewhat, since the photographic plate 
was placed as a tangent to the circle of focus 
rather than exactly on the circumference. The 
frame of the spectrograph is so insulated from 
the crystal table that its dimensions remain 
sufficiently constant under all running condi- 
tions. The entire work of calibration was per- 


formed at room temperature and so all of the - 


spacings determined at the low temperature are 
referred to the calibration at room temperature. 
This fact allows comparisons to be made between 
spacings observed at — 50°C and at room temper- 
ature and thus allows us to show the existence 
or non-existence of an expansion coefficient along 
the ¢ axis. 

The main body of the spectrograph was 
constructed of brass, with the exception of the 
floor which is of Bakelite and all the joints 
between brass and brass and brass and Bakelite 
were made vacuum tight with universal wax. 
The cover of the spectrograph and the brass 
plateholder were both cemented on with uni- 
versal wax and were removed with lifting screws 
after air had been allowed to enter the spectro- 
graph. The cover of the spectrograph was pro- 
vided with a red glass window directly above 
the crystal table by means of which the rotation 
of the crystal could be observed. The 180°- 
rotation is accomplished by loosening the set- 
screw g (Fig. 1) and applying a specially con- 
structed wrench to the cooling rod a. 

Although the point has not been tried, there 
is no reason to doubt that the spectrograph could 
be used to investigate the behavior of long chain 
compounds through a wide range of temperatures 
by substituting for the CQOb»-acetone cooling 
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mixture anything from liquid air to a thermo- 
static arrangement for temperatures not too high 
above room temperature. It would only be 
necessary in each case to find the correct lubri- 
cant for the rotating joint. The constancy of the 
temperature at the surface of the sample suggests 
that the spectrograph could be used to investi- 
gate behavior near melting points, the phe- 
nomenon of molecular rotation or in general, 
transition points of all kinds. 

The general method of calculating the c 
spacings from the films has been previously 
outlined. The distance / between corresponding 
lines on either side of the primary beam was 
measured by a micro-comparator in the case of 
the strong reflections and with a millimeter rule 
in the case of the weaker reflections. It was 
generally possible to measure the first, third and 
sometimes the fifth orders with the comparator. 
In each case, five separate readings were made 
of each order and the average taken as the final 
average /. Since plates were used in the low 
temperature spectrograph, it was not necessary 
to correct for shrinkage, but in the case of the 
Seemann spectrograph, which required films, 
this was necessary. The correction was made by 
making parallel scratches on the film before 
exposure, a known distance apart and then 
measuring the distance between these scratches 
just before measuring. The average / was then 
corrected, proportional to the amount of this 
shrinkage. 

Having /, the distance between corresponding 
lines, accurately measured for several orders, it 
remains to find the glancing angle @. The cali- 
bration curves for the two spectrographs were 
made by plotting / against @ and so @ is obtained 
simply by reading from the calibration curve. 
Now, knowing @ for a variety of orders, sin @ can 
be found and by knowing n the order of reflection, 
from observation, (sin @)/n is found which re- 
duces the sine values of the reflections from 
higher orders down to n=1. The variation in 
these values of (sin 6)/n then gives us the 
maximum error in the calculation of d, the 
interplanar distance, which is obtained simply 
by substituting the average (sin @)/n value in 
Bragg’s formula. Since (sin @) /n is averaged over 
several orders, the maximum deviation from the 
mean which is given in connection with the 
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Low temperature spectrograph 


Seemann spectrograph 


Cj,2H.;O0H (a) 


~ ¢,,H.OH (a) i ; CisHs7OH (8) Cj 3H JH (a) 


Order - --— ~ ——-— 

of ’ Devi- ’ Devi- : Devi- F Devi- 
reflec- sin @ ation sin ation s ation sin 6 ation 
tion n (%) (“%) (%) n (%) 





+0.044 — — 

0.02081 —0.144 
0.02083 —0.048 
0.02087 +0.144 
0.02085 +0.048 
0.02083 —(),048 


0.01825 —0.05 — — 0.02252 

0.01828 +0.075 0.01876 —0.107 — — 

0.01826 0.00 0.01876 —Q.107 — ee 

0.01828 +0.075 0.01882 +0.212 0.02264 +0.57 

0.01826 0.00 0.01878 9.09 0.02242 —0.38 

0.01826 . 0.00 0.01881 +0.10 0.02246 —().22 
0.01878 0.00 


were | 


“Io 


Oo Oo == —@= SS DO DO — me 


8 0.01827 +0.05 


Average 
sin 0 


0.01826 0.01878 


42.12+0.05 40.96 +0.10 


0.02084 
34.17 +0.32 36.91 +0.05 








* It must be remembered that the c spacings calculated in this manner are not to be confused with the regular c spacings 
along the crystallographic ¢ axes. The two are identical only when the c axis is perpendicular to the glass plate and hence 
to the 001 plane. These c spacings are a measure of the projection of the true ¢ spacings on a line perpendicular to the 
001 plane which appears to be parallel with the plane of the glass plate. 


determination of d is probably much greater 
than the average error or the most probable 
error. However, in this work, it was thought 
best to give the variation in terms of the maxi- 
mum deviation rather than the error, since the 
number of observations is relatively small. 

In Table I are placed the calculations of films 
obtained from C,;H3,OH and C,sH37OH, re- 
spectively, by using the Seemann spectrograph. 
They represent one of the least and one of the 
most favorable cases and are typical of the 
accuracy obtained and the type of calculation. 
Likewise, Table I contains similarly representa- 
tive calculations made from diagrams obtained 
on the low temperature spectrograph. These are 
CiHs,,OH and C,;H2OH, respectively. The 
table of calculations is self-explanatory. 

In order to show the sharpness of reflection 
with which we are dealing and the type of films 
obtained, Fig. 2 shows several representative 
films made on the two spectrographs. The low 
temperature films, while easily measurable, do 
not reproduce well. 

The final calculated values for the ¢ spacings 
of all these aliphatic alcohols are given in Table 
Il. In this table are also included Malkin’s 
values for comparison. As will be seen later, 


these values indicate two distinct modifications 
of the n-alcohols. Quite arbitrarily, that modifi- 
cation which gives the longer spacing is called 
the a-modification whereas the one with the 
shorter spacing is called the 8-modification. 


DISCUSSION 

If we compare c spacings as obtained from 
pressed samples with those obtained from crystal- 
lized or melted samples, we see that occasionally 
the agreement is not good. Of the two values, 
the spacing obtained from the pressed sample is 
likely to be in error. This is because the pressed 
sample is thicker, presents a more uneven sur- 
face, and as Slagle and Ott? have shown is likely 
to give spacings varying with the amount of 
pressure put upon the sample. As will be seen 
in Table II, the spacings observed in this work 
do not always agree with those published by 
Malkin’ even within the maximum experimental 
deviation. This is probably not due to an 
inherent error in the methods used in the two 
investigations (although it is believed that the 
accuracy has been increased in this laboratory) 
but can be attributed to the presence of im- 
purities in the alcohol in question. As explained 
before, the purity of the alcohols used in this 
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Fic. 2. X-ray diagrams taken on the two spectrographs used in this work. Diagrams (a), (b) and (c) were taken on the 
low temperature spectrograph while diagrams (d), (e) and (f) were taken on the Seemann spectrograph. 


(a) Ci4H2»COOH—room temperature 
(b) C 4H »COOH— — 50°C 


(c) CivH3;0H— —50°C—ea-modification 


work is unusually high and consequently the 
values given here are considered the more correct. 

As in the case of the fatty acids, it is evident 
that the alcohol molecules crystallize in a double 
arrangement of some kind, probably with the 
hydroxyl groups adjacent to each other. The 
lengths of the c spacings as well as the increase 
observed for consecutive members of the series 
are much too great to be explained in any other 
manner. This viewpoint is supported by evidence 
in the paper following calculating intensities of 
reflection from assumed atomic arrangements. 

Some of the diagrams show the presence of 
two modifications. This is clearly brought out if 
we plot the c spacings against the number of 
carbon atoms, in the usual manner. The graph, 
shown in Fig. 3, is made up of the c¢ spacings 
obtained both by Malkin and the present 
authors. It therefore serves as a comparison of 
the two sets of values. The modification giving 
values which fall on the upper line is the a- 
modification while values which fall on the lower 
line belong to the 8-modification. 

It is interesting to note that in all of this 
work, there have been only two modifications 


(d) CigH3;0H—s-modification 
(e) C,;H;,OH—a-modification 
(f) CigH»OH—a-modification 


observed as contrasted with the case of the 
fatty acids which show four and possibly five 
different modifications? in either even or odd 
series. In all cases the reflection from the a- 
modification is the more pronounced, when both 
appear on the same film. However, in the case 
of some preparations of the Ci, and Cys alcohols, 
the presence of the 6-modification only was 
observed. Reference to Fig. 3 will show that the 
B-modification has appeared only in the cases 
where the alcohol has an even number of carbon 
atoms. This is probably due to some difference 
in the ease of packing for even and odd molecules. 
That there is a difference in the packing of the 
a- and 8-modification is indicated by a quite 
different relation between intensity and order 
of reflection. 

Meyer and Reid‘ working with .this same 
series of alcohols have obtained evidence from 
cooling curves which they interpret as showing 
the existence of two modifications of both even 
and odd members of the series. Our failure to 
obtain two modifications in the case of the odd 
members may be due to one of two reasons; 
either in our method of preparing the sample 
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TABLE II. Observed c spacings in A of n-aliphatic alcohols.* 








No. of 
carbon Crystallized 
atoms’ from ether 


Results of Malkin 


Pressed Melted Pressed 





Room Temperature 


36.88 -+0.07 (a) 
39.67 £0.07 (a) 
42.12-+0.05 (a) 


44.62 +0.05 (a) 
36.97 +0.03 (8) 


47.16+0.14 (a) 


49.54+0.12 (a) 
40.96 +0.10 (8) 


36.89 +0.07 (a) 
39.38 40.46 (a) 
42.13+0.08 (a) 


44.70-0.25 (a) 
36.95 +£0.30 (8) 


46.99 +0.25 (a) 


41.25+0.09 (8) 


44.9 (a) 
37.4 (8) 


50.2 


(a) 
41.35 (8) 








No. of 
carbon 
atoms 


Crystallized 


from ether Melted Pressed 





— 50° Centigrade 


10 28.74 +£0.34 (a) 
22.41 +0.06 (8) 


11 31.28 -+0.36 (a) 
34.17 £0.32 (a) 

36.91 +0.05 (a) 

39.38+0.19 (a) 

42.16-+0.25 (a) 


44.84-+40.11 (a) 
36.84 +0.20 (8) 


46.93 +0.03 (a) 


49.57+0.12 (a) 
40.96 +0.08 (8) 








* See note below Table I regarding c spacings. 


the cooling was so rapid that the transition point 
was passed without any appreciable change and 
at the low temperature, the rate of transition is 
so slow as not to give the second modification in 
any appreciable quantity, or this second transi- 
tion point is due to a change in the rotation of 
the molecule only which would have no effect 
on our x-ray diagrams but which would have a 
marked effect upon the specific heat. Which of 
these two explanations is the correct one can 
only be decided after further work. 

From our series of values of c spacings, we can 
calculate the average increase in length as the 
number of carbon atoms increases. This value 
is 2.53A for the a-modification and 2.00A for the 


umber of Carbon Atoms 


. 3. Graph plotting c values against the number of 
carbon atoms in the hydrocarbon chain. 


B-modification. Careful work with the hydro- 
carbons has shown that the increase in ¢c spacing 
with the addition of two carbon atoms to the 
molecule is 2.54A. Since, in the modification of 
the hydrocarbons referred to, the c axis is 
perpendicular to the 001 plane, the correspon- 
dence between this value and our value would 
indicate that in the a-modification, the ¢ axis is 
perpendicular to the 001 face. Likewise, the 
value of 2.00A increase per carbon atom found 
for the 8-modification, when compared with 
2.54A shows that the $-modification has the c 
axis tilted with respect to the 001 face. It is 
evident that the packing at the ends of the 
molecule would have no influence on this pro- 
jected increase. Hence we can calculate the 
angle that the c axis makes with the 001 face in 
the case of the 8-modification. If y is this angle 
of tilt, simple geometrical considerations show 
that sin y=2.00/2.54=0.787 and y=51° 55’, 
provided that we assume the chains in the case 
of the 6-modification to be also parallel to each 
other. This would seem to be a reasonable 
assumption since the chains must be closely 
parallel to each other in the case of the a- 
modification. This is shown by the close check 
between the observed and calculated maximum 
increase along the c axis. 

The tilt of the B-modification, when calculated 
in this manner, can then be used to show that 
the type of packing at the ends of the molecules 
must be different in the two modifications. This 
view, borne out further by the fact that the 
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relationship of the intensities of reflection is 
different in the two cases, has not hitherto been 
proposed. We know from the c spacing, the 
length of the projection of the total ‘‘double 
molecule.” From the angle of tilt, we can easily 
calculate the value c would have for the £- 
modification, if the chain axis were perpendicular 
to the 001 face. This calculated c spacing should 
be the same as the c spacing for the a-modifica- 
tion if the packing were the same. The calculated 
c spacings for the 8-modification are 46.99A and 
52.04A for the Cig and Cis alcohols, respectively, 
while the corresponding c spacings for the a- 
modification which is perpendicular to the 001 
plane are 44.62A and 49.54A. Although the 
experimental error in the determination of the c 
spacings is magnified in the value of sin y, we 
believe that this difference is real and not 
apparent and that it tends to strengthen our 
view that the packing is different in these two 
modifications. 

If, in Table II we compare the c spacings 
obtained at room temperature with those ob- 
tained at — 50°C, we fail to detect any difference, 
outside the experimental error. This means that 
the coefficient of expansion of these alcohols in 
the c direction is negligible. This compares 
favorably with the results obtained by Miiller® 
in the case of the hydrocarbons. However, one 
acid was investigated in this manner and a 
change in the c spacing was observed. The Ci; 
n-aliphatic acid, which was measured by Slagle 
and Ott? and used in the calibration of the low 
temperature spectrograph, gave a value at 
— 50°C of 35.19+0.10A as compared with 35.65 
+0.05A at room temperature. This is well 


8 Miiller, Proc. Roy. Soc. Al27, 417 (1930). 
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outside the experimental error and gives a value 
for the coefficient of expansion of 0.00015A per 
A per °C. Although this shift in the ¢ spacing 
with temperature lies outside the experimental 
error, the calculated coefficient of expansion 
should only be considered significant in order of 
magnitude because the experimental error is so 
large compared with the observed shift. 
Because these alcohols are extremely simple 
as compared with other long chain compounds 
in regard to the possible number of modifications, 
they lend themselves very nicely to the study 
of mixtures and solid solutions. Some work along 
this line has been done by the authors but 
because of the fundamental nature of the 
problem, it was considered advisable to expand 
this work and publish it in a following paper. 


SUMMARY 


(1) The c spacings of the series of -alcohols 
C,9H2,0H to C,sH37,0H have been studied with 
improved accuracy. 

(2) The n-alcohols crystallize in “‘double mole- 
cules’”’ and the members with even numbers of 
carbon atoms exist in two distinct crystal 
modifications. 

(3) These two modifications are shown to 
crystallize with the c axis making angles of 90° 
and 51° 55’, respectively, with the 001 planes 
and with the packing at the ends of the molecules 
likewise different in the two cases. 

(4) A special spectrograph designed to permit 
such x-ray studies at temperatures below room 
temperature has been constructed. 

(5) The results presented here are used in a 
following paper to determine atomic arrange- 
ments from intensity calculations. 
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A Calculation of the Intensities of Reflection of a Series of n-Aliphatic Alcohols* 
DonaLp A. WiLsonft AND Emit Ott, Department of Chemistry of The Johns Hopkins University 
(Received February 9, 1934) 


Comparison between cbserved and calculated intensities of c spacings of the a-modification 
of a series of -aliphatic alcohols gives excellent values for the coordinates (projected on the 


c axis) of the constituent atoms. 





N 1925 Shearer! published a paper in which he 
attempted to explain the distribution of 
intensities of reflection with order from a number 
of long chain organic compounds. Because of 
later refinements in the theories of scattering, 
these calculations offer at best only a qualitative 
explanation of the phenomena. In the light of 
some of the observed reflections from long chain 
compounds used in this work, Shearer’s results 
are actually in error. For these reasons it would 
seem advisable to use a more modern conception 
of such scattering and recalculate the intensities 
of reflection for some of the alcohols used here. 
Suppose that instead of using a model such as 
Shearer proposed, we consider each atom along 
the chain as a diffracting center. Now if we draw 
a line parallel to the c axis of the crystal, we may 
project the centers of the atoms along the 
molecule on to this line and obtain an arrange- 
ment such as shown in Fig. 1. We may deal with 
this projected row of atomic centers as a row of 
diffracting centers and by finding the coordinates 
of each of these points along the line of pro- 
jection, we may use the regular expression for the 
structural amplitude and thus calculate the 
intensities for various orders of the 001 reflec- 
tions. This simplified arrangement is strictly 
rigid for such calculations. 
In applying this general method to the case of 
the alcohols, let.us see what dimensions or 


*A portion of the thesis submitted by Donald A. 
Wilson to the Board of University Studies of The Johns 
Hopkins University in partial fulfillment of the require- 
ments for the degree of doctor of philosophy. 

t Pennsylvania state fellow under the National Fellow- 
ship Plan; holder of the John Wiley and Sons’ fellowship. 

‘Shearer, Proc. Roy. Soc. A108, 655 (1925). 


coordinates are unknown. As was previously 


C 
shown,’ the projected 7 \_ distanceis 2.54A 
S ¢ 


* 
and hence the projected / distance is 1.27A 
Cc 


for the a-modification which has the molecule 
chains perpendicular to the 001 plane. This gives 
us all the necessary dimensions relative to heavy 
atoms within the projected molecule except the 


O 
projection of the / distance. Atomic radii 
ts 


tables give the distance between oxygen and 
carbon centers as 1.42A.° This gives the pro- 


O 
jection of the /_ distance as 1.42/1.54X1.27 
C . 
=1.17A for this linkage, assuming in this case 


also, the tetrahedral angle. Having these values, 
we see that there are but two unknown param- 
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Fic. 1, Showing method of projection of the atoms along 
the n-alcohol chains, 


2—D. A. Wilson and Emil Ott, preceding paper, p. 231. 
3M. L. Huggins, Chem. Rev. 10, 427 (1932). 
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eters, 5; and 62 (their definition is evident from 
Fig. 1). 6; which measures the shift of the entire 
molecule along this axis of projection is con- 
sidered unknown, then 42 is simply related to 6; as 
follows: 


d=2{5:+(m—1)1.274+1.17+6e}, 


where d is the interplanar distance and 7 is the 
number of carbon atoms in the molecule. 

Thus we see that in an aliphatic alcohol 
containing ” carbon atoms, the carbon atoms are 
located at +6,, +(1.27+6,), +{2(1.27)+61}, 
+ {3(1.27) +61}, --- +{(m—1)(1.27)+61} and 
the oxygen atom is located at +{(1.27)(n—1) 
+6,+1.17}. These values are obtained by finding 
the distance along the axis of projection of each 
diffracting center from the origin, placed at a 
center of symmetry of the “double molecule”’ for 
convenience. However, in the expression for the 
structural amplitude, it is necessary that these 
coordinates be expressed as fractions of d along 
the z axis and hence we must transform these 
coordinates to +u, +{1.27/d+u}, +{2(1.27)/d 
+u}, +{3(1.27)/d+u}, --- +{(m—1)(1.27)/d 
+u} for the carbon atoms and + {[(m—1)(1.27) 
+1.17 ]/d+u} for the oxygen atom when u=6,/d. 
In all these considerations, the hydrogen atoms in 
each —CH2— group are considered as being 
projected to the same point on the axis of 
projection. This is merely done for convenience, 
since as is readily seen, the calculation remains 
rigid for 001 reflections. The —OH group is 
treated as a single diffracting center while the 
end hydrogen on the hydrocarbon chain is 
neglected. This simplification can certainly intro- 
duce but a minor error. 

These coordinates are then substituted in the 
expression for the structural amplitude, 


S=2E F,e?*i" =(A —iB), 


where 
A=2)>. F,, cos 2ar2n, 


and 
B=? : F,, sin 2rrzn, 


n 


where 2, is the fractional coordinate of the nth 
atom along the z axis, 7 is the order of reflection 
and F, is the scattering factor for the mth atom. 
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However, it is most reasonable to assume that 
the two alcohol molecules are arranged sym- 
metrically about the origin, which is the center 
of the ‘‘double molecule’ and corresponds to a 
center of symmetry in this case, since the 
molecule beyond the origin is the precise image 
of the one before it. The intensity is proportional 
to the square of the absolute value of the 
structural amplitude and is therefore proportional 
to A?+ B*. However, now the term B?=0 and so 
the intensity, 


I«A?«(2 >) F, cos 27rz,)? 


and because B?=0 only cosine values need be 
considered and only positive values of the 
coordinates z, are necessary. Nothing quanti- 
tative is known about the scattering factors for 
the carbon and oxygen atoms at such small 
values of 6, so as a first approximation, this 
factor was taken to be proporticnal to the atomic 
number, or 8 for the —CH2— groups and 9 for 
the —OH group. 

The F, values for oxygen and carbon‘ in this 
region of small glancing angles actually do change 
and although the shape of the curve is not known 
accurately, we find that apparently to the first 
approximation it may be compared satisfactorily 
in our region with the curve of the function 
1/sin 6. With the rotating crystal method of 
analysis used in this work, it is necessary to 
introduce a correction factor sin 6 in calculating 
intensities since the area of crystal surface 
intercepting the x-ray beam increases as sin 6 and 
the energy striking the crystal surface must 
increase as that function. Consequently, it was 
thought that the value of this calculation would 
not be seriously impaired if these two corrections 
were allowed to cancel and only the relative 
heights of the two F, curves for oxygen and 
carbon be applied as the scattering factors. 

The intensity is then proportional to A? or to 
(A/2)?, since any constant factor does not need 
to be included in a relative calculation and is 
equal to f(@)(A/2)* where f(@) is a function of 8 
including the polarization, Lorentz and ab- 
sorption factors. The value of f(@) used is given 


4R. W. G. Wyckoff, The Structure of Crystals, Chem. 
Cat. Co., (1931). 
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by Wyckoff‘ and can be shown to be valid for this 
type of calculation. 

Since this calculation is based on the as- 
sumption that the alcohol molecules crystallize in 
pairs with the polar —OH groups adjacent it is 
perhaps necessary to justify this viewpoint. 
This idea has been held by organic chemists for 
some years as applicable to ali such cases of long 
chain molecules bearing a polar group on one end. 
X-ray researches with fatty acids and particu- 
larly in the case of lauric acid,> have further 
justified this point of view, at least with the case 
of the acids. That the alcohols present an 
analogous case can be seen by comparing the 
forces of molecular cohesion,® that of the acids 
being about 8970 cal. per mole and that of the 
alcohols being about 7250 cal. per mole. Shearer! 
has presented a more complete justification of 
this assumption as a preface to his calculations. 
At any rate, the success of the present calcula- 
tions gives evidence strongly in support of this 
primary assumption. 

We can now state our complete expression for 
the calculation of intensities. 


I=f(0)(A/2)’, 


1+cos? 20 
I=————_ [ F,, cos 22, f. 
sin 20 * 


(1) 


This expression can be criticized only from the 
point of view that the Darwin “extinction 
coefficient’’ of absorption has been neglected. 
The extinction coefficient does not become 
appreciable until we deal with reflections which 
lie exceedingly close to the primary beam, such as 
the first, second and possibly. the third orders of 
reflection as observed in this work. This fact is 
used later to explain slight disagreements be- 
tween observed and calculated intensities in the 
first two or three orders. 

The C,; alcohol was chosen first for calculation 
for two reasons; first because reflections from this 
alcohol have been observed to higher orders than 
for any of the others and secondly because the c 
spacing obtained for this member is “‘average”’ in 
that the value falls accurately on the curve, 


° Strukturbericht, Leipzig, (1931). 
°K. H. Meyer and H. Mark, Der Aufbau der Hoch- 
bolymeren Organischen Naturstoff, Leipzig, (1930). 
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Parameter (ud) 


Fic. 2. Calculated intensities plotted against parameter 
values for the Cy, alcohol. 


plotting c spacings against the number of carbon 
atoms.” (See Fig. 2.) 

The coordinates of each carbon and oxygen 
atom in the chain are found in the manner 
presented. The —CH2— groups can be placed at 
u, (0.02692+4x), 
(0.10771+4), 
(0.18850+- 1), 
(0.26929+-1), 


(0.35008+1), 


(0.05385 +4), 
(0.13464+4), 
(0.21543+4), 
(0.29622+4), 
(0.37701+4), 


(0.08078+4), 
(0.16157+4), 
(0.24237+4), 
(0.32315+4), 
(0.40394+-), 


(0.43087+-) and the —OH group at (0.45568 
+u). Thus u=6,/d is the only unknown param- 
eter. It was determined by making preliminary 
calculations of the intensities, by using values of u 
increasing regularly on each side of the origin. 
Having obtained in this way a series of calculated 
intensities varying with the parameter u, the 
next step was to plot these calculated intensities 
against the parameter u for each order. That 
value of « which was found to give the proper 
sequence of relative intensities for the various 
orders, as compared with the experimental data, 
was found to be a singular value at «= +0.027. 
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TABLE I. Calculated and observed intensities.* 












Order 









“ CisH3;70H Ci7H;;0H CisH330H C15H;3,0H Ci4Hx»OH Ci3H270H 
reflec- 
tion Calc. Obs. Calc. Obs. Calc. Obs. Cale. Calc. Obs. Calc. Obs. 



















































1 1614 500 1241 625 1180 320 1020 300 849 460 1161 385 
2 312 250 349 200 354 165 371 150 419 230 251 200 
3 526 500 385 500 355 320 300 300 246 460 341 350 
4 155 150 224 175 225 160 183 120 246 200 150 150 
5 255 250 200 225 161 160 81 250 127 230 176 190 
6 155 50 146 200 121 115 66 90 137 140 120 150 
7 124 125 99 100 104 60 48 90 98 70 
8 150 175 168 165 161 120 181 185 114 115 
9 72 50 54 60 
10 174 175 

















With this value, the calculated intensities are 
found to be the values shown in Table I, where 
they may be compared with the experimental 
values. 

The calculation is made simply by performing 
the designated operations in Eq. (1) for the 











P 1+ cos? 20 
™ sin 26 


It is necessary to solve this equation for each 
order since 7 is the order of reflection and @, the 
glancing angle, depends upon the order of 
reflection. 

In order to compare these with the observed 
intensities, in Fig. 3 are shown diagrammatic 
representations of the calculated and observed 
intensities, with the absolute scale of each 
changed so that they are made to coincide inso- 
far as possible. One will note that the agreement 
is very good indeed, particularly in the higher 
orders. It would be surprising if the first two or 
three orders were to coincide exactly because we 
have, as mentioned previously, neglected the 
“extinction coefficient” of absorption which 
becomes increasingly important in the low orders 
lying very close to the primary beam. 


+cos 2rr 0.13471+cos 2rr 0.16164-+ cos 227 0.18857 -+cos 227 0.21550+ cos 227 0.24243 
+cos 2rr 0.26973-+ cos 2rr 0.29629+ cos 2rr 0.32322+ cos 277 0.35051-+c0s 277 0.37708 
+cos 2rr 0.40401 -+cos 227 0.43094+-cos 2rr 0.45787) +9(cos 227 0.48268) }?. 





* The observed intensities listed here were estimated by direct visual observation. The error in such a method is about 
one part in ten and greater accuracy than this is not claimed. 


factor F,, completing the sum and squaring. 


The complete mathematical expression for the 


intensity of the C,; alcohol is 





{8(cos 2rr 0.027+cos 2rr 0.05392-+-cos 2x7 0.08085 + cos 2rr 0.10778 








Let us see what this calculation tells us about 
the structure of this a-modification. In the first 
place, we have assumed that the chain axis was 
perpendicular or very nearly so with the 001 
plane. The success of this calculation tends to 
substantiate this assumption. Since we have 
found « =0.027 in this manner and since u =56,/d, 
we can find 6; which turns out to be 1.273A. This 
value is the distance from the origin to the center 
of the first carbon atom projection, which tells us 
that the projections of two molecules along the ¢ 
axis have a distance of 2.55A between projections 
of the carbon centers on the ends opposite the 
hydroxy! groups. Solving for 52 then gives us the 
distance between the projection of the hydroxy! 
group and the point d/2. This value of 4: is 
subject to some variation due to inaccuracies in 






intensity, summing the cosine values first over all 
the carbon atoms, then over the one oxygen 
atom, multiplying by the appropriate scattering 
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Fic. 3. Calculated intensities compared with observed intensities. The shaded values indicate 
calculated intensities. 


the experimental determination of d, but when 
calculated for all the available alcohols in the a- 
modification, the results shown in Table II are 
found. 


TABLE II. Calculated values of 52. 








Alco- 
hol Cis Cu Cis Cis Cis Cis Cie Ci Cro 
62 0.737 0.877 0.807 0.837 0.877 0. 753 0.667 0.497 0.497 


Average 62=0.73A 











252 is, of course, the distance between the 
projections of neighboring —OH groups. The 
fact that all the variation due to experimental 
error in the determination of d seems to appear in 
the value of 52, is only apparent, because the 
calculation could just as easily have been made 
using the point d/2 as the origin which would 
have caused 62 to appear invariant. It is readily 
seen that the same experimental error is inherent 
both in 6; and éo. 

So far we have derived all these quantities 
from the parameter found by calculating in- 


tensities for only the one alcohol C,;H;;OH. 
Precisely the same calculation, by using the same 
parameter values, should be applicable to the 
other members of the a-modification. The 
possibility of a change in the packing of the end 
groups from one member to another is quite 
remote, because not only do the c values all fall 
quite accurately on the same straight line? but the 
observed intensities all show the same type of 
alternation (Fig. 3). Thus, if we use the parameter 
6,= 1.273 with 6.=0.73 and apply the calculation 
to other members of this series and if the 
calculated intensities again check the experi- 
mental, this would offer good evidence in support 
of the theory. Such calculations have been made 
for the Cis, Cu, C15, Cig and Cig members which 
give the calculated intensities shown in Table I 
and when the comparison of the intensities is 
shown diagrammatically as in Fig. 3, the results 
are good; in fact the agreement is even better in 
some cases than with the C,; alcohol. 

Work is now in progress in this laboratory 
attempting to calculate intensities in a similar 
manner for the 6-modification. The results are 
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promising but the calculations are complicated 
by the presence of two parameters, one measuring 
the angle of rotation of the molecule around its 
long axis as well as the one measuring the 
position of the projection of the molecule on the 
perpendicular. This new parameter arises from 
the fact that the c axis of the 6-modification is 
tilted with respect to the 001 plane and conse- 
quently the projections of the individual atoms 
along the chain change their relative positions as 
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the molecule is rotated. 

It should be possible to extend the theory here 
presented to the calculation of intensities of 
reflection for all the long chain compounds, 
although many of the others offer complications 
which are not present in the case of the alcohols 
such as tilt, double bonded oxygens and so on. 
However, if the proper assumptions are made, it 
should be possible to extend this or a similar 
calculation to all cases. 
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By measuring electrically the ion emission through an orifice in an equilibrium chamber at 
high temperature containing salt vapor the lattice energies are calculated. The values obtained 


VOLUME 2 


are 151.3 kcal. for RbBr and 181.3 kcal. for NaCl. The probable error is about 3.0 kcal. for 
both measurements, and the maximum error in the neighborhood of 7.5 kcal. The corresponding 


electron affinities are 84.2 kcal. for Br and 88.3 kcal. for Cl. 





INTRODUCTION 


ECAUSE of the interest to chemists of the 
Born cycle for the calculation of energies 
involved in the formation of ionic salts, the 
experimental determination of the lattice energies 
of some alkali halides has been undertaken to 
extend the original work of one of us on KI 
and Cs]I.! Further improvements in the theo- 
retical calculation of lattice energies”: * have made 
desirable experimental work of the greatest 
possible accuracy. 
The method used for the experimental de- 


‘termination is essentially as follows. By deter- 


mining the current carried by ions flowing out 
of a measured hole in an oven containing salt 
vapor, the equilibrium pressures of ions P, and 
P_ of the ions M+ and X~ in equilibrium with 
the gaseous salt vapor MX at a pressure Pyx 
are measured. From these data the free energy 
change, AF, and by calculation of the entropies 
the energy change, AE, of the reaction (1) 
MX gas—>M*t gas t+X~gas are calculated. By em- 
ploying known energies of other reactions the 
lattice energy, AE for the reaction MX¢er. 
> M* gas t+X~gas, (O°K) and the electron affinity 
of the halogen, AE for Xgast+€~gas—>X ~gas (0°K), 
may be calculated. The actual handling of the 
data used in this paper, although in principle 





*From the dissertation submitted by L. Helmholz to 
The Johns Hopkins University in partial fulfillment of the 
degree of Doctor of Philosophy. 

‘J. E. Mayer, Zeits. f. Physik 61, 798 (1930). 

*M. Born and J. E. Mayer, Zeits. f. Physik 75, 1 (1932). 

*J. E. Mayer and L. Helmholz, Zeits. f. Physik 75, 19 
(1932). 
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assuming the same knowns as previously,' is 
such as to reduce greatly the uncertainty in 
the final lattice energies and electron affinities. 
The experimental procedure has also been greatly 
improved. 


EXPERIMENTAL 


The experimental arrangement has been altered 
from the original set-up of Mayer to give higher 
temperatures for the equilibrium (1) above. 
This was necessary in order to extend the 
measurements to the lighter alkali halides. One 
might instead use a more sensitive device for 
measuring the currents, but because of the 
difficulty of avoiding stray ion and electron 
currents this did not seem worth while. Modifica- 
tions have also been made in the experimental 
technique to improve the accuracy of the experi- 
ments, since in Mayer’s work the currents were 
not satisfactorily reproducible, and it seems 
probable, also, in view of our recent experiences 
that some unrecognized errors were inherent in 
the earlier procedure. 

Fig. 1 shows the essential parts of the appa- 
ratus and may serve to illustrate how the 
measurements were carried out. The salt to be 
studied was placed in the lower half of the 
graphite oven (A) which was threaded into the 
upper part (B). (B) was heated directly by 
radiation. (A) was heated principally by con- 
duction, so that its temperature could be ad- 
justed by changing its position relative to the 
upper oven. The gaseous salt molecules, at a 
pressure determined by 74, the temperature of 
the lower oven, diffused into B and then were 
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assumed to reach equilibrium with the ions and 
atoms at 7, the temperature of the upper half. 

The gaseous salt, ions, and atoms, diffused out 
of (B) through a measured hole in the platinum 
plate over a somewhat larger hole in the graphite 
oven. The pressure was such that Knudsen’s 
equation could be used to calculate the numbers 
of ions diffusing out. The ions were caught by 
the Faraday cage D and the current measured 
by a sensitive galvanometer in series with it. 
A potential of a few volts was applied to the cage 
to collect either of the two ion types dependent 
on the sign of the potential. 

The water-cooled shield C served to prevent 
ions and electrons from the heating elements and 
oven from reaching the cage, and also to select 
out a definite solid angle of the ion beam diffusing 
through H. Since salt collects on this part of 
the apparatus during an experiment the tube is 
water-cooled to prevent emission of ions or 
electrons from the surface which would other- 
wise be hot enough to permit this. A fine screen 
was placed at J to prevent the field of the 
Faraday cage from reaching through the opening 
and introducing uncertainty into the calculated 
solid angle of the ion beam that is collected. 
In practice the entire tube C was insulated and a 
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potential applied to it; in this way it was possible 
at fairly low voltages (2-5 volts) to collect 
effectively the entire ion beam, since the ions 
emerge from the oven according to a cosine law 
and with velocities only as high: as 0.05 volt. 
The screen removed 40 percent of the ion beam. 

The Faraday cage was insulated by two pieces 
of quartz tubing held against a brass support 
as shown in the diagram. The rod to the cage 
was threaded and nuts tightened to hold washers 
against the quartz which in turn was held rigid 
with respect to the brass support. The support 
was insulated from ground by three glass plates 
which were kept cool during an experiment and 
were also out of the path of diffusing salt 
molecules, so could not become coated and give 
rise to surface conduction. After each run the 
cage was tested for electrical leak and always 
when cool found to be well insulated. It is 
possible, however, that occasionally when the 
system was hot the leak may have been appreci- 
able. This would only materially effect the 
reading of low currents. 

Fig. 2 shows the apparatus as a whole. The 
system is mounted on a steel plate to which a 
bell jar about 30 cm in diameter and 40 cm tall 
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had been ground. With heavy grease it was 
found possible to hold the pressure to 10-* mm 
Hg during runs with the temperature of the 
furnace as high as 1300°K. The leads are brought 
out through copper tubes soldered to the plate 
through Pyrex tubing held in place with picein. 

A brass plate G held on three brass rods 
served to support the oven A, Faraday cage B 
and protecting tube D as well as the cooling 
jacket K. The cooling jacket was made from a 
large section of copper tubing 12 cm in diameter 
around the outside of which copper tubing was 
coiled and soldered in place. 

The heating elements EF of carborundum com- 
position (Globar) were set in a framework which 
could slide on the three supporting rods and was 
capable of adjustment with respect to the oven. 
Six of the Globar rods, ground down to have a 
resistance of 30 ohms each, were mounted in 
parallel and gave fairly uniform heating at high 
temperatures. The resistance of these sticks de- 
creases with temperature and they have the 
decided advantage over tungsten resistance ele- 
ments that they are unaffected by residual 
oxygen pressure. Because of the amount of gas 
given off by the graphite furnace with rising 
temperature this fact was of great importance. 
As an example of the performance of the furnace: 
at 60 volts it carried 12 amperes and the oven 
temperature rose to 1225°K. 

Several radiation shields (F) were interposed 
between the jacket and furnace to prevent loss 
by radiation. It was found particularly important 
to have the oven well shielded from above. 
Unless these precautions were taken differences 
of as much as 40° between the middle and the 
ends of the oven resulted. 

The oven itself was turned out of a solid block 
of graphite with an outside diameter of 5 cm 
and drilled out on the inside to 1.5 cm. This 
gave a wall thickness that insured good con- 
ductivity and helped to prevent any large tem- 
perature gradient from the middle to the ends. 
In spite of this thickness and good shielding the 
differences rose as high as 10° at 1350°K. The 
influence of this on the results will be considered 
later. The lower half of the oven, also turned 
out of one block of graphite was 3 cm in diameter 
at the base with walls 1 cm thick, and served to 
keep the temperature of the solid salt fairly 
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uniform. The threaded walls of the tube con- 
necting it with the equilibrium chamber were 
only effectively 1.5 mm in diameter so that 
large temperature differences could be main- 
tained between A and B. A small cap (shown 
in Fig. 1) was made to fit over the opening from 
A into B to prevent molecules from diffusing 
directly out of the oven without first suffering 
collisions with the walls in the upper oven. 

C indicates the method by means of which the 
lower half of the oven could be raised or lowered 
during an experiment without breaking the 
vacuum. The iron rod screwed to the bottom 
of the furnace was connected to a ground joint 
fixed in the base plate by means of a surrounding 
tube in which it could slide. 

Three thermocouples were employed to meas- 
ure the temperature of different parts of the 
oven. They were C.P. platinum-platinum 10 
percent rhodium couples and were calibrated at 
five points up to 1083°C and a table of deviations 
from the Bureau of Standards table recorded. 
One at a depth 1/3 and another at depth 2/3 of 
the oven were used in the upper part to determine 
the temperature of the equilibrium and to check 
the uniformity of the temperature along the 
length of the oven. The cold junctions of these 
thermocouples were insulated and fastened to a 
brass block attached to the inlet cooling water. 
A thermometer inserted in the block in the same 
way enabled us to make a correction for the 
difference in cold junction temperature from 
the standard. This was usually 20 or 25°C. 

The system was evacuated with a three stage 
steel mercury diffusion pump which was, of 
course, in continuous operation. It reduced the 
pressure to 10-5 mm when the system was cold. 
The pressure was determined by means of a 
McLeod gauge attached to the system as shown 
in Fig. 2. 

The galvanometer was the best Leeds and 
Northrup type M 8. rated at 10-"' amperes per 
mm scale deflection at 1 meter. As set up in 
this experiment it gave a sensitivity of 2.17 
X10-" amperes per mm on calibration. The lead 
to the galvanometer was carefully insulated and 
came in contact with glass only when it emerged 
from the vacuum system. It was sealed in with 
de Khotinsky cement. The rubidium bromide 
used was Kahlbaums which had been crystallized 
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three times by Professor Jones and used in con- 
ductivity work in the department. The sodium 
chloride was used as a standard in x-ray analysis. 


CALCULATION OF RESULTS 


From the observed ion currents and the 
knowledge of the two temperatures 7,4 and T, 
of the solid salt in A, Fig. 1, and of the vapor in 
B, as well as the geometry of the holes and 
measuring system, it is possible to calculate 
directly the lattice energy. For the reaction 
taking place at A, the sublimation of the salt, 
we have 


F4°(MX gas) -_ Fx°(MX eryst.) =— RT, In P MX (1) 
and for the equilibrium at B, 


Fxz°(ion gas) — Fg°(MX gas) 
= —RT1n (P,P_/Pemx) (2) 


where Fr? is the standard free energy of the indi- 
cated substance at the indicated temperature, 
and P,, P_, and Pyx, the pressures of the posi- 
tive and negative ions and of the salt molecules 
(in atmospheres). It is assumed that due to the 
Knudsen effect Pawx = (T. /Tp) *Prux. 
Multiplying Eq. (1) by 7s and (2) by 74 and 
adding, the relatively uncertain value of Pyrx is 
eliminated, and the less uncertain difference 





TeF°r, (MX gas) —T, F°vy( MX gas) introduced. 
For the various free energies the expression 


id 
Fr? =Ey+ f CaT~-TS (3) 
0 


is used, where E is the energy at 7=0°K. The 
energy Eo(MXeryst.), the energy of the crystalline 
salt at the absolute zero, is arbitrarily chosen as 
zero, FE,(ions), the energy of the ions at the 
absolute zero, is then the lattice energy. The 
specific heat integral of the salt is calculated 
from a Debye function, RTD(@/T), 0=Avmax/k, 
and the entropy of the salt RSp(@/T) is similarly 
calculated assuming a Debye specific heat curve. 
Probable values for the ymax have been taken 
from an article by the authors.* E(sub.), the 
energy at 0°K of the gaseous molecule, has been 
taken from the same source, as have the probable 
values for the frequencies of the molecular 
vibration. These latter values for the gaseous 
molecules are all of negligible importance, since, 
as before remarked, only the difference of 
TFr°(MX gas) at two not very distant tempera- 
tures enters the equation. If ¢€ is the energy hy 
of vibration of the gaseous molecule, the in- 
tegrated specific heat due to vibration may be 
written RTP(hv/kT) and the vibrational entropy 
RSg(hv/kT) where P(hv/kT) and Sze(hv/kT) are 
appropriate well-known functions. The equation 
for the lattice energy then becomes 


Ta-T . - : ; 
E,(ions) = -—*——*zy(sub.)+RT {P(—) -P(—) i {Se(— ) —S (— )} 
A ad . . ‘ 


0 6 
+6|D(—) -5»(—) —7.314+1.5 In Wa Wx+2.0 In Tg+3.0 In T4 —In P.P-| (4) 


Ts 


A 


Wy and Wy are the atomic weights of the two ions. The ionic gas is assumed to be unimolecular 


with only the Stern-Tetrode entropy. 


The pressure of the gaseous ion is related to the galvanometer deflection by the relation 


(2eWRT>)* 2.17X10-"X6 





1.013108 AX FX96,500 


6 
= 5.07 X10-'8(WT,)? = (5) 


Here 6 is the galvanometer deflection in mm and 2.17 X10" the sensitivity of the galvanometer in 
amperes per mm. A is the area in cm? of the oven opening, and F the calculated fraction of the ions 
recorded. F is of course dependent on the voltage placed on the screen and Faraday cage. The 
pressure is given in atmospheres. W is the atomic weight. 
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The pressures may be replaced in (4) by deflections, changing to logiy we obtain 


Te—T, 1 € € 
E)(ions) = -—*—“a(sub) +4.5987y 31.36+——| P( ) -P( ) 
T 2.303 kT, kT 4 


A 


1 { € € 6 6 
-|s,(+)-s(*)] 42a08|0(")-ss(*)] 
2.303| \kTy kT 4 + T, 


+log Wu Wx +log T,+3 log T7,4+2 log A+2 log F—log i.0-| (6) 


This method of calculation, due to the partial elimination of the constants of the molecule gas, 
considerably reduces errors present in the original work on CsI and KI. 


DISCUSSION OF THE METHOD AND ERRORS 


In order to give an idea of the relative im- 
portance of the factors involved in the calcula- 
tions the numerical values for experiment (1) 
in Table II are reproduced in Table I. 


TABLE I. 








26.0 kcal. 
6.04 keal./2.303 
31.36 
3.83 
3.12 
8.80 
— 6.87 
—0.12 
—5.85 
—4.68 
29.59 
178.5 kcal. 
152.5 


((Te—Ta)/T aE (sub.) 
4.598 Tp 


numerical constant 
log Wy Wx 

log Tz 

3 log Ta 
2.6[D—Sp]r-ra 
(1/2.3)[((Ps—Pa)—(Sa—Sx)] 
2 log A 

2 log F— log 6,5_ 
sum 

4.598 TpXsum 
Lattice energy 








It is evident that an error of 5 percent in 
Ji’C,dT—TS for the crystalline salt, corre- 
sponding to the same percentage error in S(6/T) 
—D(@/T) would introduce about 2 kcal. error 
into the final lattice energy. This is possible, and 
the greatest single questionable item in the 
calculations alone. A fairly large error in the fre- 
quency of the gaseous molecule would have a 
negligible effect. 

The 5 percent uncertainty (2.5 kcal.) in the 
energy of sublimation, E)/X, introduces only 
1 kcal. uncertainty in the lattice energy. About 
3 kcal. maximum error, and probably con- 
siderably less, is then introduced by uncertainties 
in data unconnected with the experiment itself. 

The temperatures are undoubtedly measured 
correctly in the experiment to better than 5°, 
which would introduce about 0.25 percent or 





about 0.5 kcal. uncertainty. Although the de- 
flections are measured with considerable accu- 
racy the fraction of the beam measured, and the 
diameter of the hole necessarily introduce some 
25 percent error in this term, or about 1 kcal. 
in the lattice energy. 

One may then expect a maximum error of 
+3 kcal. due to uncertain values used in the 
calculations +1.5 kcal. due to inherent experi- 
mental errors, provided only that the current 
measured is due to ions emitted from the oven 
and provided Eq. (5) relating the current meas- 
ured to the pressure of ions in the oven holds 
correctly. 

It is difficult to state exactly all the reasons 
adduced by the authors for believing that 
actually only ion currents from the oven were 
measured in the experiments reported. In various 
experiments not here reported stray currents 
undoubtedly existed. In many runs when salt 
had previously coated the outside of the oven at 
H in Fig. 1 abnormal positive currents and some 
negative currents were observed at even com- 
paratively low temperatures. These currents 
dropped rapidly with time, presumably as the 
salt and other possible impurities distilled off 
the platinum plate. In other experiments when 
the shield G of Fig. 1 was absent abnormal stray 
currents were observed. Certain elementary 
reasons for believing that the currents observed 
in the reported experiments were due only to 
ions from the oven orifice may be given. 

(1) The experiments were reproducible when 
no obvious difficulty arose. 

(2) Without salt in the oven no deflections 
were observed. In one case during a run the salt 
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completely evaporated from the oven and the 
deflections fell to zero. 

(3) Over a range of 125° and sevenfold change 
in current no systematic change in the calculated 
lattice energies is observable. 

(4) Variations of the current with voltages 
applied to the screen J and Faraday cage D of 
Fig. 1 agreed qualitatively with the expected 
behavior of ions emerging from the hole H with 
about 0.05 volt energy and directions dis- 
tributed according to a cosine law. 

(5) Current with NaCl and RbBr in the oven 
differed by 35-fold in agreement with probable 
theoretical values of the lattice energies. 

With regard to the second factor, the validity 
of Eq. (5) only negative remarks can be made. 
It is to be expected that what may be called a 
contact potential will exist between the gas in 
the interior of the oven and the platinum 
surface surrounding the hole H. An ion sheath 
will clothe the surface. An ion will then, de- 
pending on its sign, either be accelerated toward 
the hole H or repelled from it. It is consequently 
not to be expected that exactly equal numbers of 
negative and positive ions, or exactly numbers in 
the ratio of the square roots of the molecular 
weights, will emerge. For one ion species the 
deflection will be larger than that corresponding 
to the existing pressure, and for the other smaller. 
This effect is actually observed. The positive ion 
currents are always considerably larger than the 
negative ion currents, and the factor is in some 
cases as large as three. However, if this interpre- 
tation is correct the approximation of using the 
product of the two deflections in Eq. (6) will 
not be greatly in error. It is to be emphasized 
here that a threefold error in the estimation of 
the pressure of one of the ion species introduces 
only a 3 kcal. error in the final lattice energy, 
an error of 2 percent. 

It is hoped that a slightly altered experi- 
mental technique and apparatus will enable us 
to completely obviate this uncertainty. 

Little has been said in the above about the 
factor F in Eq. (5) giving the fraction of the ions 
caught by the Faraday cage. The screen J 
transmits only a certain fraction dependent on 
its transparency. If the true potential difference 
between H and J (Fig. 1) were exactly zero, a 
certain fraction calculable from cosine law 
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emission and the geometry of the apparatus 
would hit the screen. As was to be expected this 
true zero of potential difference was found not 
to be zero as read on a voltmeter, but different 
(by 0.3 volt) due to the contact potentials and 
high temperature gradients existing at the con- 
necting parts. If the potential of J were raised 
with respect to 7 (in either direction) one of the 
ion species would be accelerated perpendicular 
to the oven face and the fraction of these ions 
striking within the diameter of the hole at / 
should rise asymptotically to 1. This fraction 
should approach very nearly 1 at about 1 or 2 
volts. This behavior was observed. By this 
method the authors believe that they were able 
to estimate the deflections when 100 percent of 
each of the ion species were striking within the 
hole radius 7 to within better than 20 percent. 


RESULTS 


The result of 10 experiments on RbBr and 2 
experiments on NaCl are given in Table II. 
The lattice energies of RbBr and NaCl deduced 
from the experiments, 151.3 and 181.3, respec- 
tively, agree well with the theoretical values 
153.5 and 183.1 previously calculated by us. 
The authors would be inclined to set an upper 
limit of 4.5 kcal. on the possible error of the 
experimental values, about 3 kcal. of which are 
due to uncertain data for the specific heat of 
the salts and the heat of sublimation. In addition 
there is, however, a possible error of some +3 
kcal. due to possible faulty interpretation of the 
results, which error, however, if it exists, will be 


TABLE II. Lattice energies of RbBr and NaCl. 








Lattice 
energy 
(kcal.) 


Deflections 
(6, mm) (6. mm) 


169 103 
188 
84 


Te i 
(°K) (°K) 


1316 860 
1319 868 
1262 831 
1260 837 
1251 827 
1244 830 
1230 827 
1226 827 
1211 810 
1197 824 
1264 945 ° 
1160 932 . 1.5 


Average lattice energies, RbBr 151.3, NaCl 181.3 kcal. 





152.5 
1532 
149.4 
152.2 
150.1 
152.0 
152.0 
151.7 
149.5 
150.2 
181.7 
181.0 


NaCl 














> «de> te de eee. de + 4 


NRUNOOH RRS 


LATTICE ENERGIES OF RbBr 


systematic in the sense that it has the same sign 
and closely the same order of magnitude for 
both salts. The usual method of calculating 
probable error by root-mean-square method 
leads to +2.5 kcal. neglecting the possible 
systematic error, and +4 kcal. including this. 
Owing to the small number of experiments 
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carried out on NaCl there is an additional un- 
certainty in this case. 

With the Born cycle, and the particular values 
of AE for the involved reactions previously used 
by us, the electron affinity, AE for the reaction 
X gas te gas PX ~ eas O°K, is calculated to be 84.2 
kcal. for Br and 88.3 kcal. for Cl. 
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The Mutual Repulsive Potential of Closed Shells! 


WILLARD E. BLEICK? AND JosEPH E. MAYER, Chemical Laboratory, The Johns Hopkins University 
(Received March 5, 1934) 


A method is developed for representing, in terms of a small number of integrals, the mutual 
potential of two ions or atoms having a rare gas electron configuration. The accuracy is the 
usual first approximation of the Heitler-London method. The integrals are evaluated for two 
like atoms, and numerical calculations made for two neon atoms. The total (repulsive) potential 
at R=1.8X10-* cm is 344 X10 erg; at R=2.3 X 10-8 cm, 35 X10- erg; and at R=3.2X10-8 


cm, 0.410- erg. 





INTRODUCTION 


HE calculation of the potential between two 
rare gas atoms by a generalized Heitler- 
London* method is a problem which in some 
respects is unique. First, the absence of degener- 
acy makes the solution comparatively simple in 
spite of the large numbers of electrons. Secondly 
the results are, almost without alteration, ap- 
plicable to the repulsive potential between ions 
having a rare gas electron configuration, and so 
refer to a large number of chem‘7al compounds. 
Thirdly, the experimental conditions under which 
this potential is observed: comparatively large 
distance between the ions and consequently low 
values of the potential, and in crystals sym- 
metric surroundings for the ions leading to a 
minimum of distortion; lead one to believe that 
the theoretical calculations should give fairly 
accurate agreement with observed potentials. 
The mutual potential of two helium atoms has 
been calculated by Slater* and others, but no 
adequate treatment has been given for the case 
of eight electrons. 
In this paper the potential for two atoms or 
ions having closed shell configurations of eight 


1 From a dissertation submitted by W. E. Bleick to the 
Board of University Studies of The Johns Hopkins 
University in partial fulfillment of the requirements for 
the degree of Doctor of Philosophy. 

2 J. T. Baker Chemical Company Fellow. 

3 W. Heitler and F. London, Zeits. f. Physik 44, 455 
(1927). 

4J. C. Slater, Phys. Rev. 32, 349 (1928). 


electrons is represented in terms of a com- 
paratively small number of integrals. The evalu- 
ation and numerical calculations are carried out 
for the case of two neon atoms. In a subsequent 
paper semi-empirical approximations of the 
integrals will be given which enable a simple and 
moderately accurate functional representation 
of the repulsive potential of two (rare gas shell) 
ions with known outer electron distribution. 


THE INTEGRALS 


The potential is determined as a perturbation 
energy of two atoms in a manner exactly anal- 
ogous to the method of Heitler and London.’ 
Two atoms a and b of core charge 2, and 2s, 
assumed to be point charges, have a nuclear 
separation R. There are sixteen electrons (twelve 
p and four s electrons) numbered from 1 to 16.’ 
Yaz and 7,; are used as symbols for the distances 
from the nuclei a and 6 respectively to the 
electron k. The zeroth order wave function is a 
linear combination of permutations of the 
product of one-electron hydrogen-like functions, 


5In order not to make the succeeding mathematical 
presentation too cumbersome 8 electrons and a core charge 
(+8 for neutral atoms) are used instead of the more exact 
treatment with all electrons. The final results are so 
simple that the exact equations can be written readily. 
The extra integrals which would be introduced by the 
inclusion of the electrons underlying the closed shells can 
be seen to be negligible. The integrals which are evaluated, 
however, are modified to allow for the existence of the 
inner electrons. 
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REPULSIVE POTENTIAL OF CLOSED SHELLS 


16! terms 


Wo= (1/164) Lv 


i= Wa0,a(1) . ¥a0,8(2) : Ya—a(3) ¥ va—B(4) 7 Ya0p(5) ae 


where Wao,, Ya, Wao, and a; are the ms and np 
one-electron (Hartree) functions of atom a with 
m,=0;,, —1,0, and +1, respectively, and a and 
6 represent the two possible spin functions. The 
zero m;’'s are given the subscript s or p according 
as the m, refers to an s or p state. 

P is a permutation operator, operating on the 
product of the 16 different functions of the sixteen 


Cpe p, 


Wp=P(W)), a) 


¥a+B(8) ° ¥00,a(9) eee ¥o+8(16), 





cp has the value +1 when the permutation P is 
even, and —1 when P is odd. The symbol W;, 
is used for the identical permutation of the 
product, written out as in Eq. (1), with electrons 
1 to 8 on atom a, and 9 to 16 on atom b. Y"p is 
used for this product operated on by the per- 
mutation P. 
The Schrédinger equation may be written 


electrons; 16! different permutations of the 16 
electrons among the 16 functions are possible. 


H(o) =(T+ V)(Wo) = EW. (2) 





The potential energy V has the form, 


k=16 1=16 | k=16 f 2, 25 
v=e = ; —_—— > —4+—)+ 


k=1 l=1 Tiel k=1 Tak T bie 


ZaZb 
| (3) 
R 
If this value for V is substituted in (2) we may follow Heitler and London and cancel, in every 
permutation occurring additively in Wp», the kinetic energy term plus that part of the potential 
energy which for the particular permutation occurs in the isolated atoms, operating on the permuted 
product p with the sum of the energies of the two atoms Eo.+ Fo multiplied by Wp. There remains 
an equation in which the perturbation energy, and what may be termed a perturbation potential, 
occur. The form of the perturbation potential depends, however, on the permutation of the term 4p 
which follows it. Calling AV, the perturbation potential for the identical permutation term , the 
perturbation potential AV; for the term “p is obtained by operating on AV, with P. 


i k=8 2, Za2 
+ =| 
R 


1=16 Z, 


=9 Tal 


avize| | = 


Tei k=l Lor 


AVp=P(AV)). 


The procedure outlined above leads to the equation, 


16! terms 


(1/16 13) Zz cpA VpWp=AEW>. 
a 


Multiplication by YW * and integration over the whole space r gives 


16! terms 16! terms 


(1/16!) > » ¥ coc | We* Vow pdr = AE { Weeds = AE(1—S). 
P 


P’ 


S corresponds to the S of Heitler and London and has a small value (~0.01). 


16! terms 16! age 


—S=(1/16!) Y hy CpCp: [renee dr 
a 
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in which >~’ indicates that P and P’ never take 
identical values. As will become evident later 
only permutations of one pair need be considered 
and the alternate Eq. (9) may be used, 


1 
S=200 2D Si, (} =o -1, 0, +1) (9) 
i j 
Si= | $ij|? 
su | vostdadr, 


a and j being indices that in this case indicate 
the value of the quantum number 1;. 

The sum on the left of Eq. (7) contains (16 !)? 
terms. After integration the numerical value of 
any term is, of course, independent of the 
numbers assigned to the electrons. We, therefore, 
change the notation of the electrons in each term 
taken as a whole, so that Wp, and AVp- become 
Ww, and AV;, operating on the term with (P’)~, 
obtaining 16! different terms each occurring 16! 
times. This gives 


16! terms 
} Cp [wera viwdr= aE —S). (10) 
P ~ 


Of the 16! terms in Eq. (10) less than seven- 
teen are different from each other and appreci- 
ably different from zero. 

First since AV, is not a function of the spin 
coordinates and since a and £6 are orthogonal, 
JS a*8=0, terms in which any electron in Wp* 
has a different spin than in , are zero. 

Secondly, since the perturbation potential 
contains no terms with the coordinates of more 
than one electron on the same atom in W;, and 
since J Wai*¥ajdt=0 when 7+j all terms in 
which any electron has a different place (differ- 
ent y) on the same atom in Wp* and W; are zero. 

Thirdly, since the perturbation potential con- 
tains no terms with the coordinates of more than 
two electrons, and since J Wai* jdt =5;;K1 we 
may neglect terms in which more than one pair 
of electrons has different positions in Wp* and 
W,, because these will all be multiplied by 
factors S;;= |si;|?<1. 

Odd numbered electrons (spin a) may only 
be permuted with odd numbered electrons, and 
even with even. Fourthly, for every permutation 
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of one pair of odd numbered electrons there will 
be a corresponding permutation of even num- 
bered electrons in which electrons the same pair 
of functions y,; and y,; are permuted. The two 
terms will be numerically the same. 

Of the 16! different terms Wp* we need now 
only consider the identical permutation W,* 
and sixteen others. These sixteen others are 
those in which one of the four odd electrons 
(spin a) on atom a in “&, (electrons 1, 3, 5 or 7) 
is exchanged with one of the four odd electrons 
(spin a) on atom 0 in W, (electrons 9, 11, 13 or 
15). Each of these sixteen terms must be mul- 
tiplied by a factor 2 to allow for the correspond- 
ing numerically identical term of electrons with 
spin 8 (even). The term of the identical permu- 
tation will have the plus sign for cp; the other 
sixteen terms will have the minus sign. 

Due to the analytical form of the functions y 
several of the sixteen terms will be numerically 
identical. For instance if the atoms a and 3 are 
the same, there are only seven different terms. 
If the permuted pair of electrons had both 
m,’s=0, one term results; and similarly if both 
m,'s=0,. Two terms with the m,’s=0, and 0, 
are identical. Four terms with the m,’s=0, and 
+1 or 0, and —1 are identical; and similarly if 
the m,’s=0, and +1 or 0, and —1. If both m,’s 
were +1 or both —1 the two terms are identical. 
Finally two terms with the m;’s= +1 and —1 are 
identical. 

In order to further simplify the integrals it is 
necessary to introduce twelve new functions, 
six of the coordinate 7, the distance from the 
nucleus a, and six functions of the coordinate 7». 
The functions are pa, po, Pai, prj’, Ua, Urs, Uai’ 
and U,,’. 

pa is the density of electric charge (in units 
of e) in the atom a. pq,’ is the density of charge 
in the atom a with the 7’th electron of the closed 
shell missing. U, is the potential due to atom a. 
U.i' is the potential due to a when the7'th 
electron of the closed shell is missing. pai’ and 
U.:/ are not to be confused with the correspond- 
ing functions of a singly charged ion of the 
element a, where the electron distribution has 
altered due to the altered field. They are cal- 
culated assuming the electron distribution of the 
atomic functions, neglecting one electron of the 
closed shell. Strictly speaking in the case of the 
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np electrons p.;' and U,,;’ are not functions of 7a 
alone, but also of 0, and ¢,. They depend on the 
quantum number m, of the missing electron. 
This neglect of the angular part of the function 
in the field is the same assumption that is made 
in the calculation of the Hartree functions which 
will be used in the evaluations. 

The equations for p, and p,;’ are the following: 


pata) = — Ds Hail? +2a(5(ra)/4e70”), (11) 
all electrons 


Pai (Ta) = Pa + | Raila) | 2/4 7, 


where 6(r2) is Dirac’s delta-function,® and 
Rai(’a) is the radial function (normalized) of the 
i’th electron of the closed shell. 

The potential U,(ra) is given by 


1 7% co 
Uulra)=— f Aartparhdr+ ff Amrp.(r)dr, (12) 
0 Ta 


Ta® 


which is identical with 


z Wai*(ri) "Wadi, (13) 


all electrons 


U.(ra) =Ze ta 


where r; is the distance from some point with the 
coordinate r, to the 7’th electron, and dz; is the 
volume element of the 7z’th electron. 

The potential U,;’(r.) is given by an equation 
identical with (12) in which pz,’ is substituted 
for pa. This equation is identical with Eq. (13) 
in which the sum goes over all but one electron. 

With the help of these functions, which it 
should be noted are easily calculated when the 
Hartree functions of the atom are available, we 
consider the seventeen terms of the sum on the 
left-hand side of Eq. (10). 

The first term, in which the two identical 
permutations occur, is the classical Coulomb 





interaction, ef U.pdr=ef Uspadr. We may 
designate, in accordance with Heitler and 
London, this as W;. It has a negative value in 
the range of R we are considering, occurs with 
a positive coefficient, and so leads to attraction. 

The remaining exchange energy terms are not 
so trivially evaluated. Let us consider the term 
Wp* in which the electron pair 1=i=8 and 
9 =7=16 have exchanged the places they have 
in &,. All other electrons have the same positions 
(same y’s) as in W@W). 

In this term we consider first those (additive) 
terms in AV;, Eq. (4), which contain neither 2 nor 
j, that is where k #2, | +7. The partial integration 
of these terms gives 


sue [ Uai' po;'dt = Sve [ Un! pad T. 


The summation over all possible exchanges 7 
and j gives 2e?>0 ;0 jSi; Sf’ Uai' po;'dr. Since U,;’ and 
p»; do not differ greatly from U, and ps, and 
S is very small we introduce no great error by 
writing this SW,,;. The terms all occur with a 
negative coefficient. 

In the perturbation potential we have left the 
following terms: first >>)’ (7::)~'—2,/r»:, secondly 
Dx (rj) —Za/Pa;, and thirdly (r;;)—!, where the 
symbol >’ is used to indicate that the sum over / 
does not include 7, and the sum over & does not 
include 7. The partial integration of the first of 
these terms gives —s;;e? { U2; Wo;*Wadr, and of 
the second —s;e? Sf Uai Wai*sjdr. The third 
term cannot be further simplified. All these 
terms occur with negative coefficient. ~.; and 
y»; are the functions in which the electrons 7 
and j occur in @&). 

The results of these considerations are summed 
up in Eq. (14). 


AE=Wiu—[2/(1-S)] Yd (siUi’ + Ri), (| =0,, —1, 0,, +1), (14) 
i 7 J 


Wax [ Usosdr=e { Uspadr, 


U;;/' = —é [ (Ua'+ U>;')Vai*Pojd7, 


Ryt=e { Wai Poj(1) (712) Wo;* Pai(2)d 71d 72. 


niDciemnsens 


*P. A. M. Dirac, The Principles of Quantum Mechanics, p. 63, Oxford (1930). 
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The symbol dr is used as elsewhere to indicate 
the volume element, and all integrations are to 
be carried out over all space. s;; and S are defined 
in Eq. (9), and the U’s and p’s in (11), (12), (13). 

To proceed further it is necessary to have 
analytical one-electron functions for the atoms 
a and b. That is, functions whose product satisfies 
the radial part of the Schrédinger Eq. (2) for the 
single atoms.a and 8, and which are orthogonal 
and normalized. Slater’ has given such analytical 
approximations for the Hartree® functions. It is 
obvious that if the summation over 7 and 7 in 
(14) includes the inner electrons the expression 
is that which would be obtained if these electrons 
were not neglected. The expressions (11), (12) 
and (13) for p and U are taken to include all 
electrons. Actually it is here advisable, par- 
ticularly in calculating U,,’ and U,;’, to include 
the inner electrons, and this does not add appre- 
ciably to the work of evaluation. It is, however, 
numerically unnecessary to include the inner 
electrons in the summation of (14). By the use 
of elliptic coordinates the calculation of W,; and 
U;;’ is straightforward and simple, although 





p(r) = 10(6(r) /42r?) — 595e—1% 3" — 7? (15 5e-* 5" + 3. 56e7*477), 
Pop (r) = 10(6(r) /4xr?) — 613e-1% 5" — 7? (132e-5 997 + 2. 88e--457), 
pas (r) = 10(8(r) /4-r?) — 595e-19-37 — 72(134e-% 47 4+. 2. 893-417) | 


BLEICK AND J. E. 
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somewhat long. The calculation of R;;" is 
difficult. 


NEon ATOM FUNCTIONS 


The Hartree field used for neon is that of 
Brown.’ The following analytical representations 
of the 2p and 2s radial functions were employed: 


Ro»p(r) =7(17.9e-* 89 + 2.30e-1-69"), (15) 


Ros(r) = 13.6¢e-* 2" 


—r7(14.7e-3-°9" 4+ 4.76e-? 5"), (16) 


The 2) representation differs slightly from that 
given by Brown and is a better approximation 
to his calculated field for large radii. The unit 
of length in these and the following functions is 
the Bohr atomic unit, i.e., 0.528 10-° cm. 
Brown tabulates values of the charge density 
p and of the 2 and 2s radial functions. It is easy 
by substraction of (2)? and (2s), respectively, 
from p to calculate the modified charge densities 
pop and pp,’ of (11). Approximate analytical 
representations of these densities are 


(17) 
(18) 
(19) 


The integrations of p, pe,’ and pe,’ according to (12) and (13), yield the following potentials: 
U(r) =.(20.0+ 2.087-!)e13" + (53.272 +52.87+ 26.2 45.7797 )e* 


+ (3.7377+6.467+5.59+2.157-")e* 47", (20) 


Vey (r) =1/r+(20.2+2.077-) e195" + (46.17?+ 46.27 + 23.1+5.157)e 99” 


+ (3.047?+5.297+4.60+1.787r)e*", (21) 


U2,’ (r) =1/r+(20.1+2.087-!)e—!3" + (46.17°+45.87+ 22.8+5.02r-1)e-*-04r 


CALCULATIONS 


In the calculation of the integrals it is con- 
venient to introduce the elliptic coordinates p, 
q defined by 


p=(ratre)/R, g=(ta—1s)/R. 

7J. C. Slater, Phys. Rev. 42, 33 (1932). 

8D. R. Hartree, Proc. Camb. Phil. Soc. 24, 89, 111 
(1928). 


+ (3.1377+5.507+4.84+1.897-)e# 417, (22) 





As the third coordinate the azimuth ¢ about the 
line of nuclei is taken. The volume element is 
(R/2)3(p?—@)dpdqdo. It is assumed that the 
polar axis of each atom is directed toward the 
other atom along the line of nuclei. Then @, and 
6,, the angles between the radius vectors 7. and 
ry respectively and the line of nuclei, are given by 


9F. W. Brown, Phys. Rev. 44, 214 (1933). 
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sin 64=[(p?—1)(1—¢*) ]'/(p+4), 
sin 6,=[(p?—1\(1—q*) ]}/(p—4), 
cos 04= (1+ p9)/(p+q), 
cos 6,= (1— pq) /(p—q). 


In these coordinates all of the integrals except 
R;;-' reduce to a sum of terms of the type 


p=o Q=+1 
Tf oun 
p=+1 q=—1 
= (—1)™(e8-*/Ba)A,(a)Bn(B), 
(n,m>0;a>0), 
where 


A,(a)=¥ n!a-i/(n—i) |, (23) 


Bn(B) =A m(—B)+(—1)™*1e*8A »(B), 
B,,(—B) _ (-— 1)™*1¢e?6B,,(B), 


and a and £@ are proportional to R. The poly- 
nomials A,(@) are conveniently calculated from 
the recursion formula 


A,(a) =14+NA n_-1(a@)/a. (24) 


The integrals R;;-' have been evaluated for the 
2p electrons of like atoms and are given in an 
appendix. 


NUMERICAL RESULTS 


The numerical values of the integrals (9) and 
14) for two neon atoms are tabulated in Tables 
land II. The nuclear separation of 3.2 10-* cm 
corresponds to the equilibrium separation of 
atoms in the neon lattice at the temperature of 
liquid helium. The separation of 2.310-*8 cm 
is the same as the distance of closest approach of 
unlike ions in the sodium fluoride lattice at room 
temperature. 


DISCUSSION OF THE RESULTS 


Lennard-Jones” gives several choices of values 
for the three constants A, B and n of an empirical 
tquation ¢= —A/R*+B/R" for the total poten- 
tial between two neon atoms as deduced from 





"J. E. Lennard-Jones, Proc. Phys. Soc. 43, 461 (1931). 
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TABLE I. Numerical values of the integrals for two neon atoms. 























R —Ujj'/e? —sjiUij’/e2 Rij/e? 
(cm Sij Sij (cm~! (cm~! (cm~! 
x 108) 1 j ( X10") (X 104) x 10-4) « 1074) x 10-4) 
Op Op 11.6 133 4184 483 203 
Op 0, 9.73 94.7 4174 406 116 
Os 0, 5.02 25.2 2069 104 30.8 
1.8 +1 +1 3.87 15.0 1083 41.9 13.1 
0, +1 O 0 0 0 14.2 
1 Ox 0 0 0 0 6.15 
+1 +-1 0 0 0 0 1.41 
Op Op» 5.43 29.4 1502 81.5 37.5 
Op Os 3.42 11.7 1166 39.9 13.1 
0, Os 1.31 1.71 412 5.39 1.91 
2.3 +1 +1 1.24 1.54 263 3.26 1.23 
Op +1 0 0 0 0 1.67 
+1 0, 0 0 0 0 0.403 
+1 +1 0 0 0 0 0.129 
Op Op 0.984 0.968 173 1.70 0.982 
Op 0, 0.410 0.168 102 0.418 0.162 
0; O, 0.0930 0.00865 19.7 0.0183 0.00832 
3.2 +1 +1 0.141 0.0200 21.4 0.0302 0.0138 
Op +1 0 0 0 0 0.0259 
+1 Os 0 0 0 0 0.00245 
+1 +1 0 0 0 0 0.00143 
various experimental measurements. Two sets, 


n=9, A =13.7X10-, B=3.27 X10-” and n= 12, 
A=8.45 x10", B= 36.3 X 10-™, have been used 
to calculate the “‘experimental’’ values of the 
potential given in the first two columns of Table 
III. Needless to say the extrapolation of these 
equations to R=1.8X10-* cm has not much 
value. In order to calculate the corresponding 
theoretical potential it is necessary to include a 
term which would be added by the second ap- 
proximation of the method used here. At fairly 
large distances this term can be readily approx- 
imated, as shown by London," using a develop- 
ment of the perturbation potential in R. The 
first term becomes —A/R®, A =3a°?J/4. a@ is the 
polarizability, and J an energy corresponding to 
the most intense absorption frequencies. Using 
a=0.392X10-* cm* and J=25 volts, A =4.6 
x10-° erg cm*®. By adding this term to our 
previously calculated repulsion the third column 
of Table III is obtained. 

The certainty of the experimental values of the 
potential is somewhat less than that indicated by 
the agreement between the two columns of 
Table III. Obviously the attractive potential is 
actually greater than that used to calculate the 
column marked ‘‘theoretical.’’ This is due partly 
to the rough approximation used to obtain the 
constant A, and partly to the neglect of higher 
order terms. Table III indicates, however, that 
the values of the repulsive potential calculated 


11 F, London, Zeits. f. physik. Chemie 11, 222 (1930). 
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Sums. 








R post2snU'y To 
(cm X 108) 


2:2;Ri;7! 


—2 
ics Lri(siVistRii') —-Wu 


1 
(ergs X 10") 





729 


8 283 
3 79.7 
3.2 


35.3 
0.663 


446 102 
44.4 9.45 
0.526 0.101 








in this paper are not in definite disagreement with 
experiment. A much better comparison with 
experiment could, of course, be carried out if the 
repulsive potential between an alkali and a 
halide ion were calculated. 

It is interesting to attempt to fit the values 
of Table II to a simple two constant equation. 
This has been done using the two extreme values 
to calculate the constants and comparing the 
equation with the value of Table II at R=2.3A. 
The equation 9(R)=B/R"; B=2.71X10-™ in 
ergs and cm, »=11.63, fits the two values at 
R=1.8A and 3.2A but gives 19.810-" erg at 
2.3A instead of 34.9X10-", a deviation of 40 
percent. The equation ¢(R)=6 exp (27o—R)/p, 
b arbitrarily chosen as 10-” erg, ro = 1.029 10-8 
cm and p=0.2091 x 10-® cm fits the two values 
at 1.8 and 3.2A and gives 31.5xX10-" erg at 
2.3A, an agreement to within 10 percent. Ob- 
viously the exponential equation gives a better 
fit than the inverse power. 


TABLE III. The total potential between two neon atoms. 








Potential, ergs X 10“ 
Experimental 
n=9 n=12 


Theo- 


RX108 cm retical 





124 289 331 
8.9 10.9 31.9 
—0.35 —0.47 0.00 
—0.33 —0.35 —(),2 
—0.21 —0.18 —0,1 








The exponential equation above is of the form 
proposed by Born and Mayer” for the repulsion 
between alkali and halide ions, the value of p 
for all the alkali halides being 0.345 x 10-° cm, 
and the comparable ro values for Na* ion and F- 
ion being respectively 0.875X10-% cm and 
1.110X10-° cm.!* Presumably the more diffuse 
electron distribution function of the negative 
ion determines the steepness (the value of p) in 
the repulsion between the ions. 


APPENDIX 


The analytical expressions of the R;;~' integrals for the 2p electrons of like atoms are given here 
since they have not appeared previously in the literature in closed form. It is assumed that the 
normalized 2p radial function of both atoms is mre-*" where m and k are constants. The integrals 
were calculated by using Neumann’s'‘ expansion in elliptic coordinates of the reciprocal distance 
1/ri2 between two points (1, g1, $1) and (p2, g2, ¢2). This formula for the case p;> f» is 


a $F (102 +1] ah 
R ie fF" 


112 t=0 v=-T 


|| 


)! 2 
= P(p2)Qr(b1) P1v(Gi) Prv(qe)e2-9, 


where P,, and Q,, are the associated Legendre polynomials of the first and second kinds, respectively. 
For the case pi <2 exchange the réles of p; and pz. 


12M. Born and J. E. Mayer, Zeits. f. Physik 75, 1 (1932). 


1M. Huggins and J. E. Mayer, J. Chem. Phys. 1, 643 (1933). 
4 F, Neumann, Vorlesungen tiber die Theorie des Potentials, p. 341, Leipzig (1887). 
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The results are 


-1 


Roo, = 48{147W[1(0, 5, —18, 5);0( )]+15W[1(21, —25, 27, —7); 6(7, —6, 7)] 
+ W[3(35, —30, 3, 0); 10(21, -11, 0) }}, 


Raisi=g{588WT1(0, 5, —6, 1);0( )]+15W[1(21, —25, 3, 1);6(7, —6, —1)] 
+ W[3(35, —65, 33, —3); 10(21, —32, 11) ]}, 





Roi = —140g{9WT1(0, 5, —6, 1); 2(0, 5, —2)]+W[3(7, —10, 3, 0); 2(21, —23, 5) }}, 


Ruix1=70gW[3(7, —15, 9, —1);2(21, —38, 9)], 


where 
g=3R*(m*e—*/53,760a)’, a=kR, 


W[e1(de, @4, G2, do); C2(D5, bs, b1) |=c1{ M2(C+1n a) —2L Me**Ei( —2a)+ L%e*“Ei( —4a) 
4 5 
—(1/2a) }) o2MAj(a)+y;J+(1/2a) ¥ 5;A;(2a)} 
7=0 i=0 


—¢2{MN+(1/2a@)[(5d544+3b3u2+b1y0) — (6a6r5+404r3+ 2023) |} 

L=a¢A 6( — a) +44A 4(— a) +424 2(— a) +a, 

M =a¢A¢(a) +44A 4(a) +42A 2(a) +0, 

N=);A;(a)+63A3(a)+51Ai(a), 
Mf =a¢A 64i(2a) +05A 44:(2@) +024 245(2a) +a0A ;(2a) 
Ni =b5A 54i(2a) +03A 34:(2a)+0:A 14:(2a), 
wi= Mi +ipi-s/a, bo= My’, 
v3= Nj’ +1r;-1/a, vo= No’, 
An, n(a)=An(a)—mla"-"A,(a)/n!,  (OKngm—1) 


and 








Sn=5S,(L, M) => dei[ MAo:, n(a@)+(—1)**-"LA2i, n(—a@) ], t= (2n+3+(-—1)")/4 
6,=5,(1, 1). 


The A ;(a) are the polynomials of Eq. (23) and Ei(x) is the logarithmic-integral function."® C denotes 
the Euler constant 0.5772157. 





* For definition and tables of this function see E. Jahnke and F. Emde, Tables of Functions with Formulae and Curves, 
Second edition, p. 78, B. G. Teubner, Leipzig (1933). 
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On the Dispersion of Sound. Consideration of Three Energy Levels 
Morris E. Rose, Department of Physics, University of Michigan 
(Received February 2, 1934) 


$1. 


In a recent issue of this Journal there ap- 
peared a paper by W. T. Richards! on the dis- 
persion of sound due to vibrational energy lag 
in a gas with several energy states in which the 
cases of two, three and four vibrational levels 
were treated. It is important to realize that to 
consider n levels including the ground state, and 
to neglect all higher levels is equivalent to 
neglecting all powers of the Boltzmann factor 
exp (—e/kT) higher than the (w—1) st, ¢€ being 
the energy difference between two adjacent 
states. Considering the case of three energy 
levels, the calculation of the velocity of sound 
must be carried out consistently in powers of this 
parameter of smallness so that in the final result 
no terms of higher order than the second are 
present. This point, it seems, has not been 
sufficiently recognized in the literature. In 
Richards’ case this has apparently led to an 
improper choice of reaction equations. 


§2. 


Case of two energy levels. The dispersion of 
sound for the case of two levels has been carried 
out by Kneser? and Rutgers.* The method due to 
Einstein‘ is essentially as follows. The reaction to 
be considered is 


(a) NotNomNot Mi, 


where No, N; refer to a molecule in the ground 
state and first excited state, respectively. This 
leads to the reaction equation 


V(dn;/dt) =x no? — Kanon, (1) 


in which V is the total volume, «;, x2 the reaction, 
rate constants characterizing the rate of energy 
adjustment for the forward and reverse reactions, 


1W. T. Richards, J. Chem. Phys. 1, 863 (1933). 
2H. O. Kneser, Ann. d. Physik 11, 761 (1931). 
3A. J. Rutgers, Ann. d. Physik 16, 350 (1933). 
4A. Einstein, Sitz. Ber. Akad. 380 (1920). 
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respectively, and mo, m, the total number of 
molecules in the states 0 and 1, respectively. 
By the principle of detailed balancing the two 
parameters «i, ke are not independent but are 
related by 
ki/kg=exp (—e/kT) =u. 


The calculation of the velocity of sound f may 
be based on the characteristic acoustical equation 


f?=real part (p/p)(1—ApV/pAV), — (2) 


p and p referring to the pressure and density, 
respectively, the displaced equation of state 


ApV =nkAT (3) 


and the thermodynamic equation for adiabatic 
change 
—pAV=nC,AT+«An, (4) 


in which C,, is that part of the total specific heat 
per molecule which is independent of the fre- 
quency, and 1 is the total number of molecules in 
the gas. 

Taking the time dependence of m, about its 
equilibrium value to be of the form exp iwt, « 
being the frequency of the sound wave in radians 
per second, the variation of (1) yields 


tw VAn= komo| MoweAT / RT? — (uw +1) An, | 
= xon(npeAT/kRT?—Any), 


where use has been made of the fact that 
no+n,=n= constant, and therefore Ano= — Am; 
further m= m(1+4). 

Eliminating Am, from (4) and (5) and com- 
bining the result for AT/AV thus found with (2) 
and (3), the final expression for the frequency 
dependent velocity of sound is found to be 


p k é 1 
flw)="(1+ a 1). (6) 
p C, Ca2T? 1+? V?/k2?n? 








The result when expressed in this forra puts in 
evidence the relative order of magnitude of the 
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contribution to the sound velocity due to 
vibrational transitions and the part independent 
of such transitions. The presence of powers of u 
only up to the first is an expression of the fact 
that only one excited level is considered. This 
seems to the writer to be the consequent form in 
which to write the result of Kneser and Rutgers.® 
If the vibrational specific heat is calculated 

from 
C;= (d/dT)(eu/1+y), (7) 


then the expected result in the limiting case of 
small frequency 


F?(0) = (p/p) A +k/C+ Ci) (6’) 
when expanded in powers of yu is 
f?(0) = (p/p)(1+R/C— He /C.°T?) (6) 


in agreement with (6). Similarly from (6) for w 
very large 


f?( 2) = (p/p)(1+k/C.), 
as is to be expected. 


§3. 


Case of three energy levels. The importance 
of an order of magnitude consideration in terms 
of the parameter » becomes more evident when 
we consider transitions to two excited states for it 
is this factor which is instrumental in determining 
the choice of reactions to be considered. The 
reactions to be selected are of the type N;+N; 
2N,+N). Now the probability of such a reaction 
taking place is determined by two factors: (1) the 
probability of finding a molecule in state 7 and 
another in state 7; (2) the probability that this 
pair of molecules have sufficient relative kinetic 
energy (plus vibrational energy) to effect the 
transition. Both these probabilities depend on 
the Boltzmann factor uw. With this in mind, we 
find beside the first order reaction (a) the follow- 
ing second order reactions to be considered: 


b) Not MoeNit+Mi, 
Cc) Not M224, 
d) Not+Mi1=2No+Nz2, (8) 
e) Not NomNot Nz, 
f) Not N22Ni4+ Mi, 


*The further transformations by which these authors 
ntroduce explicitly into (6) the C;, see (7), are inconsistent 


( 
( 
( 
( 
( 





the last being a ‘‘resonance’”’ reaction.°® 

Assuming the principle of detailed balancing 
for each of the six reactions (8) and (a) we can 
write two reaction equations for dn,/dt and 
dn2/dt analogous to (1), in which six independent 
reaction rate constants appear. Variation of these 
two reaction equations together with (4) modified 
by the additional term (e€;+¢€2)Am2 yields three 
linear equations in An,, Ang and AT, from which 
AT/AV may be found. ¢; and e refer to the 
energy differences between the levels 0-1 and 1-2 
respectively. The resulting expression for the 
velocity of sound is then of the form 


Pap tayw*?+aow! 
fP(w) =fe+ufer+— “pw, (9) 
pP Bot Biw? + Bow*+ B3w® 





where fo?+u/:’ is just the expression (6). The a’s 
and #’s are rather involved combinations of T, p, 
€1, €2 and the reaction rate parameters which need 
not be given here. 

The limit for large frequency remains, of 
course, unchanged. In the limit of small fre- 
quency (9) becomes just the expression (6), 
with n=, plus the second order term 


(p/p) Let/RC8T*+ (1/C,2T*) 
X (2€:? — (€1+€2)?m2/ mi) Jur®. (9’) 


Here wi:=exp (—e:/kT). The result (9’) agrees 


with the above-mentioned order of magnitude considera- 
tions. Their result, however, when expanded is identical 
with (6) to terms of first order in yu. 

6 Richards has considered the reactions N;+N@Nj+N 
where N is No, Ni, Ne successively, i=0, 1, j=1, 2. 
Written out in detail these are 


(a) M+M=2Ni+No, 
(b) M+M2Ni+Mi, 
(c) No+N2=2Ni4+ No, 


(d) Nit No N2+No, 
(e) NitM, N2+M, 
(f) NitN2=2N2+N2, 


so that one could expect to find in his results unwarranted 
contributions from (c), (e) and (f) while contributions 
from (8) (b), (e) and (f) will be absent. Reactions (c) and 
(e) are third order and (f) is a fourth order reaction. 
Richards has further assumed the reaction rate constants 
of each group of three reactions (a), (b), (c) and (d), (e) (f) 
to have separately a common value. This assumption 
seems doubtful; moreover a corresponding one does not 
seem possible in the case of reactions (8) and (a) due to 
the absence of symmetry in these as contrasted to the 
reactions of Richards. Therefore it appears necessary to 
introduce different rate constants for each reaction. 
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with the expansion in powers of u of the expected 
result (6’) in which now 


d eyuit(ert+e2) mime 
"aT 





1+yit wipe 


Richards’ result for this limiting case when 
expanded similarly, with u:= ye, is in accord with 
(9’). However this is to be expected since at low 
frequencies when the gas is in equilibrium, the 
particular reactions taking place can have no 
influence on the velocity of sound. 

If we take all the x’s, except x; the resonance 
rate constant, to be of the order of the experi- 
mental value of ks, 3.67X10~™, (Kneser), then 
for any reasonable value of x; the terms in (9) will 
decrease rapidly with higher powers of w so that 
we may neglect all but the first term in the 
numerator and the first two terms of the 
denominator. The first and second order contri- 
butions may then be combined to give a fre- 
quency dependence for f? of the Kneser-Rutgers 
form with, however, another constant x’ replacing 
ke of (6). We can then form an estimate of the 
second order effect by comparing x’ with xo. We 
find 

x? = «2(1—y)/1—y(2?/KP +2) 
vy =pu(e/RC.T?—2). 


The center of the dispersion region (inflection 
point) is shifted by an amount 


Aw = (x’ —x2)n/3)V. 


Assuming for x; the extreme values x;~1, ks~xa 
numerical calculation for CO, at T=300°A 
yields for these two cases «’~3.85X10-" and 
4.04 10-, respectively, corresponding to a shift 


E. ROSE 


Aw~0.28 10° and 0.58105 sec.—!, and repre- 
senting deviations of 5 and 10 percent. For 
intermediate values of x;, x’ and Aw will assume 
intermediate values. The relative change of the 
velocity of sound in the dispersion region becomes 
of the order of present experimental accuracy so 
that at higher than room temperatures in CO, 
the second excited level should be considered.’ 

The author wishes to express his indebtedness 
and thanks to Professor G. E. Uhlenbeck for his 
advice and interest in the calculation. 


7As Dr. Richards points out in the following note, 
which he so kindly permitted me to see before publication, 
there is in addition to the above considerations a further 
factor instrumental in determining the @ priori probability 
of the reactions. Because of the finite mass ratios involved 
and the particular directional character required of the 
collision in order that excitation take place, the effective- 
ness of the translational-vibrational energy transfer will be 
decidedly limited. If we suppose that relative kinetic 
energy Xe is necessary to excite vibrational energy e, the 
probabilities of the six reactions (a) and (8) become of the 
order p*, pw, uit, wit, Ww and uw? respectively. The 
parameter A, which enters in the x«’s implicitly, may be 
calculated from impulse considerations using the experi- 
mental value of x2. The result indicates that \ may easily 
be greater than 2, in which case the reactions (8) (b) and 
(e) may be justifiably omitted. This can be taken into 
account in the frequency dependent expression by setting 
the corresponding rate constants equal to zero. The a’s 
and @’s, however, are not sufficiently simplified to permit 
of their interpretation, the frequency dependence re- 
maining as in (9), which, as pointed out in the text, is 
only approximately of the form (6). None of the numerical 
results given above will be affected appreciably. If we 
exclude all reactions involving excitation of more than 
one vibrational quantum, p» may still be used as a develop- 
ment parameter since the probability of second and third 
order reactions, for example, (~u'*, 2+, respectively) 
differ by a factor yu as before. 
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Comment on the Paper by M. E. Rose, Entitled ‘On the Dispersion of Sound. 
Consideration of Three Energy Levels’”’ 


WILLIAM T. RICHARDs, Department of Chemistry, Princeton University 
(Received February 26, 1934) 


HERE can be no question that the ultimate 

second order theory of the dispersion of 
sound should be written only in terms of trans- 
ition probabilities which are of the first and 
second orders, and Mr. Rose has found a 
compact and logical method for presenting this 
viewpoint. It is clear, therefore, that in a gas in 
which only the 0, 1 and 2 states need be con- 
sidered collisions of the type Ne+ Ne will be so 
rare per second per unit volume that their 
effect on the rate of adjustment of the vibrational 
energy will be negligible in comparison with 
encounters of the type V,+N;. The appearance 
of third and fourth order reactions in my 
treatment of the second order theory is, as Mr. 
Rose has pointed out, illogical in this sense, but 
was required by the simplifying assumption that 
the collision of an N; molecule with any other 
molecule of the same chemical composition 
having a given relative translational energy could 
be equally effective in bringing about the 
transition V;—>N;. This assumption was the only 
one I could find which gave a result of sufficient 
simplicity to permit its comparison with con- 
temporary measurements. It would be surprising, 
however, if it were justified, and hence it is very 
probable that the expressions which I have given 
must ultimately be modified in order to describe 
the behavior of any actual gas. 

I cannot agree with Mr. Rose, however, that 
the Boltzmann factor is alone sufficient to 
determine a priori the probability of a given 
transition. In all cases which have been investi- 
gated to date a considerable activation energy of 
collision has been found necessary for the 


conversion of vibrational into relative trans- 
lational energy. There are, moreover, reasons of a 
theoretical character for believing that only a 
portion of the relative translational energy may 
become converted into vibrational energy on 
collision. When this is the case transitions of the 
type No+No—Ne2+ No become far less probable 
than those of the type No+Ni>N,+No2, al- 
though the order of probability of the two is the 
same according to Mr. Rose’s criterion. Con- 
versely the reaction No+N2—Ni+ Ni is more 
probable than any other “‘second order’’ process. 
The rate of this reaction can hardly be detected 
by acoustical measurements, however, since, in 
the ideal case of an harmonic oscillator, it 
involves no conversion of translational to vi- 
brational energy, and, in any actual case, in- 
volves only a very small amount of conversion. 

It appears at present that experiment may 
adequately be described, at least in many cases, 
by very simple expressions. If, as has been found 
with carbon disulfide and ethylene, groups of 
states corresponding to several vibrational modes 
leave the acoustical cycle with a common 
characteristic transition frequency, the original 
treatment given by Einstein may be retained 
almost without alteration. Reactions of the type 
No+N2—Ni+M:1 may materially assist in re- 
distributing the vibrational energy, but this may 
be obtained from the rotational-translational 
energy by a single outstandingly efficient collision 
process. It will then be only the identification of 
the rate-determining process by acoustical means 
which presents serious difficulties. 
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The Raman Spectra of Some Compounds Containing Carbon-Bromine Bonds. 
Possible Changes of the Vibration Frequencies in Ionic Solutions 


HERVEY H. Voce, Department of Chemistry, University of California 
(Received March 6, 1934) 


Raman spectra of methyl bromoacetate and methyl phenylbromoacetate are reported. The 
carbon-bromine bond line of the former is similar to that for ethyl bromide, while in the latter 
it does not appear in its normal position at all. The spectra of methyl alcohol and ethy! bromide 
have been accurately redetermined; some minor additions result. Shifts of the methyl bromo- 
acetate and methyl alcohol lines in solutions of bromides are found. 





INETIC data suggest that in certain chemi- 

cal systeins at least a portion of the reacting 
molecules are associated. Association should 
affect the Raman spectra, either producing new 
lines or shifting old ones. An attempt to obtain 
such evidence of complex formation in reactions 
of the Finkelstein type, whose mechanisms have 
been discussed by Olson! and by Bergman, 
Polanyi and Szabo,? was made but the results 
were negative. However, as most of the com- 
pounds and solutions studied have not been 
previously observed, it is thought worth while to 
report the results. 


EXPERIMENTAL 


Spectra were observed with a prism spectro- 
graph of dispersion 34 cm™!/mm at 4358A and 
54 cm™'/mm at 5025A. Three Pyrex mercury 
arcs with appropriate reflectors served as the 
source of illumination. Use of saturated aqueous 
NaNO, or 4 percent m-dinitrobenzene in benzene 
as a filter restricted excitation to 4358 in most 
cases. The frequencies of all lines were obtained 
by linear interpolation between wave numbers of 
the nearest lines of an iron arc, photographed on 
the same plate. This procedure is liable to some 
error. Therefore the weak mercury lines appear- 
ing on the plates were measured in a manner 
identical to that used for the Raman lines, and a 
comparison with the known frequencies gave the 
desired correction. This correction ranged from 0 


1 Olson, J. Chem. Phys. 1, 418 (1933). 
2 Bergman, Polanyi and Szabo, Zeits. f. physik. Chemie 
B20, 169 (1933). 


to 2 cm”!. The error in the values listed depends 
upon the nature of the individual lines, but is 
probably less than 1 cm™ for all definite lines. 


ETHYL BROMIDE 


Eastman Kodak Company’s product was 
refractionated from P,O;. An exposure of 25 
hours with slit 0.04 mm was made. The lines 
obtained were: 290.5 (4), 560.3 (10 b), 959.4 
(2 b), 1062.8 (2), 1244.7 (2), 1252.9 (2 b), 1444.0 
(2 b), 2867.2 (2), 2924.4 (4), 2969.2 (2 bb or d). 
Here 6=broad, d=double. All these lines have 
been observed before* but the line at 1245 has 
been resolved into a doublet for the first time. 

A 1:1 mole mixture of ethyl bromide and 
methyl alcohol showed only a superposition of 
the lines of the two components. 

In asolution of anhydrous lithium bromide and 
ethyl bromide in methyl alcohol having the 
composition (mole fract.) LiBr 0.072, C.H;Br 
0.301, CH;0OH 0.627, no new lines appeared and 
it was found that the 290, 560, 959, 1062, 1445 
and 2924 lines of ethyl bromide were unshifted 
within the limits +1 cm™. 


METHYL BROMOACETATE 
The Raman spectrum of methyl bromoacetate 
is given in Table I. It has not been observed 
before.* 


3Cleeton and Dufford, Phys. Rev. 37, 362 (1931); 
Dadieu and Kohlrausch, Monatshefte 55, 70 (1930). 

* Added in proof: The Raman spectrum of methy! 
bromoacetate has recently been published by Cheng (Hua- 
Chih Cheng, Zeits. f. physik. Chemie B24, 293 (1934))- 
The agreement with the values listed above is fairly good. 
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TABLE I. CH,BrCOOCH; liquid. Prepared from bromoacetic, 
acid and methyl alcohol by the method of Lepercq.4 








Twice redistilled at 10-20 mm. Slit 0.07 mm. Exposure 
48 hours. 





205.8 (2) 845.9 (2) 
278.4 (2) 882.4 (3) 
301.9 (2) 909.9 (2) 1454.3 (5) 
374.1 (4 6 and 1003.2 (46) 1395 to (1 band) 

band to 420) 1116.4 (2) 1465 
550.3 (10 db) 1153.4 (2 6) 1739.7 (7) and band to 1762 
670.2 (7) 1189.6 ? (3) 2845.9 (2) 
710.5 (7) 1215.5 (2 b) 2956.9 (5) 
766.5 (1 b) 1287.8 (4) 2965 to 3050 (2 band) 


1404.8 (2) 
1421.0 (3) 








A solution of methyl bromoacetate and lithium 
bromide in methyl alcohol with the composition 
(mole fractions) LiBr 0.123, CH2BrCOOCH; 
0.215, CH;OH 0.662 gave the following results. 
No new lines were observed. The lines 278, 374, 
550, 710, 845, 882 and 1739 were unshifted 
within the limits +1 cm™, 670 increased by 2 
cm™', 910 appeared to decrease 3.5, but as it is 
quite broad this is not certain. 


METHYL PHENYLBROMOACETATE 


The ester was prepared from phenylbromo- 
acetic acid (melting point 81°-82°. Kindly 
furnished by Miss M. J. Young) and methyl 
alcohol in the presence of sulfuric acid. Ex- 
posures of 15 to 48 hours were made on methyl 
alcohol solutions with slit widths of 0.05 to 0.07 
mm. The lines obtained, averages over four 
exposures, are given in Table II. bd=band. 


TABLE II. CsH;CHBrCOOCH; in CH;0H. 








616.2 (3) 
675+15 (1 bd) 
707.8 (1 db) 
782.1 (3) 
846.9 (2 b) 
898.5 (2) 


1003.3 (7) 
1156.3 (2 b) 
1184.5 4 bb) 
1218.3 (4) 
1602.4 (7) 
3064.7 (3) 








Notable in the spectrum of methyl phenyl- 
bromoacetate is the absence of the C-Br line 
from its usual position at 550, where it has been 
found for all other aliphatic bromo compounds.® 

In view of the thorough investigation of 
vibration spectra of carbonyl compounds now 


*Lepercq, Bull. Soc. Chim. (3) 11, 297 (1894). 
*Dadieu, Pongratz and Kohlrausch, Monatshefte 61, 
409 (1932). 


being carried out by Kohlrausch and his co- 
workers no attempt had been made to interpret 
the spectra of methyl bromoacetate and methyl 
phenylbromoacetate, observed for the first time 
in this work. 


METHYL ALCOHOL 


Methyl alcohol, used as a solvent, was rather 
carefully investigated. The interpretation given 
its spectrum is shown in Table III together with 
the results of others. The alcohol used was a 
commercial product twice fractionated with 
retention of only the center third. No lines of the 
common impurities were observed. The values 
listed are means from three plates with a slit of 
0.05 mm and exposures to 41 hours. 


TABLE III. CH30H. 








Kohlrausch in 
Der Smekal-Raman 
Effekt® 


Observed Collins? 





1032.3 (6 b) 
1062.7 ? (1) 
1109.8 (1 d) 
1166 (1 b) 


1034 (5) 1028 (10) 
1107 (2) 
1362 (0) 


1462 (5 db) 


1452.2 (4 b) and 
1430 to 1490 band 
center = 1470 


2834.4 (5 Db) 
2919 (2 side band) 
2943.7 (4 b) 
2980 (2 side band) 


1457 (10) 


2729 (1) 


2835 (6 b) 2831 (10) 


2943 (5 b) 2942 (10) 


3388+85 (band) 








In an exposure at low dispersion 1360 was 
positively identified. 1063 is doubtful as it may 
possibly have been 2834 excited by the 4047A 
mercury line. The unsymmetric band at 1460 has 
been rather arbitrarily interpreted as a combi- 
nation of line and band. 


TABLE IV. Shifts of the CH;OH lines. 








13 mol 
% HBr 


Solu- 7 mol 
tion % LiBr 





Line 
1032.3 
2834.4 
2943.7 








6 Kohlrausch, Der Smekal-Raman Effekt, p. 309. 
7 Wood and Collins, Phys. Rev. 42, 386 (1932). 
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The methyl alcohol lines were slightly dis- 
placed in concentrated LiBr and HBr solutions. 
The changes are shown in Table IV and agree in 
direction with those observed by Hibben® for 
ZnCl, solutions, and by Goubeau® for LiClO, 
solutions. The alterations of the methyl alcohol 
lines are roughly proportional to concentration. 
They indicate a weakening of the C—O bond 
(1032 line) and a strengthening of the C—H bond 
(2900 lines). LiBr and HBr form complexes with 


8 Hibben, Proc. Nat. Acad. Sci. 18, 523 (1932). 
9 Goubeau, Naturwiss. 21, 468 (1933). 


HERVEY H. VOGE 


methyl alcohol’: " and it is to be expected that 
the intermolecular forces involved will effect the 
binding forces. The weakening of the C—O bond 
is in accord with the important part played by 
the oxygen in compound formation.” ¥ 

In conclusion I wish to thank Professor A. R. 
Olson for suggesting this topic. 


10 Turner and Bissett, J. Chem. Soc. 105, 1777 (1914). 

11 Maas and McIntosh, J. Am. Chem. Soc. 34, 1284 
(1912). 

12 McIntosh, J. Am. Chem. Soc. 28, 588 (1906). 

13 Archibald and McIntosh, J. Chem. Soc. 85, 919 (1904). 
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The Para-Ortho Hydrogen Conversion by the Hydrogen Iodide Reaction 
and by Iodine Atoms 


E. J. ROSENBAUM AND T. R. HoGness, Department of Chemistry, University of Chicago 
(Received March 5, 1934) 


The rate of conversion of para to ortho hydrogen was determined in the presence of equi- 
librium concentrations of hydrogen iodide and iodine. From these data, velocity constants 
and collision efficiencies for the conversion by iodine atoms were calculated. The results are 
in agreement with a theoretical equation derived by Wigner. 





N the basis of experiment there is no 

evidence of any appreciable difference be- 
tween the velocity constant of a _ reaction 
involving ortho hydrogen and that of the corre- 
sponding reaction involving the para modifica- 
tion. Therefore, in any reversible reaction 
involving hydrogen in which the hydrogen is not 
in equilibrium with respect to its two modifi- 
cations, the para and ortho modifications are 
removed by the reaction in the proportion 
existing in the reacting mixture, while the 
hydrogen formed by the reverse reaction is 
composed of para and ortho hydrogen in the 
equilibrium ratio, which is 1 : 3 if the reaction is 
taking place at a temperature above 250°A. If 
the hydrogen present at the start of a reaction 
contains a fraction of para hydrogen larger than 
the equilibrium proportion of }, the rate at 
which the fraction of para hydrogen changes can 
be related to the velocities of the forward and 
reverse reactions. 

A method for determining the velocity con- 
stants of a reaction at equilibrium is thus made 
possible, at least in principle. A reaction suitable 
for experiment must satisfy several conditions: 
the equilibrium concentration of hydrogen must 
be large enough to allow the withdrawal of 
samples for analyses; the reaction must proceed 
with a moderate velocity in a temperature range 
below 750°A, where a thermal reaction becomes 
important; and the conversion reaction must not 
be catalyzed directly by any component of the 
reacting mixture. 

At the inception of this investigation the 
hydrogen iodide reaction seemed an excellent 


example which satisfied these conditions. More- 
over, reliable data on the velocity constants of 
this reaction were available and a simple and 
reasonable mechanism had been proposed to 
explain these data. 

When the results of the experiments here 
reported had been analyzed, it was found that the 
values of the velocity constants of the hydrogen 
iodide reaction which were obtained were so 
much larger than those which had been obtained 
by Bodenstein! that it was necessary to look for 
an additional conversion mechanism. Experi- 
ments showed that the hydrogen iodide and 
iodine molecules had no catalytic effect. It was 
known that hydrogen atoms could effect the 
conversion of para to ortho hydrogen, but the 
equilibrium concentrations of hydrogen atoms at 
the temperatures at which the experiments had 
been performed were far too small to account for 
the large velocity constants. Similarly, any 
possible catalytic effect of the walls of the 
reaction tubes was not large enough to account 
for the observed rate of conversion. 

The iodine atoms, the only remaining com- 
ponent of the reacting mixture, were tentatively 
considered to be responsible for the catalysis. 
This consideration was greatly strengthened by 
the appearance of a series of papers which 
showed, both theoretically and experimentally, 
that paramagnetic solids, gases and solutions are 
effective in catalyzing the para-ortho hydrogen 


1 Bodenstein, Zeits. f. physik. Chemie 13, 56 (1894); 
22, 1 (1897); 29, 295 (1899). 
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conversion.? Since iodine atoms are para- consider a mechanism involving two processes in 
magnetic, it seemed reasonable to ascribe to accounting for the observed rate of conversion of 
them the additional rate of conversion beyond para to ortho hydrogen: 
that calculated from Bodenstein’s data. The 
experimental results are consistent with this H2’+1,>2HI>H2’ +1, 
assumption. and 

With this assumption, it is necessary to H,?+1->H,°+1. 





The rate at which para hydrogen is removed by both processes is 
—d[H.")/dt = (K,/2)[H2" [Ie] +K-[H2 [1]. 
K, as used here is the velocity constant of the formation of hydrogen iodide: 
dL HI ]/dt = —2d[H: ]/dt. 
The rate at which para hydrogen is formed is 
d(H.” ]/dt = (K,'/8)[HI }?+ K.’[H," ][1]. (2) 


The factor } appears in the first term of the right side of Eq. (2) because only } of all the hydrogen 
molecules formed are of the para modification. The net rate at which para hydrogen is removed is 


— d(H»? ]/dt = (Ki/2)(He” (12 ]— (Ki'/8) CAT P+ (K2[H2” ]— Ke’(H2")) (1). (3) 


When this is divided by the total concentration of hydrogen, and the following relationships are 
introduced : 


CHI} Ky K,’ 1 


———_—_ =—— and —= 
[He j{I2] Ky’ K.+K 4 


(since Kz=3K,’ above room temperature), then Eq. (3) becomes 
—dN,/dt=(K,/2)[12](Np—4)+4Ke'(N,—4)L1], (4) 

in which N, is the mol fraction of para hydrogen. When Eq. (4) is integrated, there is obtained 
Np=it(Np,—4) exp { —((Ki/2)[T2]+4Ke'[1])¢}. (5) 


With N, and ¢ measured, N,, known from the temperature of the charcoal catalyst, K, obtained 
from Bodenstein’s data, and [I] calculated from [I.] by the use of the equilibrium constant of the 
reaction I,—2I, this equation can be solved for Ke, the velocity constant for the para to ortho hydro- 


gen conversion: 
3 (In(4N,-1) Ky, 
— +—[I,] . (6) 
4[T] t 2 


At the highest temperature used an additional process, due to the equilibrium concentration of 
hydrogen atoms, is effective in converting para to ortho hydrogen: 


H+H,.?>H+H.’; —dN,/dt=K;3(N,—}3)(H]. 





The inclusion of this process leads to the integrated equation 

N,=i+(N>,—4) exp {—((Ki/2)(12]+4K2'[1]+K;s[H])é}. (7) 
The values of K;[H] were determined from blank experiments in which 1 : 1 para-ortho hydrogen 
only was heated. 


? Farkas and Sachsse, Zeits. f. physik. Chemie B23, 1 (1933). Wigner, ibid. 28 (1933). Taylor and Diamond, J. 
Am. Chem. Soc. 55, 2613 (1933). Taylor and Eyring, Proc. Am. Phil. Soc. 72, 255 (1933). 





PARA-ORTHO 


EXPERIMENTAL 


Mixtures of equilibrium proportions of 1:1 
para-ortho hydrogen, hydrogen iodide and iodine 
were sealed into Pyrex reaction tubes, heated at a 
constant temperature for a measured length of 
time, and then rapidly cooled. The hydrogen 
iodide and iodine were frozen out with liquid 
nitrogen, each tube was opened and the hydrogen 
was analyzed for its ortho-para composition in a 
thermal conductivity cell. 

Commercial electrolytic hydrogen was purified 
by passing it through a heated palladium tube, F 
in Fig. 1. C.P. iodine was resublimed and stored 
over phosphorus pentoxide. To obtain dry 
hydrogen iodide, a colorless concentrated solution 
of hydrogen iodide was dropped on phosphorus 
pentoxide in an evacuated generator J, and the 
gas was passed through a phosphorus pentoxide 
trap. 

The form of the thermal conductivity cell 
adopted was essentially that used by other 
investigators. The principal difference lay in the 
fact that the spreader for the heated wire was 
held by a glass rod attached at the top of the cell 
to the electrode seal. The resistance of the cell 
wire (a 0.01 mm Wollaston wire) was measured 
by a very sensitive Wheatstone bridge. The 
wires used had resistances in the range 80-115 
ohms at operating temperature. 

The furnace was of the resistance type, with a 
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Fic. 1. Vacuum line. A is the thermal conductivity cell, 

isa mercury manometer, C is a mercury piston, D is the 
activated charcoal catalyst, E is a trap, F is a palladium 
tube, G is a Toepler pump, H is an interchangeable ground 
glass joint and J is the hydrogen iodide generator. 
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large heat capacity. An iron jacket surrounded 
the reaction tube. Because of its high thermal 
conductivity this metal lining made the tempera- 
ture of the reaction zone of the furnace uniform 
to within 1°. The temperature was measured 
with a calibrated chromel-alumel thermocouple 
and a potentiometer. It was possible to hold the 
temperature constant to within 1°. 

The form and dimensions of the reaction tubes 
were chosen to give a volume large enough for 
duplicate analyses of the hydrogen, together 
with a small length to insure uniformity of 
temperature throughout the tube. A vacuum- 
tight connection was made between each of the 
reaction tubes and the vacuum line by means of 
an interchangeable ground joint H. The neck of 
each tube was closed off by a thin glass bubble 
which could be broken by releasing a weight 
with a small solenoid. A piece of capillary tubing 
in the side arm which was used to fill each 
reaction tube made negligible the effect of 
diffusion to the heated area when the tubes were 
sealed off. 

The calculated equilibrium amounts of iodine 
were contained in capsules with capillary tips. 
To prepare one of these capsules an open 
capillary tip was drawn on a piece of 2 mm 
tubing. After this piece of tubing was weighed, an 
approximate amount of iodine was introduced 
into it and it was purged thoroughly with dry 
hydrogen while the capillary was kept in contact 
with solid carbon dioxide to prevent excessive 
loss of iodine. With the hydrogen flowing, the 
capillary and the other end of the capsule were 
successively sealed a few centimeters from the 
iodine. The capsule was weighed together with 
the rest of the original piece of tubing, and the 
weight of the iodine was obtained by difference. 
A capsule was sealed into each reaction tube 
(except those used for blanks); the tubes were 
sealed on the line and were evacuated to a 
pressure of 10-* mm or lower. 

A 1:1 para-ortho hydrogen mixture was 
obtained by placing liquid nitrogen around the 
tube D which contained degassed activated 
charcoal, and allowing the charcoal to adsorb 
hydrogen until the pressure above the charcoal 
was about 20 cm. After an interval of an hour, 
some of the hydrogen was pumped away; then 
the reaction tubes were filled to a desired pressure 
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of 1:1 para-ortho hydrogen by means of the 
Toepler pump G. The stopcocks on the tubes 
were closed, the hydrogen remaining in the line 
was pumped out, and hydrogen iodide was added 
to each tube until the desired total pressure was 
reached. An effort was made to have each stop- 
cock attached to a reaction tube open for as short 
a time as possible to prevent diffusion of hydro- 
gen. The reaction tubes were sealed off at a 
previously prepared constriction and the capil- 
laries of the iodine capsules were broken by 
shaking the tubes. The tubes were then ready to 
be placed in the furnace. 

After a tube was heated for a measured length 
of time, it was removed from the furnace and 
cooled as quickly as possible in the air current of 
an electric fan. The grind on the tube was fitted 
into the one on the line H and the line was 
evacuated while the reaction tube was immersed 
in liquid nitrogen to condense the iodine and 
hydrogen iodide. The breaker was released to 
open the tube and a sample of hydrogen was 
pumped with the Toepler pump into the thermal 
conductivity cell A through the trap E immersed 
in liquid nitrogen to prevent traces of iodine and 
hydrogen iodide from reaching the cell. The 
pressure of hydrogen in the cell was adjusted 
with a mercury piston to a value of approxi- 
mately 40 mm, which was reproducible to 0.02 
mm. The thermal conductivity cell was always 
covered with liquid nitrogen, the level of which 
was kept constant. After the current had been 
regulated, the resistance of the cell was measured. 


RESULTS 


With the final composition of the hydrogen N, 
determined this value could be substituted in 
Eq. (6), and since the other quantities were 
known, Ke could be calculated. By using the 
average value of Ke for each temperature, the 
mol fraction of para hydrogen could be recalcu- 
lated and compared with the experimental value. 
The data obtained are compiled in Table I, from 
which the agreement between the calculated and 
observed values can be seen. 

No changes were noticed in the blank experi- 
ments carried out at a temperature of 716°A or 
lower. 

The errors in the pressures of hydrogen and 
hydrogen iodide were less than 0.5 percent. The 
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TABLE I. 








I; Time Np 
(moles/cc) 107 (sec.) Calc. Obs. 





Temperature 753°A 
Mean constant 9960 cc/mol sec. 


2.54 7860 0.289 
3.33 7500 .281 
1.95 5120 .340 
1.14 4380 .383 
1.88 5280 339 
2.84 4500 329 
3.84 6500 .284 
0.91 5480 .376 

Blank 7200 

Blank 7200 

Blank 7200 


Temperature 716°A 
Mean constant 3690 cc/mol sec. 
7020 0.407 
7220 Re 
7770 343 
8080 374 
9200 .348 
7730 343 
9800 358 
6260 409 
7770 343 


_ 
Now 


DEMS SKS Oe 
WIR Rae wick 


_ 


_ 
Cw 


Temperature 693°A 
Mean constant 2600 cc/mol sec. 
9660 0.424 
13610 391 
5770 448 
8030 428 








volumes of the reaction tubes were measured to 
within 1 percent. The iodine was weighed by 
difference; the weighing error varied from 0.3 to 
4 percent. An error in the concentration of iodine 
did not cause an equal error in the observed value 
of N,, however, because of the logarithmic 
relationship between the iodine concentration 
and N, and the fact that the equilibrium 
concentration of iodine is established during the 
heating period. The intense heat required to seal 
off the tubes caused a slight amount of con- 
version which was reduced to a negligible value 
by the use of capillary tubing in the side arm of 
the reaction tubes. The temperatures of the 
reaction tubes were within one degree of the 
reported temperatures. 

A possible source of error is the catalytic effect 
of the Pyrex walls of the reaction tubes. H. A. 
Taylor® found a relatively small variation in the 
decomposition velocity constant with surface- 
volume ratio, using crushed glass to increase the 


3H. A. Taylor, J. Phys. Chem. 28, 984 (1924). 
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surface. In the present work, however, the 
surface-volume ratio was smaller than any in 
Taylor’s work, approximately 1.4 cm™, and the 
walls were smooth. That the condition of the 
wall had no marked effect on the reaction 
velocity was shown by a comparison between 
reaction tubes which had not been cleaned 
carefully and those which had been cleaned with 
hot K,Cr,O;—H»SO, solution just before sealing 
off. 

DIscussION 


It can be seen that the use of the mean velocity 
constant leads to calculated values of N, which 
are in agreement with the data. The consistency 
of the data can also be judged from the deviation 
from the mean of the individual values of the 
velocity constant which is less than 15 percent 
with the exception of two values, one of which 
deviates by 30 percent, the other by 25 percent. 
In view of the fact that each determination 
involves a considerable number of manipulations 
and that the desired constant is obtained as a 
difference, this degree of concordance is as high as 
could be expected. 

From the values of the velocity constants, the 
collision efficiency can be calculated for each 
temperature. In particular, the collision efficiency 
for the para to ortho hydrogen conversion can be 
obtained by dividing the product of the corre- 
sponding velocity constant, the concentration of 
para hydrogen molecules and the concentration 
of iodine atoms by the total number of collisions 
between para hydrogen molecules and iodine 
atoms. The calculation of this total number of 
collisions depends on an assumed value for the 
closest distance of approach between a hydrogen 
molecule and an iodine atom, since this quantity 
has not been determined. With an assumed 
distance of approach of 2A, the collision efficiency 
for the para to ortho hydrogen conversion is 
found to be 9.410-”, 13 10-" and 3410-” 
at 693, 716 and 753°A, respectively. 

It is of interest to compare these values with 
some obtained from a theoretical equation of 
Wigner* and with experimental data obtained by 
Farkas and Sachsse* for other reactions. By an 
application of quantum mechanics, Wigner has 
derived expressions for the probability of a para 





‘See reference 2. 
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to ortho conversion in those collisions in which 
the required energy is present: 


we j+i 
P 54541 1.1«10-" — ere 
a’T 2j7+1 


2 


M 
P34j-1=1.1K10-" 
a’T 2j7+1 








in which yw is the magnetic moment of the 
paramagnetic particle (in this case the iodine 
atom), a is the effective minimum distance of 
approach, 7 is the absolute temperature, 7 is the 
rotational quantum number of the hydrogen 
molecule, and the factor 1.1 x 10~'” contains only 
universal constants and the square of the 
magnetic moment of the proton. 

The normal state of the iodine atom is *P3,/.. Its 
magnetic moment, calculated from the expression 
u=([g?7(j+1) ]! (g, the Lande splitting factor) is 
2.58 Bohr units or 2.38X10-* gauss cm. It is 
unnecessary here to consider any excited state of 
the iodine atom, since the fraction of atoms in the 
first excited state, ?Pj,2, is only 310-7 at 750°A. 

In the complete absence of any information 
about the quantity a, we may assume, with 
Farkas and Sachsse, a value of 1A. It should be 
observed that this distance need not coincide 
with the kinetic theory distance of closest 
approach, usually obtained from viscosity data, 
which was used to calculate collision efficiencies 
from experimental data. 

When these values are introduced into 
Wigner’s expressions and the effective value of 
the factor involving j is obtained by summing 
over five® rotational states of the hydrogen 
molecule, the probability, and thus the collision 
efficiency, of para to ortho transitions is found to 
be 6.2X10-", 6.0 10-" and 5.8X10-" at 693, 
716 and 753°A, respectively. 

It is evident that so far as the order of magni- 
tude of the collision efficiency is concerned, the 
theory is in agreement with the experimental 
data. However, the theoretical equation of 
Wigner requires a very small negative tempera- 
ture coefficient, while the temperature coefficient 
of collision efficiency obtained experimentally has 


a small positive value. It is not surprising that the 


5 The fraction of molecules with a j value of 5 is only 


0.035. 
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agreement is less good at the higher temperatures 
since the theory was designed to fit a lower 
temperature range. 

The results of Farkas and Sachsse for a 
number of paramagnetic substances are also in 
agreement with Wigner’s equation. In particular, 
the order of magnitude they obtained for the 
collision efficiency of NO(2X10-") is approxi- 
mately the same as that found in the present 
work with iodine atoms. They, too, found a 
positive temperature coefficient in the tempera- 
ture range 600—750°A. 

It might be worth while to mention at this 
point a confirmation of a result obtained by a 
previous investigator. A. Farkas,® in the study 
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of the thermal conversion of para to ortho 
hydrogen, determined values for the velocity 
constant for this reaction which have a small 
temperature coefficient and from which the 
extrapolated value at 750°A is approximately 
9x10" cc/mol sec. From the blank experiments 
at 753°A, K;{[H]=2.6X10-*% The hydrogen 
atom concentration at this temperature calcu- 
lated from the equilibrium constant is 2.3 10-" 
mols/cc. From these data the value of K; is 


(2.6 X 10-5) /(2.3 X 10-7) = 11 X10" 


which confirms the values obtained by Farkas. 


6 A. Farkas, Zeits. f. physik. Chemie B10, 419 (1930). 
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By using the principle of microscopic reversibility, 
unimolecular decompositions are considered from the point 
of view of the reverse reaction, which is a bimolecular 
association. In cases where the immediate product of the 
decomposition is a pair of free radicals, it is assumed that 
the probability of these free radicals recombining upon 
collision is independent of their energy. An explicit 
expression for the activation energy in terms of the energy 
of reaction and the thermal properties of the molecules 
involved is thus obtained. From this one may calculate 


the energy of activation at absolute zero and the energy 
necessary to break the bond which falls apart in the 
reaction, quantities which may differ considerably from 
each other and from the energy of activation at the 
temperature of reaction. The results are applied to several 
cases of interest, and are critically compared with those of 
the earlier theories of Rice and Ramsperger, Kassel and 
Rice. A brief discussion is given of the calculation of the 
rate constant as a function of pressure, and its modification 
according to the point of view of the present paper. 





I. INTRODUCTION 


F the theories that have been put forth to 

explain the experimental data on homo- 
geneous first order gas reactions, those that give 
best agreement with the data are theory II of 
Rice and Ramsperger, the classical and quantum 
theories of Kassel and the quantum theory of 
Rice. All these involve arbitrary assumptions 
about the nature and mechanism of dissociation. 
It is assumed that a necessary (but not sufficient) 
condition for reaction is that the total energy of 
the molecule be equal to or greater than a con- 
stant quantity €9. A second necessary condition 
is that this energy be localized in the molecule. 
According to the theory of Kassel, the energy 
must be localized in a bond (two energy terms), 
that is, the bond that breaks. This permits an 
identification of .¢9 with the energy of the bond.‘ 
It can be shown that the translational energy of 


1Presented, in part, at the Boston meeting of the 
American Association for the Advancement of Science, 
December, 1933. 

* Thayer Fellow. 

3 (a) O. K. Rice and Ramsperger, J. Am. Chem. Soc. 40, 
1617 (1928). (b) Kassel, J. Phys. Chem. 32, 225, 1065 
(1928). (c) O. K. Rice, Proc. Nat. Acad. Sci. 14, 114, 118 
(1928). For a more detailed analysis of these theories 
than will be given in this paper and for a more complete 
list of references to original papers, see (d) Kassel, Kinetics 
of Homogeneous Gas Reactions, Chem. Cat. Co., 1932, 
especially Chapter V. 

* More exact definition of the various energy quantities 
will be given below. 


the molecule does not contribute to the energy 
necessary for dissociation, while usually rota- 
tional energy contributes a trivial amount, so 
that only those molecules that have vibrational 
energy greater than or equal to ¢€ 9 can react. 
Sometimes, to get agreement with experiment, 
it is necessary to assume that only a few of the 
vibrational degrees of freedom take part. 

A consequence of the classical theories of Rice 
and Ramsperger and of Kassel, and of the 
quantum theory of Kassel, as well, is that for 
high pressures the energy of activation, as de- 
termined by the temperature coefficient of the 
rate constant, is equal (or at least very approxi- 
mately equal) to €9; and it has been more or less 
customary to consider that the energy of activa- 
tion gives an experimental determination of €». 
According to the quantum theory of Rice, 
however, €9 may be very considerably less than 
the energy of activation, and hitherto there has 
been no satisfactory way of determining which 
of the quantum theories gives the most nearly 
correct extension of the classical theories.® 

The difference between the theories is due to 
different arbitrary assumptions made in calcu- 
lating the probability that in a molecule having 
energy €>€ 9, €9 or more will be localized in a 
given part of the molecule. It would, therefore, 


5Q. K. Rice, L’ Activation et la Structure des Molecules, 
Réunion Internationale de Chimie Physique, Paris, 1928, 
p. 309. Kassel, reference 3 (d), pp. 105-6. 
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be of great advantage if we could replace these 
assumptions by some of a more rational char- 
acter. We believe that we have been able to do 
this by making use of the principle of microscopic 
reversibility and considering the reverse reaction, 
for decompositions in which the reverse reaction 
has no activation energy or in which the activa- 
tion energy of the reverse reaction consists only 
of the mutual translational energy of the re- 
combining parts; and on the basis of more 
careful definitions of the various quantities in- 
volved we have considered the calculation of the 
energy of the bond and the energy of activation 
at absolute zero from the observed activation 
energy. 

Recently, with a different goal in view, 
Rodebush* has given a quasithermodynamniic deri- 
vation of the rate of reaction for a simple type of 
unimolecular reaction by a somewhat similar use 
of the reverse reaction. In this paper we shall 
go more deeply into the statistics involved in the 
treatment of the direct and reverse reactions, 
and shall at first restrict our discussion to the 
energy of activation at high pressures, where the 
Maxwell-Boltzmann distribution is maintained. 


II. GENERAL CONSIDERATIONS 


According to the principle of microscopic 
reversibility the reaction in equilibrium with a 
unimolecular decomposition is a bimolecular 
association.’ Let us consider the reaction AB 
=A +B. At equilibrium the number of molecules 
AB that are dissociating is equal to the number 
of molecules AB that are being formed by the 
bimolecular association of A and B. Furthermore, 
the average energy of the molecules AB that are 
dissociating is equal to the average energy of 
those being formed. Therefore we have the 
relation 


éaptean =taten° testes +€ (1) 


where éapt EAB’; éa ten’, and és ten? are re- 
spectively the average energies of the molecules 


6 Rodebush, J. Chem. Phys. 1, 440 (1933). La Mer (J. 
Chem. Phys. 1, 289 (1933)) has also considered energy of 
activation as a function of temperature but has applied 
his results only to reactions in solution. 

7The possibility of bimolecular association has been 
discussed by Kassel, J. Am. Chem. Soc. 53, 2143 (1931). 
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of AB, A and B that actually react,’ and € is 
the energy of the bond that is formed.’ The 
quantities €4,°, «4° and e,° are the respective 
zero point vibrational energies, while the other 
symbol designates the average energy of those 
molecules which react, in excess of the zero 
point energy. 

We wish now to consider the activation energy 
of the direct and reverse reactions. Activation 
energy has been defined in many ways by 
previous investigators.!° For cases in which the 
Maxwell-Boltzmann quota of activated molecules 
is maintained, i.e., the cases in which equilibrium 
statistics can be applied, the most rigorous 
definition of activation energy is that given by 
Tolman,® as the difference between the average 
energy of the molecules that react and the 
average energy of all the molecules. Therefore for 
the activation energy of a unimolecular reaction 
we have 

€a = EaB— EAB, (2) 


where é4g is the average energy of all the 


8 See Tolman, Statistical Mechanics, pp. 261-2, Chem. 
Cat. Co., 1927. Our definitions of é4z8, €4, and ég are the 
same as Tolman’s definition of the corresponding quanti- 
ties.* But see reference 11. 

All average energies include the translational energy. 

*In previous published work the symbol ¢« with two 
bars above it has been used for this quantity. Throughout 
this paper the symbol é is being substituted. 

9 The exact significance of the quantity € is perhaps 
best understood by considering the corresponding quantity 
for diatomic molecules, namely, the energy from the 
minimum of the lowest potential energy curve representing 
the bond to its asymptote (if potential energy curves 
cross, to the lowest asymptote). It is therefore identical 
with the quantity that spectroscopists call D, (Jevons, 
Report on Band Spectra, Chap. 1X, Cambridge University 
Press, 1932). It must be distinguished from Do, the 
quantity that is experimentally determined spectroscopi- 
cally, and which differs from it by the amount of the zero 
point energy. In this paper we shall call © the energy of 
the bond, and distinguish it from AEo, the energy change 
of the reaction at absolute zero, which gives the difference 
between the energy of the associated molecule in its lowest 
energy state and the energy of the dissociated particles in 
their lowest energy states. The difference between © and 
AE, will be given by the difference between the zero point 
energy of the associated molecule and the sum of the zero 
point energies of the dissociated particles. In the literature 
the concept of bonding energy for the case of polyatomic 
molecules has not been unequivocally defined. 

10 See for example Rice and Ramsperger, J. Am. Chem. 
Soc. 49, 1618 (1927); and Kassel, reference 3 (d), p. 21- 
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molecules AB in excess of their zero point 
energy. Similarly, for the bimolecular reaction 
the activation energy will be given by" 


én’ = ta tép—Ea— Ep. (3) 
Combining Eggs. (1), (2) and (3), 
€a=€atep—Eapten +E+er°t+en°—enn®. (4) 


If now we let AE be the energy of reaction and 
AE, the energy of reaction at absolute zero,” 
then according to the usual definitions of these 
quantities 


€a= Ea tEn—Eapten +ALs, 
€a = €q' + AE. 


(Sa) 
(Sb) 


We shall need a special symbol for the activa- 
tion energies of the unimolecular and bimolecular 
reactions at absolute zero; they will be denoted 
respectively by €a, 9 and e€’a, o. 

In many recombinations that involve heavy 


Jn the case of the bimolecular reaction our definition 
differs slightly from that of Tolman. He defines the energy 
of activation as the average energy of all pairs of reacting 
molecules minus the average energy of all pairs of colliding 
molecules. This is not very clearly brought out in his book, 
but is clearly stated by Kassel, reference 3 (d), p. 26. 
We define the activation energy as the average energy of 
all pairs of reacting molecules minus the average energy of 
all pairs of molecules. The latter quantity is 3k7 less than 
the average energy of the colliding pairs (i.e. we have an 
“energy of activation of collision’), so that in this way 
we avoid the term 37 which appears in Tolman’s Eq. (614). 
If we write Na, Ng, and Naz for the concentrations of the 
respective molecules and define K’ by the equation 
dNap/dt=K’NaNp, where dN4z/dt represents the rate of 
increase of Ng due to reaction between A and B, we have 
simply 

kT°d log K’'/dT =é4+ép—€4— Ez. 
If we combine this with the equation 
kT*d log K/dT =tap—€ap 


for the unimolecular reaction, it is readily seen from the 
immediately following development of this paper that we 
get the correction thermodynamic relation for the temper- 
ature dependence of the equilibrium constant K’/K. It 
must be noted, however, that if the bimolecular rate 
constants are defined in terms of partial pressures rather 
than concentrations a slight correction must be made. 
Such a correction has in fact been applied in the case of 
the data on the association of ethylene and hydrogen 
discussed below. 

“For thermodynamic quantities we use the notation 


of Lewis and Randall, Thermodynamics (McGraw-Hill, 
1923), 
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REACTIONS 


free radicals it seems very probable that no 
activation energy will be required (we neglect 
the ‘‘energy of activation of collision’’"). This is 
the case that will be considered first. For this 
special case Eq. (5b) reduces to 


é,=AE (6a) 
and we also have 


&., = AE. (6b) 


The equations of this section are not new in 
themselves,"* but the terms which enter into 
them do not seem to have been carefully defined ; 
when this is done their range of applicability is 
much increased. 


III. Discussion OF THE ACTIVATION ENERGY IN 
SEVERAL SPECIAL CASES 


The derivation of Eqs. (6a) and (6b) is almost 
on a thermodynamic basis. Nevertheless, it may 
not be entirely superfluous to mention several 
cases in which it may be shown that they are 
substantiated by the experimental facts. 

An example in which (6b) is applicable is 
afforded by measurements recently made in this 
laboratory by Dr. R. A. Ogg, Jr., on the 
unimolecular decomposition of methyl iodide 
into methyl radicals and iodine atoms.’ The 
experimentally determined activation energy is 
equal to AE, for the C—I bond as calculated from 
thermochemical data. In this case the change of 
the vibrational energy (indeed, the whole vibra- 
tional energy) will be negligible. 

All the information necessary to test (5b) is 
supplied by data available on the reaction 
C:H,@C.H,+He. The decomposition of ethane 
has been studied by Marek and McCluer” who 
find it to be a first order reaction with an 
activation energy of 73.0 kg cal. per mole. The 
association has been studied by Pease'® and 


13 For example, they have been used by Bodenstein and 
Jung (Zeits. f. physik. Chemie 121, 129 (1926)), who got 
them from Trautz (Lehrbuch d. Chemie 1, 406 (1922); 
3, 81 (1926), Berlin). The relation used here was not 
exact being Q=«.— ea’ where Q is the heat of reaction. 

4R, A. Ogg, Jr., J. Am. Chem. Soc. 56, 526 (1934). 

18 Marek and McCluer, Ind. Eng. Chem. 23, 878 (1931). 

16 Pease, J. Am. Chem. Soc. 54, 1876 (1932). We have 
made our calculations for Pease’ temperature, about 
800°K. Marek and McCluer worked at slightly higher 
temperatures, but this difference can be neglected. 
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found to be bimolecular with an activation 
energy of 44.75 kg cal. per mole. (This has been 
corrected so that it applies to the rate constant 
expressed in terms of concentrations rather than 
pressures.'!) The heat of hydrogenation at room 
temperature (AH3o.) has been determined by 
v. Wartenburg and Krause!’ as 30.6 kg cal. per 
mole. AC, as a function of temperature is given 
by the equation 6.5—0.0047 cal. per mole per 
degree.'® If we make the approximation AC, 
=AC,—R we may replace AE by AH-—RT. 
We have therefore 


€q, s00 = AH go — RT + €’ a, 800 


800 
=AHawt f AC,dT —R-800+ ' a, 800- 
300 


=76.0 kg cal. per mole, 


the subscripts denoting the absolute tempera- 
tures at which the quantities are evaluated. 
This compares favorably with Marek and 
McCluer’s value of 73.0 kg cal. 

Pease has already made a calculation to show 
that the ratio of the reverse and direct reaction 
rates gives the equilibrium constant, which in- 
volves essentially the above calculation, except 
that it does not take account of the AC, term. 
This term is small in this case, involving as it 
does only small molecules. Since (5b) differs 
from equations in common use only in that 
these terms involving the vibrational and other 
internal energy are explicitly taken into account, 
it would be of more interest if it could be verified 
for a case in which these are large, but data for 
this purpose are not available. 

The real importance of the equations of Section 
II lies in the fact that they give us the possibility 
of calculating © and e., 9 in various cases and 
comparing with the quantum theories of Kassel 
and Rice. This comparison is, however, rendered 
difficult by the fact that the zero point energy 
was not carefully considered in the earlier work. 
It was treated as a more or less non-distributable 
and constant residue of energy. On this basis it 
would seem that the quantity called Ey by 


17y, Wartenburg and Krause, Zeits. f. physik. Chemie 
151, 105 (1930). 

18 Parks and Huffman, The Free Energies of Some Organic 
Compounds, p. 80, Chem. Cat. Co., 1932. 
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Kassel®» and ¢) by Rice** corresponds more 
nearly to our €, 9 than to our &. Accepting this 
identification, then according to the theory of 
Kassel €., 9 should always be just equal to the 
activation energy at the temperature at which it 
is measured, while according to the theory of 
Rice it should be less than the energy of activa- 
tion at the temperature of reaction and in- 
creasingly so for larger molecules. We may thus 
consider €, 9 in attempting to form an opinion 
as to which of the older theories is the better 
approximation, but © would seem to be in- 
trinsically the more interesting quantity, and we 
shall consider how to calculate both of them 
from €q. 

The calculation of «€,, 9 and © rests upon our 
ability to find out something about the reverse 
reaction. In cases where dissociation is into free 
radicals it seems reasonable to assume that the 
bimolecular association requires no activation 
energy,'® which implies that the steric factor is 
essentially independent of the energy. This very 
important assumption will be shown in the next 
section to take the place of the assumption made 
in the previous theories of unimolecular de- 
composition that the rate of reaction of a 
molecule of a given energy is proportional to the 
fraction of such molecules which at equilibrium 
would have sufficient energy localized in a given 
bond. This latter assumption was purely ad hoc; 
the one which replaces it seems very natural, 
especially if a large number of degrees of freedom 
are involved in the reaction; for at the energy 
at which a molecule is activated the energy 
levels are fairly closely and uniformly spaced in 
an actual molecule.2° The probability of two 


19 A consideration of this assumption, and how our 
results would be affected if it does not hold, will be given 
later on in the paper. We should remark that there is some 
evidence that free methyl radicals do not combine in the 
gas phase (Ogg, reference 14; Bates and Spence, J. Am. 
Chem. Soc. 53, 1689 (1931)). Since the recombination 
does not occur over a considerable temperature range it is 
unlikely that need of activation energy is the cause but 
rather some quantum-mechanical selection rule concerning 
rotation. Work is now being done on this problem. 

20 To verify this point the spectrum of the vibrational 
energy levels of the cyanogen molecule was plotted as far 
as the level representing 27.43 kg cal. per mole. Between 
26.01 kg cal. and 27.43 kg cal. there are ten separate energy 
levels. As the total energy of the molecule increases the 
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radicals recombining also depends upon the 
amount the energy level, into which they are to 
go, is broadened. But all the energy levels are 
presumably equally broadened, since the broad- 
ening depends upon the lifetime of the molecule 
in the given state, and the latter quantity de- 
pends only upon the collision frequency, when 
the pressure is sufficiently high. 

In order to see how our assumption works out 
it is desirable to treat a specific case which is as 
simple and exact as possible. To avoid arbitrary 
assignation of vibrational frequencies an actual 
molecule with known frequencies is preferable. 
For this purpose the cyanogen molecule was 
chosen. Although there are no data on the 
kinetics of the reaction, the equilibrium C.Ne 
@2CN had been studied by G. B. Kistiakowsky 
and one of us,”' and the vibrational frequencies 
are known. Hence it is possible to assume that 
the decomposition is unimolecular, and to calcu- 
late €,, 9 and © in terms of ¢,. Since we assume 
é =0, it follows that the average energy of the 
radicals that combine upon collision is equal to 
the average energy of all the radicals. We can, 
therefore, set up our equilibrium statistics and 
calculate the energy distribution at equilibrium. 

The molecules CN and C.Ne each have, of 
course, three translational degrees of freedom, 
and since the molecules are linear they have only 
two degrees of rotational freedom. CN radicals 
have one vibrational frequency, while C.Ne2 has 
seven normal modes of vibration, only five of 
which are distinct, as due to the symmetry of 
the molecule the transverse vibrations are double. 
We thus get, for the average energies of CN and 
C:No, respectively, 


€on = (5/2)RT+hvo/Lexp (hvo/kT)—1] (7) 


where v9 is the vibration frequency of CN, and 


foone = (5/2)kT+¥ hy./Lexp (hv;/kT)—1] (8) 


where the »; are the vibrational frequencies of 


_ 


energy levels become closer together. Since the energy of 
dissociation for this molecule is about 77 kg cal. the energy 
levels in this region will be exceedingly close together. 
xe also, Kassel, reference 7. 

* Kistiakowsky and Gershinowitz, J. Chem. Phys. 1, 
432 (1932). 
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C.Ne. If we now insert the necessary numerical 
values of the vibrational frequencies” and let 
T =1500°K, the temperature at which the equi- 
librium was studied, we obtain 


(9a) 
(9b) 


€cn= 8.47 kg cal. per mole 


€c,n, = 19.87 kg cal. per mole. 
Hence 


€x, o= AEg=€,+2.9 kg cal. per mole. (10) 


The zero point energy of cyanogen is 11.2 kg cal. 
per mole (being 1/2 hv for each oscillator) and 


that of cyanide radicals 2.9 kg cal./mole, so 
using Eqs. (4) and (10), we get 


= e,+8.3 kg cal. per mole. (11) 


Considering the high temperature involved it is 
seen that the effect of the internal energy is not 
large in this case which deals with a rather small 
molecule, though the difference between © and €, 
is nevertheless appreciable. For molecules of this 
small size neither the theory of Rice nor that of 
Kassel would demand much difference between 
€.4,0 and ¢€. Since, however, such deviation as 
exists is in the opposite direction to that pre- 
dicted by the theory of Rice it would appear 
that as far as may be indicated by this specific 
case the theory of Kassel is to be favored. 


Before considering the application of the equations of 
Section II to more complicated molecules it is desirable to 
examine more closely the relation between the concepts of 
“binding energy’ and energy of dissociation of organic 
compounds as used in the literature. The calculations of 
Pauling*? have shown that to each bond in an organic 
compound can be assigned an energy practically inde- 
pendent of the compound, such that the sum of all the 
“binding energies” gives the energy of formation of the 
compound from its component atoms. In the application 
of these binding energies to chemical data it has been 
assumed (at least implicitly) that they should be equal to 
the energy of breaking the molecule into two parts at 
that bond. Thus in order to explain the low dissociation 
energies of the highly substituted ethanes it has been 
necessary to assume that the bond is weakened, in the 
cases where phenyl groups are the substituents by reso- 
nance,*4 and in others by steric hindrance; more careful 
investigation, however, indicates that there is no reason 


22 Reference 21, p. 436. 
23 Pauling, J. Am. Chem. Soc. 54, 3570 (1932). 
* Pauling and Wheland, J. Chem. Phys. 1, 362 (1933). 
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to expect the dissociation energy to be equal to the bonding 
energy. The quantity measured in chemical reactions is 
AE, the difference between energies of the products and 
the reactant, both with their zero point energies. Even 
though the energy of all the carbon-carbon bonds in 
heptane were the same we should expect that at absolute 
zero the energy required to break off an end CH; group 
would be different from that required to break the chain 
in the middle, simply because the products have different 
zero point energies in each case. In calculating energies of 
formation, however, the zero point energy is of a so much 
smaller magnitude than the total energy of formation that 
it can be neglected. But when we consider ordinary 
chemical decompositions, involving only one bond and 
try to compare the energy of dissociation with a so-called 
“binding energy”’ it is first necessary to investigate the 
AC, change in the reaction, as well as the change in zero 
point energy. 

Thus there are two distinct quantities that have been 
called binding energy in the literature, one of which 
corresponds to our AX», the other to ©. The quantity 
calculated by quantum mechanics,” at least for diatomic 
molecules, is ©. This, therefore, is the quantity that 
might be expected to be constant for a given bond in 
various compounds if anything is. On the whole, however, 
there is little evidence either experimental or theoretical 
for the constancy of either quantity. Quantum-mechanical 
calculations by Eyring”* show in fact that the hydrogens 
on ethane are bound less firmly than those on methane. 
A wave-mechanical analysis of methane using group theory 
indicates that it takes more energy to pull off the first 
hydrogen than the rest.?” Heitler and Schuchowitzki,”® 
investigating bonding energy in organic molecules, find no 
way of explaining why it should be constant to more than 
20 percent at best. As we have stated above, these binding 
energies are not the same as those considered by, Pauling.” 
The equations derived in this paper afford a means of 
computing energies of dissociation and bonding energies 
from activation energies, provided the mechanism of the 
reaction and something about the properties of the 
reactants and products is known. 


Let us now consider the decomposition of the 
azo compounds. Azomethane,?** methylazoiso- 
propane?*» and _ isopropylazoisopropane?** are 


25 We refer to calculations by the variation method 
which manipulate the electronic wave functions holding 
the distance between nuclei fixed. 

*6 Eyring, J. Am. Chem. Soc. 54, 3202 (1932). 

2? Eyring, Frost and Turkevich, J. Chem. Phys. 1, 783 
(1933); also see Van Vleck, ibid. 2, 24 (1934); Mulliken, 
ibid. 1, 500 (1933). 

28 Heitler and Schuchowitzki, Phys. Zeits. d. Sow. 3, 
241 (1933). 

2? Ramsperger, (a) J. Am. Chem. Soc. 49, 912, 1495 
(1927); (b) 51, 2134 (1929); (c) 50, 714 (1928). The 
energies of activation used are those calculated by Kassel, 
reference 3 (d). 
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homologues that decompose unimolecularly with 
activation energies of 52,440 cal., 47,480 cal. and 
40,900 cal., respectively. Recent investigations 
have shown that the most probable mechanism 
for these reactions involves the breaking of a 
C—N bond with the formation of free radicals,*” 
so we assume that the reverse reaction requires 
no activation energy and Eqs. (6a, b) are 
applicable. 

Although there are no data on the vibrational 
frequencies of the azo compounds, some estimate 
of the frequencies involved can be made by con- 
sidering the frequencies in compounds containing 
equal numbers of heavy atoms. Examination of 
data obtained from Raman and infrared absorp- 
tion spectra*' indicates that the larger the 
molecule the lower the frequency of the trans- 
verse vibrations, and hence the more temperature 
vibration and the less zero point energy per 
oscillation in the molecule. Since the reaction 
products ‘are smaller than the decomposing 
molecule, it is seen that decomposition will 
necessarily be accompanied by a decrease in 
temperature vibration (more or less counter- 
balanced by the increase in translational and 
rotational energy) and an increase in vibrational 
energy. 

As a rough first approximation we might 
assume that on decomposition all the frequencies 
remain unchanged except four transverse vibra- 
tions, one longitudinal vibration, and one free 
(classical) rotation that are converted into six 
classical degrees of freedom, three of translation 
and three of rotation. For a compound with 
four heavy atoms these frequencies will be on 
the order of 600 cm~ as in cyanogen (or perhaps 
a little lower, since there are triple bonds in 
cyanogen), while for eight heavy atoms they 
may be as low as 200 cm~. Using these figures 
we will have for methylazomethane at 600°K, 


€A +é€p—€spr~0.3 kg cal. per mole, 


€a° +en°—€4n°~ —4.3 kg cal. per mole, (12a) 
€a, p= AEy~ €.—0.3 =52.1 kg cal. per mole, ’ 


€~e,+4.0=56.4 kg cal. per mole, 


39 Leermakers, J. Am. Chem. Soc. 55, 3499 (1933); 
F. O. Rice and Evering, ibid. 55, 3898 (1933). 

31 Andrews, Phys. Rev. 36, 548 (1930); Kohlrausch, 
Der Smekal-Raman Effekt, pp. 304 ff. Springer, 1931. 





ACTIVATION 


and for isopropylazoisopropane 

ea t+é€s—€an~ —1.7 kg cal. per mole, 

ex° + en® — €an°~ —1.4 kg cal. per mole, 

€z, o= AEo~ €a+1.7 = 42.6 kg cal. per mole, 
€~e,+3.1=44.0 kg cal. per mole. 


(12b) 


It is thus seen that AE» and &€ are different 
for the two azo compounds, and methylazoiso- 
propane will be intermediate. It does not seem 
possible that any reasonable assumption as to 
the frequencies could bring the values of © 
together. If, however, there is a somewhat more 
profound alteration of frequencies, i.e., a general 
shifting to higher values, when the molecule 
decomposes (which is not beyond the realm of 
possibility) the values of AE») may be closer 
together, as we shall show. We still assume that 
the longitudinal vibrations of the hydrogens are 
more or less independent of the size of the 
molecule. As to the rest of the vibrations, we 
assume in the case of azomethane that we have 
18 vibrations of average frequency 1000 cm~ 
(which may be taken to include the free rotations 
in this case) going into twelve vibrations of 
average frequency 1100 cm and six classical 
degrees of freedom, which gives us €4+€s—€4z 
~1.3 kg cal. per mole at 600°K. If for isopropyl- 
azoisopropane we assume that we have 46 vibra- 
tions of average frequency 700 cm! going into 
40 vibrations of average frequency 1000 cm~! 
and 6 classical degrees of freedom we get é4+€, 
—€4p~ —6.3 kg cal. per mole. We get for methyl 
azomethane €4°+¢s°—e41n°~—6.8 kg cal. per 
mole and for isopropylazoisopropane €,°+€;° 
—e,n°~11.1 kg cal. per mole; therefore for 
methylazomethane 


to, o=AEo~eg—1.3=51.1 kg cal. per mole 
4 (13a) 
€~e.+5.5 =57.9 kg cal. per mole 


while for isopropylazoisopropane 
&, o= AEo~ €.+6.3 =47.2 kg cal. per mole 


13b) 
€~e,—4.8 = 36.1 kg cal. per mole. 


So with these assumptions we have made AE, 
the same for the two compounds, almost within 
the limits of experimental error, while the 
difference in the © values is much increased. 
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We feel, however, that the assumptions made, 
though by no means out of the question, may be 
a little extreme in the direction of making the 
average frequencies change too much on re- 
action. The assumptions resulting in Eqs. (12a, b) 
may also be too extreme, but tending in the 
opposite direction. We may thus feel reasonably 
certain that the true values for AE, and € lie 
between those we have obtained, provided the 
reaction mechanism we have taken is correct. 

The matter of the reaction mechanism is one 
of some importance, and, unfortunately, one 
that cannot be definitely settled at the present 
time. We might, for example, have the azo 
compounds breaking into three pieces—a nitro- 
gen molecule and two alky! radicals—in one step. 
The reverse reaction would then be trimolecular, 
but in this case the method of calculation would 
not be substantially affected. On the other hand, 
we might suppose one bond to break, as we have 
hitherto assumed, but with a simultaneous re- 
arrangement of some kind in the fragments.” 
Such a rearrangement might well differ as 
between a methyl and an isopropyl compound, 
and one might perhaps invoke something of 
this sort to explain the difference in the & values 
we have found. And if, in addition, some of the 
energy of activation is thrown back into the 
molecule in the process of this rearrangement, 
and left localized at or near the place which 
breaks, then of course our calculations would 
not be justified. For then the recombining 
radicals would, by the principle of microscopic 
reversibility, also require an energy of activation 
and we should be barred from making any very 
effective calculation by the same difficulty which 
arose in connection with the earlier theories, 
namely the difficulty of calculating the fraction 
of molecules with a given total energy which 
have a certain amount of energy localized in 
some particular way. 

Be that as it may it seems reasonable to con- 
clude from the examples we have given that, 
in cases where a bond is broken and other 
complications do not come in, the quantum 
theory of Kassel is to be preferred as a first 


32 Such a rearrangement is postulated in theories of 
chemical binding advanced by Heitler and Schuchowitzki, 
reference 26. 
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approximation over that of Rice, and we may 
judge the order of magnitude of the correction 
(to get €., 0 from ¢.) which is needed, it being 
much smaller than demanded by Rice’s theory. 
Furthermore, in cases where the activation 
energy of the reverse reaction is small or neg- 
ligible, we see how to calculate the important 
quantity ©, which was not considered at all in 
the older theories. 


Though difficulties might arise similar in nature to the 
possibilities discussed in the last paragraph but one, we 
nevertheless believe that it is very probable that consider- 
ations like those applied to the azo compounds may be 
used in discussing the decompositions of hydrocarbons 
into hydrogen and an unsaturated compound, in spite of 
the fact that it obviously does not involve simply the 
breaking of a single bond, but rather the splitting off in 
one step of a hydrogen molecule. In this case, of course, 
the reverse reaction requires an activation energy but this 
is probably largely the relative translational energy 
necessary to get the hydrogen molecule over a potential 
energy hill. It may* also involve an excitation of the 
vibrational energy of the hydrogen, but if it does not 
involve any localization of energy in the unsaturated 
molecule in order for the hydrogen to unite with it, it 
seems very reasonable to suppose that the probability 
that the hydrogen will stick will be very nearly independent 
of the amount of energy the system may have, provided 
only that the relative translational energy and the vibra- 
tional energy of the hydrogen be at least equal to the 
minimum. Since no carbon atoms split off in this reaction, 
it is probable that the average vibrational frequency is in 
this case the same before and after the reaction, and it is 
very probable that éa, 9 and € are both essentially equal 
to the value of «. at the temperature of reaction. Thus in 
this case, Kassel’s quantum theory would be correct. The 
values of e. for the dehydrogenation of propane and 
n-pentane have been determined by Marek and McCluer 
to be 61.8 and 46.5 kg cal. per mole, respectively, which 
are to be compared to 73.0 for ethane. In this case while 
varies from one compound to another, it is seen from the 
heats of combustion, which are known, that AE» is constant 
from compound to compound.* Thus, though the activa- 
tion energies of the forward and reverse processes in the 
hydrogenation of the hydrocarbons differ from one com- 
pound to another, their differences remain constant. Just 
why this should be the case is not at present clear. 

In this connection two series of compounds studied by 
Coffin,* in which the activation energy €a seems to remain 


33 From the data of Kharasch, Bur. Standards J. Re- 
search 2, 359 (1929), the heats of hydrogenation of both 
propylene and butylene can be calculated to be 32 kg cal. 
as compared with the 31 kg cal. experimentally determined 
for ethylene. 

* Coffin, J. Chem. Phys. 2, 48 (1934); Can. J. Research 
6, 417 (1932). 
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constant from compound to compound, are of 1nterest. 
To explain these data, Coffin has adopted the assumption 
made on different grounds by Rice and Ramsperger and 
by Rice*® that a molecule may be divided into reactive 
and nonreactive sections (the latter having little or 
nothing to do with the part of the molecule which actually 
reacts) and he has offered a new definition of activation 
energy as ‘‘the difference between the average energy of 
the reactive sections which react and the average energy 
of all the reactive sections.’’ This definition necessitates 
such a complete separation of the sections that absolutely 
no exchange of energy between reactive and nonreactive 
sections can take place; nor can the various sections affect 
each other’s vibration frequencies. It seems doubtful to us 
whether as strict a separation as this can exist, and the 
parts still be tied together at all. Further, in view of the 
considerations of this paper, it is our feeling that the data 
do not warrant a reformulation of the definition of activa- 
tion energy, for given this new definition it would again be 
necessary to go through an analysis such as has been given 
in Part I of this paper, treating the average energies of 
the reactive parts instead of the whole molecules, and éq, 
and & for the series of reactions would have to be calculated 
before any comparisons could be made. In the cases 
discussed by Coffin, for example the decompositions of 
butylidene diacetate and ethylidene diacetate, the fact 
that the absolute rates are different at the same tempera- 
tures, while the activation energy is the same, indicates 
to him that the respective homologues have different 
numbers of degrees of freedom in their active parts. If 
this were so the average energies of the reactive parts of 
the homologues would be different and it is difficult to see 
offhand why the differences between these quantities and 
the average energies of the activated parts of the molecules 
that actually react should be the same for the homologues. 
For these reasons it is our opinion that the data of Coffin 
should be treated on the same basis as the other reactions 
we have discussed. There is so little basis at present for 
estimating the vibrational frequencies of these molecules 
that we have made no attempt to find ¢€¢, 9 and &, but we 
feel that ultimately this must be done to establish the 
validity of any interpretation of the constancy of the 
observed activation energy. It seems very likely, however, 
in these cases, in which the group of atoms added as one 
goes from one member of a homologous series to the next 
member is small compared with the whole molecule, that 
the values of €a, »>—€. and €—, will not differ much from 
one member of the series to the next, so that the fact that 
€, remains constant is really significant. In fact, it seems 
more likely to be significant if we disregard Coffins 
suggestion that the molecules are so completely divided 
into separate portions, than if we take it into consideration. 
These examples of the application of the equations of 
Section II, will serve to illustrate what we hope may be 
done with much better success when more experimental 
data become available. As our knowledge of vibrational 
frequencies and specific heats of complex molecules grows, 


% (a) Rice and Ramsperger, J. Am. Chem. Soc. 50, 619 
(1928); (b) Rice, Zeits. f. physik. Chemie B7, 226 (1930). 
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the usefulness of our equations will continually increase. 
It is to be noted that in the derivation of Eqs. (6a) and 
(6b) no assumptions have been made regarding the 
mechanism of energy exchange within the molecule or the 
number of degrees of freedom involved in the reaction. 
It is hoped, therefore, that these equations may serve as 
testing blocks for theories of unimolecular reaction based 
upon the statistics of the dissociating molecule, and 
perhaps suggest new approaches to the solution of this 
important problem. 


IV. DEPENDENCE OF THE RATE CONSTANT ON 
THE PRESSURE 


In the foregoing portions of the paper we have 
confined ourselves to considerations involving 
only the average energy of the reacting molecules. 
The theories of Rice and Ramsperger, Kassel and 
Rice were, however, chiefly concerned with de- 
termining the change of the rate constant with 
pressure, and this involves a detailed considera- 
tion of the energy distribution in those molecules 
which react, and the determination of the 
probability b(¢) that a molecule of given energy, 
¢ (or €4z, but we will omit the subscript for 
simplicity), will react. To determine b(e) it was 
necessary, as we have noted, to make some 
rather artificial assumptions; we shall now show 
how the considerations of this paper allow b(e) 
to be determined if recombination of the parts 
in the reverse bimolecular reaction is independent 
of the total energy of those parts. This last 
condition is equivalent to the assumption that 
& is zero (but if it is necessary for the colliding 
particles to have a certain relative kinetic 
energy to get over a potential hill, and once they 
have this energy the probability of their re- 
combining is independent of the excess energy, 
the indicated modification is. readily made). 
With ¢’., 9 equal to zero, if a molecule of energy 
«decomposes (€, of course, does not include the 
zro point energy) it will go over into a pair of 
molecules whose total energy, above their zero 
point energy, is €—€a, 0. Now we let the equi- 
librium fraction of all undissociated molecules 
with energy between ¢ and e+de (i.e., the 
distribution function of the molecules) be W(e)de, 
and we let the corresponding equilibrium fraction 
of such pairs of the colliding dissociated molecules 
a have a combined energy between e—€,, » and 
~€a, o+de be F(e—€a, o)de. Then by the assump- 
tions we have made and the principle of micro- 
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scopic reversibility, we must have 
b(e) W(€) =kiF(e—€a, 0), (14) 


where &; is a constant which contains a steric 
factor and a collision number. From this 0(e) 
may be calculated and the result may then be 
inserted into the general equation of Rice and 
Ramsperger.*® 

This differs from the result of the old theories 
only in the nature of the distribution function 
on the right. According to Kassel’s classical and 
quantum theories W(e—e,,9) would replace 
F(€—€a, 9). This is also approximately true of 
the classical theory of Rice and Ramsperger, 
while the quantum theory of Rice would demand 
quite a different distribution function. 

There is some evidence, however, as mentioned 
in the last section, for the supposition that in 
certain cases only a limited number of degrees 
of freedom are directly concerned in the re- 
action,® that is only these particular degrees of 
freedom transmit energy freely to the bond which 
is to break in the reaction, and the rest of the 
molecule exchanges energy only slowly with this 
which may be called the active part of the 
molecule. At first sight then it might appear 
that the characteristics of the reaction, at 
reasonably high pressures (so long as only a 
negligible portion of the activation of the active 
part is due to the slow transfer of energy from 
the inactive part), would depend only on the 
properties of the active part. That this cannot be 
strictly true, however, is evident from the fact 
that the general equations of this paper contain 
only quantities which refer to the whole molecule. 
That nevertheless no contradiction can arise 
will appear from the following calculation. 

In the case under consideration it will be seen 
that the quantity b(¢) W(e)de will be replaced by 

€ 
def W‘(e—e.) W°(e.)b(e-)de., (15) 
0 
where W’(e.)de, is the probability that the active 
part of the molecule contain an energy between 
€- and ¢e,.+de., while W‘(e—e.)de is the proba- 
bility that the inactive part of the molecule 
have an energy between e—e. and e—e.+de, 
and b(e.) is the probability that the molecule 


react provided the energy of the active part 


% Reference 3a, p. 1622, Eq. (8). 
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is €-. (This last quantity, it will be observed, is 
entirely independent of the condition of the 
inactive part of the molecule which is of course 
the important characteristic of any theory which 
distinguishes between an active and an inactive 
part.) That the above expression is the one which 
should replace W(e)b(e)de is readily seen from 
the fact that W‘(e—.) W*(e,.)de.de can be shown*? 
to be the probability that the active part have 
energy between ¢, and ¢.+de,. and the whole 
molecule have energy between ¢ and e+de. 
The integral will therefore give the entire con- 
tribution to the reaction which comes from 
molecules whose total energy lies between e 
and e+de. The expression (15) therefore may be 
set equal to ki F(e— eq, o)de, and b(e.) may be 
obtained by solution of the integral equation 
thus set up. 

Let us now compare this with what would be 
expected from the simplest of the older theories, 
namely the classical or the quantum theory of 
Kassel, on the assumption that in the calculation 
of b(e.) the inactive part is to be treated as 
though it were non-existent, and the calculation 
made on the assumption that the active part of 
the molecule is the whole thing. According to 
these theories of Kassel we should have on this 
basis 


b(e.) =k2aW°(€-— €a, 0)/W(e), 
= 0, 


if €.>€a, 0 
(16) 
if Ec < Ea, 0 


where kz is a constant. 
Substituting (16) into (15) we get 


kde f W°(€-— €a, 0) Wi(e—€-)de-. (17) 
0 


But since W°(e.-—€a, 0) is zero for €-<é€a, 9 we 
have from (17), if we let €.=€a, o+x, 


hele f  “(e) Wile—x— ea, o)dx. (18) 
0 


This can be seen to be equal to keW(e— a, o)de, 
for the integrand (if we include the factor de) is 
equal*’ to the probability that the active part 
of the molecule have energy between x and 
x-+dx and the whole molecule have energy be- 


37 For the treatment of a similar problem see reference 
3a, p. 1622. 
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tween €—€,, 9 and e—e,, o+de. Therefore, if the 
two theories were to agree, we should have 


koW(€—€a, 0) = ki F(e—éq, 0) (19) 


This shows that the amount by which the results 
of the two theories will differ will depend only 
upon the difference in the energy distribution 
for the direct and the reverse reactions, as 
before, and that it is the entire energy of the 
molecule, not of any part, which comes into 
consideration, and insofar as (19) is a good 
approximation, (16) may also be taken as a 
good approximation. This conclusion has an 
application in forthcoming work by Rice and 
Sickman, for their work is based on an assump- 
tion equivalent to (16). 

The above deductions are, as stated above, 
based on the assumption that activation and 
deactivation of the active part is by collision 
with other molecules, rather than due to the slow 
leakage of energy from or to the inactive part. 
At very low pressures, however, the latter mode 
of activation will be dominant, and may cause a 
second apparently unimolecular stage of the 
reaction to appear at pressures so low that the 
true high pressure unimolecular rate constant 
has already dropped off. This has been discussed 
by Rice*> and will be discussed by Rice and 
Sickman in more detail. What we wish to note 
at this place is that a very peculiar effect comes 
into play in this low pressure region. If two free 
radicals unite, the active part becomes de- 
activated by leakage of energy to the inactive 
part (from which it is eventually removed by 
collision). But if the energy of the two free 
radicals is high we may expect that enough 
energy will find its way back to the bond which 
breaks in the reaction before it is removed from 
the inactive part by collision, so that the 
molecule will decompose again. Such a molecule 
cannot truly be said to have been formed, and 
should not be counted when we reckon the rate 
of recombination. Thus the probability will 
depend, and depend strongly on the energy of 
the colliding molecules, i.e., k: will not be 
independent of ¢, and the considerations of the 
preceding paragraphs will not hold. Thus we can 
draw no conclusion here about the temperature 
coefficient of the low pressure unimolecular 
stage. 
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Determination of the Speed of Flames and the Temperature Distribution in a 
Spherical Bomb from Time-Pressure Explosion Records! 


BERNARD LEWIs? AND GUENTHER VON ELBE,’ U. S. Bureau of Mines 
(Received February 2, 1934) 


A method has been developed for determining the 
velocity of flame relative to the mass movement of the 
gases, in a closed spherical bomb from an analysis of the 
time-pressure record of the explosion. The speed of the 
flame can be evaluated at any moment during its progress 
from the center to the periphery of the bomb, as well as 
the temperature existing in the unburned phase, the 
temperature immediately behind the flame front, the 
temperature gradient from the latter point to the center 
of the bomb, and the pressure in the bomb at the same 
moment. Given a certain fraction burned of the total 
amount of gas, the volume occupied by the products can 
be determined for three conditions: (1) Before it has 
expanded against the rest of the unburned gas; (2) after 
it has expanded; and (3) when combustion is complete 
and it has been compressed by subsequent burning of gas 
nearer the periphery. Calculations of flame speeds, temper- 


atures, etc., have been made and tabulated for explosions 
of mixtures of ozone and oxygen. The speed of flame 
increases from the center of the bomb to the wall. At the 
same time the pressure and temperature of the gas about 
to be burned increases. The temperature gradient in the 
gas from the center to the wall has been calculated when 
the combustion is complete. A complete diagrammatic 
description is given for one explosion. It is shown that the 
temperature gradient actually existing in the bomb does 
not affect the specific heat results obtained by the usual 
method of calculating the final temperature from the 
maximum pressure by means of the gas law. It is pointed 
out that the practical coincidence of the pressure curve 
with the zero line in the first part of the time-pressure 
record is not due to a time lag between passage of the 
spark and ignition, but to the small fraction of gas which 
has burned during this time. 





INTRODUCTION 


METHOD of measuring the velocity of 

flames was successfully developed by 
Stevens.* It consists of igniting a combustible 
gas mixture at the center of a soap bubble and 
photographing the spherical flame as it expands 
at constant pressure. The method is simple and 
direct. It eliminates the serious disadvantages 
of flame-speed measurements in tubes, but it 
suffers from certain limitations. First, it cannot 
be used for gases that attack or are affected by 
a soap film. Second, the gases are necessarily 
always contaminated with water vapor. 

On the other hand, time-pressure diagrams 
obtained from explosions in a closed spherical 
vessel can be analyzed successfully for the 
determination of flame speeds. The advantages 


‘ Published by permission of the Director, U. S. Bureau 
of Mines. (Not subject to copyright.) 

* Physical chemist, Pittsburgh Experiment Station, 
U. S. Bureau of Mines, Pittsburgh, Pa. 

* Member of staff, Coal Research Laboratory, Carnegie 
Institute of Technology, Pittsburgh, Pa. 

*F. W. Stevens, J. Am. Chem. Soc. 48, 1896 (1926). 


of the soap-bubble method are preserved, while 
its disadvantages are eliminated. Experimentally, 
the method is both simple and applicable in 
general. It will be apparent from the following 
that a single explosion record provides a whole 
series of data on the explosive mixture under 
consideration, from which can be determined 
the flame speeds at different temperatures and 
pressures, the specific heats of the burned® and 
unburned gases* separately, the temperature of 
the unburned gases, and the distribution of 
temperature in the burned gases at any time 
during the explosion up to and including the 
time the flame reaches the wall. 

The experimental procedure has been de- 
scribed in detail elsewhere.’ Briefly, it consists 
of measuring the increase in pressure during the 
explosion of a gas mixture ignited at the center 
of a closed spherical bomb. The change of 
pressure with time was recorded on a rotating 
film by the movement of a light point resulting 


5 Lewis and von Elbe, J. Am. Chem. Soc. 55, 511 (1933). 
6 See following paper, p. 291. 
7 Lewis and von Elbe, J. Am. Chem. Soc. 55, 504 (1933). 
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from the deflection of a specially constructed 
diaphragm manometer. 

When a gas mixture is ignited at the center of 
the bomb the fiame spreads spherically and 
concentrically toward the wall. Since the flame 
movement is slow compared to the velocity of 
sound, pressure equilibrium is considered to 
prevail at all times. The gases burning initially 
at the center of the vessel expand to a pressure, 
P, compressing thereby the unburned phase 
ahead. The work done on the unburned phase 
raises its temperature from its initial to some 
higher value. The burned phase loses an amount 
of thermal energy equivalent to the work done. 
As the flame progresses the pressure increases, 
further compressing the unburned phase as well 
as the burned phase; in consequence the temper- 
atures of both phases are increased. 

The principal object of this report is to 
determine the velocity of the flame relative to 
the moving unburned gas at any point in the 
enclosed sphere or, in other words, the velocity 
of propagation of the reaction zone into the 
unburned gas if the latter is considered to be 
at rest. This is the true speed of the flame. 
It is smaller than the flame velocity relative to 
the fixed coordinates of the bomb by an amount 
equal to the mass movement of the gases. 


DEVELOPMENT OF THE METHOD 


The first step toward solution of the problem 
is to develop the relation between the amount of 
gas burned up to a certain time and the pressure 
established at that time. From this equation and 
the experimental time-pressure record of the 
explosion the amount of gas burned during any 
time interval can be calculated. This amount of 
gas will be understood to occupy an annular 
shell. If the time interval is the unit of time, 
then the thickness of this shell just before 
burning represents the distance the reaction 
zone penetrates into the unburned gas in unit 
time, that is, the desired flame velocity. The 
thickness of the shell can be calculated only 
when the amount of gas contained in it, the 
prevailing pressure, and the radius of the shell 
and its temperature are known. The latter is 
the temperature of the unburned phase. 
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The relation between pressure and the amount 
of gas burned has been developed by Mache.®: ® 
In reproducing this equation we shall employ 
more familiar symbols relating to mols instead 
of unit masses, as used by this author. Let 


V =the total volume of the bomb. 

v =the volume of the burned gas at pressure P. 

vj=the volume that the mass of gas enclosed in 
occupied before ignition at initial pressure P; and initial 
temperature 7;. 

P,.=the final pressure (maximum) established when the 
entire contents of the vessel have reacted. 

C, and C, =the mean molar heat capacity of the burned 
and unburned gases, respectively, at constant pressure 
between temperatures 7); and 7, and 7; and 7,, respec- 
tively, where 7, is the temperature of the unburned phase 
and 7; is the temperature bordering directly on the burned 
side of the reaction zone. T; is identical with the tempera- 
ture the gases reach when allowed to burn at constant 
pressure at temperature 7,,. Likewise 7); is the temperature 
reached when the gases burn at constant pressure at 
temperature 7;. 

vo and y,=the ratios of heat capacities at constant 
pressure and constant volume of the gases in the burned 
and unburned phases, respectively, for the same tempera- 
ture ranges as C, and C,. 

nm, and m)»=the number of mols per unit mass of un- 
burned and burned gases, respectively. 


In general, indices 7 (initial) and e (end) refer 
to the states before ignition and after complete 
combustion, while indices « (unburned) and ) 
(burned) refer to the unburned and burned 
phases during the progress of the combustion 
from the center to the wall of the bomb. 

Furthermore Mache introduces the relation 


(np/Nu) CT, — Ce 2 =K, (1) 


where K is a constant. The only assumption is 
that mean molar heat capacities C, and C, are 
practically constant over a narrow range of 
temperature. For C, the maximum range of 
temperature extends from 7; to the temperature 
of the unburned phase at its maximum com- 
pression, namely, 7... Correspondingly, the 
maximum range of temperature for C, extends 
from 7}; to 7;,, the latter of which is the temper- 
ature reached when combustion (at the end) 
takes place at constant pressure at temperature 
T.... Neither C, nor C, are exactly constant in 


8 L. Flamm and H. Mache. Wien. Ber. 126, 9 (1917). 
*H. Mache, Die Physik der Verbrennungserscheinungen, 
1918. 
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this temperature span (100° to 200°), but for all 
practical purposes Eq. (1) appears to be fulfilled, 
as will be shown later. 

Then if v;/V is the fraction of the total gas 
burned at the time pressure P is established, 


' RT¥(P./P;:—P/P:) 
RT .(Y0— Yu)/(Y¥u-1) + (9-1) K 


where R is the gas constant. 

The fraction v;/V must equal 1 when P=P, 
and zero when P= P;. The first condition is met 
by Eq. (2). The second condition allows the 
determination of K by Eq. (2) using pressure 
and temperature data of an actual explosion 
experiment. One must demand that this value 
of K agree reasonably well with the value of K 
calculated by Eq. (1) using band spectra specific 
heats. Furthermore, the practical constancy of 
C, and C, must be demonstrated by the agree- 
ment in the values of K calculated by Eq. (1) 
for different combinations of 7, and T,. These 
comparisons will be made later in the application 
of the method to gaseous ozone explosions. 

T,, is given by the adiabatic law: 


a = Ti(P/P) oe ™, 





(2) 


(3) 
T, is given by the thermodynamic relation: 
T,= (dp /(no)Cpru™) +T., (4) 


where g, is the heat liberated per unit mass of 
gas at constant pressure and Cpr,,7° is the mean 
molar heat capacity of burned gas at constant 
pressure between 7), and 7;. 

It may be desirable to show briefly how Eq. 
(2) was obtained. 

Consider the mass of gas at pressure P 
enclosed in an elementary volume dv, in the 
form of an annular shell. Before the gas is burned 
it is compressed by the rising pressure, according 
to the adiabatic law, 


P' rudy, = constant = P;!!/ idy,, 


(2a) 


where dv; is the shell volume occupied by the 
same mass of gas before ignition occurred. It 
appears from Eq. (2a) that the change of state 
of the unburned phase is unambiguously given 
by the initial conditions for any pressure P. 
When the mass of gas in the element considered 
burns, it suddenly changes its volume and 
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temperature. Its subsequent changes of state 
with increasing pressure in the vessel are 
governed by the same adiabatic compression law, 


namely, (2b) 
2b 


However, the conditions existing when this 
subsequent compression starts are different for 
the various burned shells, considering that the 
innermost shells, for example, burn practically 
at pressure P; from an initial temperature close 
to T; while the shells near the wall burn practi- 
cally at pressure P, from an initial temperature 
close to Tx. 

In order to apply the adiabatic compression 
law to a shell after it has burned, its temperature, 
T,, immediately after burning must be known. 
For this purpose Mache introduces Eq. (1). 
Eqs. (1), (2a) and (2b) in conjunction with the 
condition of constancy of the total mass and 
volume, transform to Eq. (2) by a series of 
mathematical operations. 

It can readily be shown that at the moment 
maximum pressure is established the gas is 
hotter in the center than it is near the wall." 
An elementary mass of gas at the center burns 
by expanding at practically constant pressure P,. 
It is subsequently compressed to nearly its 
original volume by the combustion of the rest 
of the gas in the vessel. The latter work of com- 
pression exceeds the former work of expansion 
because the compression of the elementary mass 
takes place at a steadily increasing pressure 
from P; to P., while the expansion takes place 
at the lowest pressure P;. An elementary mass of 
gas near the wall, on the other hand, is first 
compressed under a pressure which rises from P; 
to P. and then expands (on combustion) to 
approximately its original volume at the maxi- 
mum pressure P,. The latter work of expansion 
is obviously larger than the former work of com- 
pression. Thus the elementary mass of gas near 
the wall loses some of its energy while the 
elementary mass of gas in the center gains 
energy in excess of the chemical energy released 
within it. 

The radius of any shell is calculated as follows. 
According to the gas law, the total number of 


P"'/7+dy, = constant. 


10 Experimental proof of a temperature gradient was 
first given by Hopkinson, Proc. Roy. Soc. A77, 387 
(1906). 
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mols in the bomb before ignition is given by 
Ni= VP;/RT;. (5) 


The number of mols left unburned at any 
instant is given by N,(1—7;/V). The volume 
occupied by this number of mols is (V—v»). 
By using the gas law once more, 


(V—v)P=N,(1—2;/V)RT,. (6) 
Substituting Eq. (5) and transposing, 


Up Fils vi 
MPT) g 
V PI; V 


Since v, is the volume of a sphere containing the 
expanded burned gas, its radius 7, is given by 


ef2E(2) 


where 7 is the radius of the bomb. 

The flame speed is derived from these relation- 
ships as follows: From Eq. (2) find 7;, the radius 
of the burned phase when no expansion whatso- 
ever occurs, that is, 7;/ro = (v;/ V)*. Plot r; against 
time, the latter being known from the time- 
pressure record. Draw three more graplis in 
order to correlate 7:,, P, and ratio r;/r,, plotting 
these variables against 7;. From the first plot 
the differential quotient dr;/dt is obtained as a 
function of r;. As indicated above, dr;/dt repre- 
sents the thickness of a shell about to be tra- 
versed by the flame front in unit time under the 
condition that no expansion of the adjacent 
burned phase occurred; that is, the temperature 
and pressure of the shell are 7; and P;. But by 
expansion the radius of the burned phase has 
increased from r; to 7, the pressure to P, and 
the temperature of the unburned phase to 7,. 
The volume of this shell of thickness, dr;/di 
before expansion of the burned phase, was 





(dr;/dt)4x(r;)*, (9) 
while after expansion the volume of the shell is 
Sy am 
dt T; P 
This volume also equals 
Aares, (11) 


where S is the actual thickness of the shell 
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which is just about to burn and whose com- 
bustion will be complete in unit time. S therefore 
represents the desired velocity of the flame: 


’ _ oF TP; 
MOVES 
dt 


T,,, P, and r;/ry are obtained from the last three 
plot ; for various values of dr;/dt. 

Tie velocity of the flame in any part of the 
ves: -l, the initial temperature 7, of the gas 
abcut to undergo combustion, the temperature 
T, immediately behind the flame front (from 
Eq. (4)), and the existing pressure P are now 
known. 


APPLICATION OF THE METHOD TO GASEOUS 
OZONE EXPLOSIONS 


We have applied the above system of equations 
to explosions of mixtures of gaseous ozone and 
oxygen. 

Ozone decomposes according to the reaction 


Os (gas) = 1202+AH, 


where A/J = 33,920 calories at 18°C and constant 
pressure. If a spark is passed in gaseous ozone 
or a mixture of ozone and inert gas, an explosion 
takes place whose intensity is dependent upon 
the degree of dilution with inert gas. In order to 
obtain flame speeds of the order of 10? to 10* cm 
per second the gaseous ozone was diluted with 
pure oxygen, the ratio oxygen /ozone being varied 
from about 1 to 4. 

Designating the ratio mols O2/mols O; by m 
and the energy liberated at constant pressure 
per mol of ozone by AH, Eq. (4) becomes 


Th= (AH/(1.5+m)Cpru?)+Tx (4a) 
and Eq. (1) becomes 
[(1.5+m)/(14+m)]C,7,—C.T.=K. (1a) 


We shall now present calculations of the 
constant K, using, on the one hand, Eq. (1a) 
with band spectra heat capacities and different 
combinations of 7, and 7, and, on the other 
hand, Eq. (2) with experimental pressure and 
temperature data from our explosions. For all 
practical purposes we may identify the mean 
molar heat capacities C, and C, with the true 
molar heat capacities of the burned and unburned 
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gases at the mean temperatures between 7%; 
and 7, and 7; and 7,,, respectively. Accurate 
values of the true heat capacity of oxygen -are 
known.® The heat capacity of ozone is known 
from band spectra with sufficient accuracy for 
this purpose. (These were kindly furnishe. us 
by Dr. L. S. Kassel.) , 
We shall employ the results of a typ:cal 
explosion, experiment 32, series 3 (3:32). _ , 


7; =302°K, P,.=4814 mm Hg;.- 
Tuc=494°K, P;=492 mm Hg, 
Tp = 1919°K, m=1.016, 

Toe = 2069°K, 


m 1 
Cu rie = ——C (02) rie ° +———-C(o;) Ti a ( 1 3) 
1+m 1+m 
The heat capacities are at constant pressure. 


C,4% = (1.016/2.016) X7.187 + (1/2.016) x 10.50 


"302 
= 8.83 cal./deg. mol, 
C269 = 9.207 cal./deg. mol., 


From Eq. (1a) 


K =(2.516/2.016) X9.207 X 2,069 —8.83 x 494 
= 19,410 calories, 


K=(2.516/2.016) X9.207 X 1,919 —8.83 x 302 
= 19,370 calories. 


From Eq. (2), 
RT,((P./P:i) —-1] 


{= 
RT(¥6— Yu)/(Yu- 1)+(yo- 1)K 
= 5266/(—33+0.275K), 





K = 19,280 calories. 


The good agreement of the K values shows 
that Eq. (1) is fulfilled for all combinations of 
T, and 7, and that Eq. (2) meets the require- 
ments of initial conditions. Therefore values for 
C, and C, based on our explosion data may be 
evaluated by a simple mathematical transfor- 
mation of the above equations. This procedure 
provides an independent method for determining 
the heat capacity of ozone." In all the following 
calculations these values of C, and Cy, are 
employed; all values are thus based on data 
obtained from the experiments under consider- 
ation. 


ee 


" Data on the heat capacity of ozone are presented in 
the second paper following, p. 294. 


TABLE I. Calculation of flame speed, S, in ozone explosions 
° . 
from time-pressure records. 








Time, P, P, S,cm 
sec. mm mm ri ’b Tims Th, per 
x10 Hg Hg Yb cm °K °K sec. 





Experiment 3:13 
m =3.054, P; =624 mm Hg, P,=3572 mm Hg, 7; =300.0°K 


0 624 624 0.599 1.02 300 1239 55.0 
20.70 665 630 601 2.99 301 1240 60.7 
27.63 712 650 603 4.96 303 1241 61.5 
34.52 795 700. ~—«615 7.04 309 1247 63.6 
37.98 862 785 .637 8.86 317 1253 65.5 
41.42 991 885 666 10.15 327 1259 74.3 
44.88 1141 1070 702) 11.33 343 1274 86.1 
48.35 1353 1324 745 12.41 361 1289 100.0 
51.79 1661 1622 791 13.20 380 1304 114.0 
55.25 2155 2024 845 13.83 401 1322 134.2 
56.97 2499 2560 907 14.38 427 1343 158.0 
58.72 2948 3220 970)» «=14.80 451 1362 178.0 
59.41 3192 
59.75 3324 
60.44 3572 





Yu=1.328, Cpruc™®*=8.18, Cpri™*=7.935, Ty.=1373, 
Toi = 1239, Tye =462, ri/ro(for P = P;) =0.599. 





Experiment 2:09 
m= 1.497, P;=367 mm Hg, P.=3029 mm Hg, 7;=301.1°K 


0 367 370 = 0.535 1.83 301 1680 141 

8.66 445 389 546 4.96 305 1684 165 
9.70 507 400 555 6.54 307 1685 187 
10.38 564 445 .568 8.20 314 1689 203 
11.08 626 575 610 10.54 333 1704 224 
11.77 698 690 638 11.48 347 1714 237 
12.98 910 820 612 1225 361 1726 257 
13.57 1045 980 708 12.90 376 1738 282 
14.34 1313 1185 750 13.46 392 1750 314 
14.89 1582 1420 791 13.87 408 1763 352 
15.31 1845 1695 833 14.22 425 1777 395 
15.62 2110 2000 878 14.48 441 1790 435 
16.17 2627 2825 978 14.92 477 1818 563 
16.30 2812 
16.52 3029 


yu=1.291, Cppuc™’*=8.445, Cpri7*'=8.21, T,.=1825, 
T»i = 1680, Te = 484, ri/ro(for P = P;) =0.531. 








Experiment 3:32 
m=1.016, P;=492 mm Hg, P,=4814 mm Hg, 7;=301.5°K 
492 495 0.502 0.92 302 1922 160 


0 

5.61 663 510 508 4.11 304 1923 219 
6.30 777 575 526 7.28 312 1928 287 
7.00 958 730 = .560 3610.05 317 1930 309 
7.52 1164 970 607 11.76 350 1958 379 
7.79 1314 1355 675 12.95 376 1976 495 
7.97 1469 2042 765 13.92 411 2003 615 
8.28 1692 2735 840 14.32 438 2021 693 
8.66 2199 3200 .880 14.55 453 2035 733 
8.94. 2700 3760 925 14.76 468 2044 747 
9.11 3134 4440 975 14.92 486 2063 780 
9.32 3702 

9.49 4252 

9.59 4542 

9.67 4814 





Yu=1.276, Cprue™®*=8.565, CpriT?*=8.34, T,.=2069, 
To; = 1919, Tue 494, r;/r»(for P=P;) =0.500. 
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The flame speeds for a considerable number of 
explosions were calculated, but only the results 
of three representative ones will be presented, 
namely, a moderately slow, a medium, and a 
moderately fast explosion. The results are tabu- 
lated in Table I. The radius of the bomb is 
150 mm. 

The first column for each experiment gives the 
time in milliseconds at which the absolute 
pressure in the second column was reached. 
These values were obtained from the experi- 
mental time-pressure record of the explosion. 
The next six columns show the results calculated 
from the first two columns by the method given 
above. The sixth column of the group contains 
S, the flame speed relative to the moving gases 
at various radii, 7, given in column 3 of this 
group. In column 5 of this group the corre- 
sponding 7,, calculated from Eq. (4a), is given 
for convenience. The Cpr,7° in this equation 
may be found readily by interpolating between 
the two limits given for Cpru.7** and Cpr;7", 
which are indicated under each experiment in 
Table I and which were obtained from earlier 
work.5 

r;/r, at the beginning cannot be determined 
easily from the above equations because both 
values approach zero. It is evaluated in the 
following way. Assume that a very small number 
of mols, dn, burns at the beginning. The volume 
it occupies, if no expansion occurs, equals 


vi=dn(RT;/P.), (14) 
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where 7; is the temperature calculated from the 
gas law and the maximum pressure P, (see 
below). After expansion the temperature de- 
creases to 7»; and the pressure to P;. Then 
the new volume equals 


v,=dn(RT>;/P;). - (15) 
Combining, 
v;/0»= (P;i/P.)(Tz/T») (16) 
and 
ri/to= ((Pi/P.)(Tz/T>) }. (17) 


TEMPERATURE GRADIENT IN THE BOMB AND ITs 
EFFECT ON THE DETERMINATION OF 
SPECIFIC HEATS BY THE Ex- 

PLOSION METHOD 


Before developing the above treatment we 
attempted to devise a simplified method for 
determining flame speeds by assuming that the 
temperature is equalized throughout the burned 
phase at any time. However, the values of 7; 
calculated and the values of 7, found from 
Eq. (4a) varied widely, indicating a rather steep 
temperature gradient from the center of the 
bomb to the wall which could not be neglected. 

Mache? developed a formula for the tempera- 
ture gradient in a spherical bomb at the instant 
combustion is complete. If Te(:;v) is the temper- 
ature at the surface of the sphere which encloses 
a fraction, v;/V, of the total gas burned, then 








vo—1 1+m } (yo-D/ yor XK Yu Py (uD! 14 
Te(es/v) = r(—) | + (—) | (18) 
- Vb 1.5+m P RT; Yu-—1 P; 


in which all the symbols are known. 

In the explosion method for determining the 
heat capacities of gases the maximum pressure 
developed is used to calculate the temperature 
of the gases, 7;, from the gas law under the 
assumption that the temperature is equalized 
throughout the bomb. As already pointed out by 
Hopkinson” this procedure would be strictly 
justified if the heat capacity of the gases did 
not vary with temperature. In the normal case, 
however, where one deals with polyatomic 
molecules and their dissociations the maximum 





pressure recorded is a little less than the maxi- 
mum pressure which would have been recorded 
had the temperature throughout the bomb been 
equalized without heat loss. We have made 
extended calculations of this effect for the 
present type of explosion including the effect of 
dissociation and find that the difference between 
these pressures is from less than 0.1 percent to 
less than 0.3 percent in the extreme case of the 
highest temperature explosions. This work wil 
be published soon. 

For other types of explosive mixtures one may 
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DETERMINATION OF THE SPEED OF FLAMES 


estimate the effect to be very small except in 
explosions where large dissociation occurs. But 
even in such cases the pressure difference can be 
calculated and therefore this effect does not 
constitute an objection to the explosion method. 

Fig. 1 shows graphically the conditions in the 
bomb while combustion is in progress and before 
and after combustion for explosion 2:09. The 
hemisphere represents one-half of the bomb. The 
abscissa is the radius of the bomb, 15.00 cm. 
The upper quadrant shows the positions reached 
by the flame front after equal time intervals. 
The time interval in this instance is 1.65 x 107% 
sec., which is exactly 1/10 of the total explosion 
time. The pressure existing in the bomb for 
these positions of the flame front are indicated 
on the left. Temperatures 7°, and 7, in front of 
and behind the flame front, respectively, are 
given to the right of these. Thus the upper 
quadrant provides a complete picture of condi- 
tions during the progress of the flame. It should 
be emphasized that temperature 7°, in front of 
the flame is equalized throughout the entire 
unburned phase, whereas temperature 7), behind 
tne flame is localized in the immediate vicinity 
of the flame front. From the rear of-each flame- 
front position the temperature rises toward the 
center of the bomb. The arrows represent flame 
speeds S relative to the moving gas when the 
flame reaches the positions indicated. 

The solid curves in the lower quadrant desig- 
nated by numerals with the subscript 7, which 
denotes pre-ignition conditions, mark the bound- 
aries of the spheres occupied by the same mass 
of gas contained in the corresponding spheres in 
the upper quadrant whose boundaries are marked 
by numerals with the subscript 6. The difference 
between the corresponding i and 0b spheres 
results from the expansion which occurs on 
combustion. The two-shaded portions represent 
this for the mass of gas burned during the first 
1/10 of the total explosion time. 

The dotted curves designated by numerals 
with the subscript e mark the boundaries of the 
spheres occupied by the same masses of gases 
Corresponding to the z and 6 states after the 
explosion is complete or when the maximum 
Pressure is reached. They show the extent to 
which the b spheres have been compressed by 
the subsequent combustion. The radii in centi- 





Fic. 1. Experiment 2:09. Showing the hemisphere of the 
bomb. Abscissa =radius of bomb in cm. Total explosion 
time 16.5X10-* sec. Upper quadrant: Contains speed of 
flame in various parts of bomb, Curves marked 13, 23, etc., 
are volumes occupied by gas burned in equal time intervals, 
in.this case 1/10 of total explosion time. Pressure in bomb, 
temperature 7, in unburned phase and temperature 7) 
immediately behind flame front when flame reaches points 
indicated after successive tenths of total explosion time, are 
given to the left. Lower quadrant: Solid curves, 11, etc., are 
volumes occupied by same mass of gas enclosed in 1, etc., 
for condition that no expansion of the gases occurred. 
Dotted curves le, etc., are final volumes occupied by same 
mass of gas in 11, etc., or 1b, etc., after 1b, etc., has been 
compressed due to subsequent burning of the outer shells. 
Temperature gradient existing in bomb when combustion 
is complete, is given to the left. 


meters of the dotted spheres were calculated from 
the formula given by Mache, which by a simple 
transformation may conveniently be written 


re=rel(P/P.)!/ 7}. (19) 


The temperatures existing at the boundaries of 
these spheres after combustion is complete may 
be calculated from Eq. (18) for the same values 
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of P used in Eq. (19). These temperatures, 
given to the left of the lower quadrant, may be 
compared with temperature 7;=2072°K, which 
would exist if the temperature were equalized. 

The temperatures of the inner portions of gas 
are somewhat too high because dissociation of 
oxygen has been neglected in these calculations, 
and therefore the value of C, used is somewhat 
too small. Radii r., however, would be affected 
only slightly. We believe the picture, as pre- 
sented, is nevertheless substantially correct. 

Table II contains the numerical data shown 
in Fig. 1. 

Fig. 2 is a photograph of the time-pressure 
record of explosion 2:09 from which Fig. 1 was 
derived. 


TABLE II. Numerical data for experiment 2:09 used in Fig. 1. 











Time, ; S,cm 
sec. ri, 4) Y, Tu. Tb, P, per Te, Ye, 
X<103 cm rb cm °K °K mmHg sec. °K cm 
0 0 0.531 0 301 1680 367 (140) 2670 0 
1.65 .79 534 1.48 302 1680 369 146 2667 -86 
3.30 1.62 538 3.01 303 1681 374 151 2657 1.75 
4.95 2.52 544 4.63 305 1683 385 163 2644 2.71 
6.60 3.52 554 6.36 307 1687 400 182 2633 3.76 
8.25 4.65 569 8.17 314 1692 447 202 2571 4.97 
9.90 5.91 595 9.94 329 1701 530 218 2496 6.32 
11.55 7.36 -640 11.50 347 1716 685 237 2380 7.81 
13.20 9.10 -708 12.85 376 1738 980 281 2225 9.58 
14.85 11.45 815 14.05 416 1776 1560 370 2013 11.76 
16.50 15.00 1.000 15.00 484 1825 3029 (650) 1825 15.00 
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PRESSURE, ABSOLUTE 





ee ae Tal ee: eee 
9.9 13.2 16.5 19.8 22,1 25.4x10-% sec, 


TIME—> 





| 
0 33 6.6 


Fic. 2. Time-pressure record, explosion 2:09. 


It is of interest that in the first part of the 
record the curve practically coincides with the 
zero line. It has been stated in the literature 
that this circumstance is due to the existence of a 
time lag between the passage of the spark and 
ignition. However, this conclusion is not ad- 
missible.”2 In the present case, for example, ref- 
erence to Figs. 1 and 2 shows that when 1/5 of 
the total explosion time has elapsed the pressure 
rise becomes just readable at 7 mm Hg. In this 
same time the flame front has progressed 3 
centimeters or 1/5 of the radius of the bomb. 


22 Payman (Fuel 1, 36 (1922)) also realized this. 
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A method is described for determining the heat capacity 
of explosive gases at temperatures up to several hundred 
degrees above room temperature, by exploding the gas 
whose heat capacity is to be determined in a spherical 
bomb. The same data are required as are needed to 
determine the heat capacity of the burned gas by the 
explosion method, namely, the maximum pressure devel- 
oped in the explosion, the initial temperature and initial 
pressure, the composition of the gas mixture, and the heat 


of reaction. Equations have been developed, and the 
method is described in detail. The method is particularly 
useful for the determination of heat capacities of easily 
decomposed gases, such as O; and ClO, for which the 
flow method is unsuitable. An interesting consequence of 
the method is that it permits simultaneous determination 
of the heat capacity of carbon dioxide at high and low 
temperatures by exploding mixtures of gaseous ozone and 
carbon dioxide. 





ITHERTO, experiments on explosions of 

gases at constant volume in a spherical 
bomb have been utilized to determine the mean 
heat capacity of the products of the explosion 
between the initial temperature of the gases and 
the final mean temperature calculated from the 
maximum pressure by means of the gas law. 
The present paper deals with a method for 
employing the same explosion experiments to 
determine the mean heat capacity of the 
reactants between the initial temperature and 
the highest temperature reached by the unburned 
gases due to their adiabatic compression by the 
expanding burned gases during the explosion 
process. 

The heat capacities of gases between room 
temperature and several hundred degrees above 
usually are determined by the flow method; 
however, there are gases, such as O; and C1,O, 
for which this method is not suitable due to the 
ease with which they decompose. Therefore, the 
method we shall describe is especially useful 
in such cases. One particularly interesting con- 
sequence of the method is that explosions of 





‘Published by permission of the Director, U. S. Bureau 
of Mines. (Not subject to copyright.) 

* Physical chemist, Pittsburgh Experiment Station, 
U.S. Bureau of Mines, Pittsburgh, Pa. 


mixtures of ozone and carbon dioxide allow 
simultaneous determination of the heat capacity 
of carbon dioxide at high and at low tempera- 
tures. 

The method is based on the principle that the 
work of compression done on the unburned 
phase and its temperature rise due to compression 
may be determined independently, yielding 
directly the heat capacity of this phase. If the 
original gas is composed of a single component, 
the heat capacity so obtained is that of this gas. 
If it is composed of a mixture of two gases the 
heat capacity of one should be known to deter- 
mine the heat capacity of the other. In practice 
it is usually necessary to dilute the explosive 
mixture. 

All equations needed are contained in the 
preceding paper.‘ The equation relating the 
fraction of the total gas burned with the pressure, 
P, in the bomb, is: 


; vi RT (P./P;—P/P;) 
——— . (1 
V RT u(yo— ¥u)/ (Yu - 1) + (0 - 1K 





3’ Member of staff, Coal Research Laboratory, Carnegie 
Institute of Technology, Pittsburgh, Pa. 
4 See preceding paper, p. 283. 
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where v; is the volume of the initial gas mixture burned 
and V the total volume of the bomb. 

v;/V is the fraction of the total gas which has burned 
when the pressure, P, is established. 

P. is the maximum pressure developed in the explosion. 

P; and 7; are the initial pressure and temperature of 
the gas mixture. 

T,, is the temperature in the unburned phase when the 
fraction v;/V has burned. Thus 7,,>7;. 

Yu is the ratio of mean heat capacities of the unburned 
phase at constant pressure and constant volume between 
the temperatures 7; and 7,., where 7,. is the highest 
temperature reached by the unburned phase just prior to 
complete combustion or near the wall of the bomb. 

yo is the same ratio for the burned phase between the 
temperatures 7); and 7;., where 7); is the temperature 
reached when the gas burns at constant pressure from the 
initial temperature 7;, and 7», is the temperature reached 
when the gas burns at constant pressure from the initial 
temperature Tx. 


It may be remembered from the previous 
paper that the temperature 7), or 7;, is reached 
when the first or last infinitesimally small 
amount of gas is burned. Since the combustion 
of this amount of gas increases the total pres- 
sure only infinitesimally, the process takes place 
at constant pressure. 

K is defined by Eq. (2) and is a constant. 


K= (ns/nu) CT. — CuTu, (2) 


where m, and m, are the number of mols per unit mass of 
burned and unburned gases, respectively. 
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Cy, and C, are the mean molar heat capacities per mol 
at constant pressure of the burned and unburned gases 
between 7); and 7}-, and 7; and Ty., respectively. 

Since the energy C,(7..-—7;) must equal the 
energy (/nu)Co(Toe— 7»), it follows that, 


K= (n»/Nu) CToe— Cul ue 


=(nv/nu)CoTvi—CuTi. (2a) 


The temperature 7;; is given by the thermo- 
dynamic relation: 


Toi=Qp/neCpri™’i+Ti, (3) 


where gq, is the chemical energy liberated at constant 
pressure per unit mass of gas, and Cp7;7*' is the mean 
molar heat capacity of burned gas at constant pressure 
between 7; and 7},. 


Likewise, the temperature 7,, is given by 


Toe=Qp/n eCPrue’*+ a (4) 


The temperature 7,,, is given by the adiabatic 
law, 


Tuc™ T (P./P Vu, (5) 


On imposing the condition P=P; and there- 
fore v;/V=0 in Eq. (1), dividing by RT;, and 
substituting the right side of Eq. (2a) for K, 
one obtains, 





~ C Ti 
1-(5-1) /[*= - ~(y- DEH “| (6 











Since 
¥ C %  C. 
—=— and =—, (7) 
yo—1 R Yun 1 R 
then 
Yo Yu u 
—(y—-1)—=—-Y, (7a) 
Yu-1 R 
and 
(yo—1)C,/R= vo. (7b) 


Eq. (6) then becomes 





R Nu T; 





P./P;-1 
(»/nu)(Ti/T3) - 1 





Eq. (8) provides a method for determining 
the mean heat capacity of the burned gases 
between 7}; and 7). which, over this narrow 
range of temperature, is practically identical 
with the true specific heat at the mean between 
these temperatures. 
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Now repeat the operations performed in 
obtaining Eq. (6) but substitute for K the left 
side of Eq. (2a), transpose, and divide by y,—1. 
Then 


Yu— Yb 


(yu—1)(ve—1) 


Np P, 1 
-~csts,—RT(—"-1) . (9) 
P; 


Nu vo—1 


RT;+ on 5 as 








With the relations (7), 
(vYu- 0) /(Y¥u-1)(yo—1) = (Cu— Cs) /R, (10) 
so that Eq. (9) becomes, 


C.A(Tue— Ti) +OQ.7T; 
=(n/mu)CoToe—RT(P-/Pi—1)/(ye—1). (11) 


Eq. (8) may be written 


Yb c ny Ti 
—_—_ =—=(—- —————p, (12 
are VEE 


and also 


1 m Toi e my Toi 
i -1)/(- th ~). (13) 
— n, T ; 5 Ny i 


1 P; T; 














Substituting (12) and (13) in (11) and trans- 
posing, 


R(P.—P3i) Toe— To; 
~ (n,/m)Po —(Toi/Ti) Pi tu 





(14) 





C,, is the mean molar heat capacity of the un- 
burned gas at constant pressure between 7; 
and Ty.. 

If, for example, the original explosive mixture 
is composed of two components, / and m, having 
mol fractions x and 1—x, respectively, 


Cy=xCi+(1—-«)Cn, (15) 


where C,; and C,, are the mean molar heat 
capacities at constant pressure between 7; and 
T..2 of the components / and m, respectively. 
Practically, over this narrow range of tempera- 
ture, these mean heat capacities are very nearly 
the true heat capacities at the mean temperature 
(Ti+ Teall a 

The temperatures 7;;, 7>., and 7, are deter- 
mined from Eqs. (3), (4) and (5) in combination 
with Eq. (14) by the method of trial and error. 
It is necessary to know, therefore, the tempera- 
ture dependency of the heat capacity of the 
burned gas, which is obtained by the explosion 
method in the usual way and for which the same 
experiments may be utilized. In practice, there- 
fore, this information is available. In working 
out specific examples it will be found convenient 
to employ a plot of the mean heat capacity of 
the burned gas between some low standard 
temperature such as absolute zero, and 7;. 
Chrue?”, for example, may be found from the 
plot by employing the relation 


Cphrue™”* 


= (TeeCpro?®* — TueCpro’“*)/(T be | y am Fe (16) 
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The mean heat capacity of ozone between 300 and 476°K 
has been determined from explosions of mixtures of ozone 
and oxygen in a spherical bomb by the method the authors 
described in a previous paper. The value obtained is 
10.94 calories per mol of ozone. The value obtained from 
band spectra is 10.39 calories per mol. The difference 
cannot be accounted for by incorrect assignment of the 
optical frequencies of ozone, and it is larger than the 
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probabie experimental error. It is proposed that the 
difference may be due to a contribution to the heat 
capacity ofan electronic level lying close tothe ground state 
of the ozone molecule. It is estimated that this supposed 
level has an energy of 0.1 to 0.2 volt. The existence of 
low-energy electronic levels probably implies an isosceles 
triangle structure and paramagnetic properties of the ozone 
molecule in the normal and the low excited states. 








[N the preceding communication‘ a new general 

method for determining the heat capacity of 
explosive gases was described. The present paper 
will discuss application of the method to the 
determination of the heat capacity of gaseous 
ozone. 

Mixtures of gaseous ozone and oxygen in 
varying proportions were ignited at the center 
of a spherical bomb and the maximum pressures 
P,, developed were accurately measured on the 
time-pressure records of the explosions, as de- 
scribed previously.’ Additional data needed are 
initial pressure P;, initial temperature 7;, the 
composition of the original mixture m, which is 
the ratio of the concentrations of oxygen and 
ozone, the heat of decomposition of 1 mol of 
ozone at constant pressure and room temperature 
equal to +33,920 calories, and the mean heat 
capacity of oxygen at constant pressure from some 
convenient reference temperature (say 0°K) to 
the temperatures of the explosions. Other sym- 
bols used are: 


Tue=the temperature of the adiabatically compressed 
unburned gas just prior to complete combustion or near 
the wall of the bomb; 


1 Published by permission of the Director, U. S. Bureau 
of Mines. (Not subject to copyright.) 

? Physical chemist, U. S. Bureau of Mines, Pittsburgh 
Experiment Station, Pittsburgh, Pa. 

3 Member of staff, Coal Research Laboratorv, Carnegie 
Institute of Technology, Pittsburgh, Pa. 

4 See preceding paper, p. 291. 

5 Lewis and von Elbe, J. A. C. S. 55, 504, 511 (1933). 
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Ty; and T,.=the temperatures established by burning 
a mass of gas at constant pressure from initial temperatures 
T; and T,,., respectively ; 

C,=mean molar heat capacity of unburned gas at 
constant pressure between the temperatures 7; and 7,,,; 

Cprue™** and Cpri7**=the mean molar heat capacities 
of oxygen at constant pressure between the temperatures 
indicated. 


The reaction which takes place is: 
O; (gas) = 1.502+33,920 calories. 


All the equations of the preceding paper remain 
the same, except for modifications indicated in 
Eqs. (1), (2) and (3) below. 


n,/n,=(1.5+m)/(1+m), (1) 
Ti = (33,920/(1.5-+m) Chri?) +T;, (2) 
Ty. = (33,920/(1.5+m)Cpruce? )+T ue. (3) 


Table I summarizes the data of a group of 
experiments and the results of the calculations. 

The last column shows the calculated mean 
molar heat capacity of ozone. The average of 
the results between 300 and 476°K is 10.94 
calories per mol. The average deviation from 
the mean is 0.30 calorie per mol, with the value 
given as the most probable value. 

The above result may be compared with the 
heat capacity of ozone calculated from existing 
optical data. Kassel® has calculated the free 
energy of ozone, based on Gerhard’s? frequency 


6 L.S. Kassel, J. Chem. Phys. 1, 414 (1933). 
7 Gerhard, Phys. Rev. 42, 622 (1932). 
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TABLE I. Determination of heat capacity of gaseous ozone, Cp(QOs). 








Experiment 


No. P P; T; Toi Tre 


Tue CP TueT CoriT* Cy Cp(02)riT Cp(0s) riT¥ 





367 1825 
414 
436 
474 
624 
513 
561 
639 
368 
366 


301 
301 
302 
300 
300 
301 


2.09 
2.41 
3.01 
3.12 
3.13 
3.14 
3.15 
3.16 
3.22 
3.23 
3.32 


1294 
1618 
1373 


1402 


301 
302 


1646.5 


1490.5 


11.2 
10.7 
11.3 
10.2 
10.7 
10.4 
12.1 
10.9 
11.2 
11.5 
11.2 


8.81 
8.35 
8.11 
8.16 
8.03 
8.07 
8.18 
8.19 
8.70 
8.80 
9.18 


484 
481.5 
464 
483 
462 
473.5 
460 
467.5 
481 
479.5 
494 


8.45 
8.34 


8.21 


~J ss) ss) SS) +) I 
push pum buh ums Gem emt Gumi eb feu fem pen 
OOK” D Ms 31-1100 








assignment 528, 1033 and 1355 cm™, and has 
kindly furnished us with values of the heat 
capacity of ozone determined from them. Be- 
tween 300 and 476°K the mean heat capacity 
is 10.39 calories per mol. Any reasonable change 
in these frequencies would hardly affect the heat 
capacity determined from band spectra by more 
than 0.05 calorie. The difference between the 
experimental and optical values is too great to 
be accounted for by experimental error. We 
suggest that this is due to the contribution of 
an electronic level lying close to the ground 
state of the ozone molecule. 

If a reasonable assumption is made regarding 
the statistical weight of this level with respect 


to the ground level, its energy can be estimated® 
to be 0.1 to 0.2 volt. 

The existence of one or more low energy 
electronic levels probably would imply that the 
ozone molecule has an isosceles triangle structure 
in both the normal and the low excited states 
and that the ozone molecule may be para- 
magnetic.’ There are no measurements of the 
latter, but it seems to be entirely feasible to 
determine the magnetic susceptibility of gaseous 
ozone. Concerning the structure of the molecule, 
the most recent spectroscopic investigations’ 


indicate that it is isosceles in form. 


8 Lewis and von Elbe, Phys. Rev. 41, 678 (1932). 
9 Private correspondence with Professor R. S. Mulliken. 
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This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing Editor 


Sufficient data are available to calculate the zero-point 
energy of some simple molecular crystals. These values 
may be compared with the corresponding heats of sublime- 


TABLE I. 








E(0°K) / 


M.W. @(max) E(O°K) AH(O°K) M.P. AH(0°K) 





He? 4. 9 20 24(4°K) 0.9 84% 
Ne! 20.2 63 142 425 24.4 33. 
A! 39.9 85 191 1830 83.9 10. 
H23 2. 91 205 216(20°K) 14. 95. 
No‘ 28. 68 153 1636 63.3 9. 
cos 28. 80 179 1845 68.1 10. 
NO*® — 30. 119 268 3680 109.5 7. 
O2? 32. 91 204 2083 54,7 10. 
HCI& 36.46 128 290 4845 161. 6. 
HBr® 80.92 92 207 5428 184.6 4. 
HI" = 127.9 73 164 6250 222.3 3. 








1 Landolt-Bornstein Physical Tables. 

2 Keesom and Keesom, Proc. Acad. Sci. Amst. 35, 736 
(1932). 

3 Fowler and Stern, Rev. Mod. Phys. 4, 697 (1932). 

* Giauque and Clayton, J.A.C.S. 55, 4882 (1933). 

5 Giauque and Clayton, J.A.C.S. 54, 2619 (1932). 

6 Giauque and Johnston, J.A.C.S. 51, 3208 (1929). 

7 Giauque and Johnston, J.A.C.S. 51, 2300 (1929). 

8 Giauque and Wiebe, J.A.C.S. 50, 116 (1928). 

® Giauque and Wiebe, J.A.C.§. 50, 2200 (1928). 

10 Giauque and Wiebe, J.A.C.S. 51, 1447 (1929). 


When light from a capillary discharge tube filled with 
hydrogen at a pressure of 1-2 mm of Hg passes for some 
time through a fluorite window into air at atmospheric 
pressure there is formed on the outside of the window a 
light film showing metallic luster. The film can be removed 
by swabbing the window with absorbent cotton moistened 
with alcohol. 
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tion calculated to 0°K. The zero-point energy may be 
calculated from the relation: 


E(0°K) = J ” (2%) F4°( hv /2)dv = 2.250 max) calories. 


The heats of sublimation were calculated to 0°K from the 
best recent available data. The results are given in Table I. 

In order for the crystal state to exist the sum of the 
zero-point energy, the attraction potential and the repul- 
sion potential must be negative. From the table it is seen 
that the ratio of the zero-point energy to the heat of 
sublimation is greater for the crystal of lighter mass 
particles. Hence a crystal of the lighter mass particles is 
the less stable one, in virtue of its high zero-point energy. 
Further it is reasonable to say that the more negative the 
above-mentioned sum, the more stable is the lattice to a 
change of state. As evidence for this the melting points of 
the monatomic solids in the table increase with the heat 
of sublimation. Likewise do the melting points of the 
diatomic solids. Of the latter the value for nitric oxide 
suggests the existence of (N2O,) units. The hydrogen halide 
crystals have an additional factor entering into the 
attraction potential, namely the electrostatic potential 
between the permanent dipoles. 

V. Deitz* 
University of Illinois, 
Urbana, Illinois, 
March 10, 1934. 


* National Research Fellow 1933-34. 






A more prolonged use of the tube gives rise to a deposit 
of white powder on top of the film which can easily be 
brushed off with dry cotton. Both deposits tend to prevent 
the passage of wave-lengths 1250-1500 Angstroms through 
the window. An imperfection in the fluorite casts a shadow 
in such a manner that no deposit is formed at a point 
where a straight line from the end of the capillary through 
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the imperfection cuts the surface of the fluorite. 

The powder has been found to be insoluble in distilled 
water but soluble in dilute hydrochloric acid. It forms 
after a shorter period of time if the air is moist than if it 
has been carefully dried. It is possible that light of ex- 
tremely short wave-length activates the oxygen of the air 
so that either (1) calcium oxy-fluoride (perhaps) is formed 
directly, or (2) hydrogen peroxide is produced first and 
then reacts with CaF, to form the deposit. In the latter 
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case it might be expected that a very strong solution of 
HO: applied to the window would produce some effect 
upon it, but upon making this test not the slightest effect 
could be detected. As yet nothing further is known about 
the nature of the deposit or of the reaction which causes it. 
FREDERIC PALMER, JR. 
Haverford College, 
Haverford, Pennsylvania, 
March 20, 1934. 


Note on the sp’ Configuration of Carbon, and Correction to Part III on CH, 


In part III'-? of my recent series on CHy,, I over- 
looked Edlén’s* new measurements and classifications for 
the spectra of C I, N II and O III, which appeared just 
as my article was being completed. The previously un- 
known sp?*S state has been found by Edlén in N II and 
O III, and may now be located in carbon I at approxi- 
mately 4.1 volts above the ground state s*p? *P by extrapo- 
lation from N II and O III with the aid of the irregular 
doublet law. A more accurate, in fact quite remarkable, 
method of extrapolation based on theoretical considerations 
has subsequently been developed by Bacherand Goudsmit.* 
It places 5S at 4.3 volts in carbon. Thus 5S is not quite as 
high above the ground state as my previous estimate 5-8 
volts based on the Slater-Condon-Johnson F-G relations, 
but is much closer to this estimate than to the value 1.6 
volts so often quoted in the literature. It is somewhat 
accidental that my estimate was as close as it is, for in 
applying these relations, I used the experimental values of 
the sp? 'D and sp? 'P terms, viz., 9.65 and 9.66 volts, usually 
listed in spectroscopic tables, such as for example Bacher 
and Goudsmit’s book. 

Edlén finds that the terms at these positions are certainly 
due to some other configuration, probably s*pd. The true 
sP’'D and sp* 'P levels have not yet been detected spectro- 
scopically, but may be located by Bacher and Goudsmit’s 
new method at 12.9 and 14 volts, respectively. (Because 
*S of N II was obviously out of line with 'P and 'D of CI 
in the old classification, I suggested in part II] that #S 
was due to some other configuration than sp’; it is now 
evident that the old 'D, 'P were the real culprits.) With 
the new classification and for more extensive data, the 
Johnson formulas fit experiment very much better than 
previously, in fact unusually well for F-G relations. The 
revised and extended version of Table I of part III is 
given in Table I. 


TABLE I. Energies, in volts, relative to the ground state s*p? *P. 








Observed 
I N Oll 


Calculated 
I NII Olll 


7.4) 

24.4) 
15.2 
23.7 

17.6) 
26.1 





5.8 
19.2 
11.8 
18.5 
13.5 
20.2 








The ‘“‘observed’”’ values in square brackets, unlike the 
others, are not really spectroscopic, but have been com- 
puted by Bacher and Goudsmit with their new extrapola- 
tion method. Their results are reliable to a few tenths of a 
volt, and better than the estimates 12.1 and 14.8 for 
1D, 1P which Edlén obtained with the irregular doublet 
law. The calculated values have been obtained by the 
Johnson formulas given in part III, with the parameters 
chosen so as to yield the proper positions for 5S, *S, *P. 
The corresponding magnitudes of the parameters are 
those labelled ‘‘empirical’’ in Table II; they would not be 
much altered if three other states were fitted instead of 
5S, 3§, 3P. We make this selection for the three because it 
makes the maximum deviation between the calculated and 
observed values in Table I less than with most other 
choices and so represents a closer approximation to the 
method of least squares. The Bacher-Goudsmit procedure 
is sufficiently accurate that it is not important that one 
of the three states fitted has been located in C I by this 
procedure rather than spectroscopically. 

The values in parentheses are those obtained in part III 
for carbon by fitting 'P, 3D, *P with the old classification. 
We have also listed for comparison the results of Beards- 
ley’s theoretical calculations by means of approximate 
wave functions obtained by a variational method. It is 
seen that the new classification largely removes the 
previous qualitative disagreement between the empirical 
values of the Johnson parameters and those computed by 
Beardsley. 

With the empirical values of the parameters listed in 
Table II, one calculates from Eqs. (7) and (8) of part III 


1J. H. Van Vleck, J. Chem. Phys. 2, 20 (1934). 

? The following minor errata may be noted in part II 
(ibid. 1, 219, 1933): add $(1—cot? 6) sin? 6N,,+cot? ON, 
to the expression (41); in the line before (41) read N2=0 
for N=0. 

3 B. Edlén, Zeits. f. Physik 84, 746 (1933); I am indebted 
to Professor Bowen for calling my attention to this 
reference and to the use of the irregular doublet law in 
connection with it. 

‘Bacher and Goudsmit have as yet published only a 
preliminary abstract on their method; Phys. Rev. 45, 
559 (1934). They have, however, very kindly informed me 
of the results of the application to carbon in advance of 
their detailed paper. 
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TABLE II. 
—Enppirical——_——_—_ Beardsley 
OF NIl O Ill Cl 
G; 2.18 (0.9) Sino 4.25 3.44 volts 
Fy 0.188(0.2) 0.290 0.400 0.29 
Wo 13.7 (41.1) 20.2 26.1 14.6 








that the energy of the ‘‘valence state” of carbon explained 
in section 2 of part III should range from 7.1 volts for the 
tetrahedral model n=} to 8.5 volts for the model p=1 
without hybridization. The effect of the L, S structure is 
thus not quite as trivial as we previously estimated, and 
favors the tetrahedral modal model slightly. Because of 


THE EDITOR 


the successful application of Johnson's formulas to the new 
data, our estimate of 7 volts for the anergy required to 
promote carbon to the tetravalent condition appropriate 
to CH, is now based on a much more certain calculation 
than formerly and is probably accurate to within half a 
volt. This estimate is fortuitously the same as that which 
we obtained in part III with the old values of the Johnson 
parameters. Hence sections 3 to 6 inclusive of part III, 
and so all the chemical applications in III, are mercifully 
unaltered. 
J. H. Van VLEcK 
Department of Physics, 
University of Wisconsin, 
Madison, Wisconsin, 
March 21, 1934. 


A Slight Difference in the Isotopic Composition of Oxygen Made by the Fractionation of Liquid Air and Oxygen of 
Ordinary Air! 


In an investigation of the fractionation of the isotopes 
of hydrogen and oxygen by the electrolysis of water, a 
series of measurements has been made on the difference in 
density of ordinary water and of the water formed by: 
(1) combining the electrolytic hydrogen with commercial 
oxygen; (2) combining the electrolytic oxygen with com- 
mercial hydrogen; (3) recombining the electrolytic gases 
as they come from the cell. The commercial oxygen, made 
from liquid air, was assumed at the start to have the 
isotopic composition of ordinary oxygen as it exists in air 
and water. The commercial hydrogen, made electrolytically 
in cells which have been in operation for years, was shown 
to be normal by combining it with ordinary atmospheric 
oxygen and comparing the density of the water thus 
formed with that of ordinary water. The difference in 
density was only 0.3 part per million. 

However, in every case the sum of (1) and (2) failed to 
check with (3) by amounts ranging from 1 to 3 ppm and 
always of the same sign. The discrepancy was finally 
traced to the assumption that the commercial oxygen was 
the same as the oxygen in normal air and water. This 
assumption was tested by preparing two samples of water, 
one composed of the dried commercial oxygen combined 
with dried hydrogen from a given cylinder, and the other 
composed of ordinary atmospheric oxygen and dried 
hydrogen from the same cylinder. After purification by 
distillations from alkaline permanganate and from phos- 
phoric acid followed by simple distillation, with the 
rejection in each case of small first and last portions, 
the water from commercial oxygen was found to be 


heavier by 1.4 and 2.4 ppm. After an additional distillation 
the difference was 2.8 ppm, giving an average of 2.2+0.5 
ppm. Assuming adequate purification, this difference can 
result only from a difference in the two samples of oxygen. 

The commercial oxygen was manufactured by the 
fractionation of liquid air in a column which produces 
nitrogen of 95.5 percent purity at the top, and apparently a 
slight fractionation of the oxygen isotopes occurs, so that 
the commercial oxygen has a density greater by about 2 
ppm and therefore contains approximately 2 in 10° more 
atoms of O'8, than ordinary oxygen. This difference may 
vary somewhat with the efficiency of different fractionating 
columns and may be of significance in other investi- 
gations. 

Before these measurements were completed Klar and 
Krauss’ found differences in the isotopic composition of 
different fractions of technical oxygen. However, the 
results reported in this letter may be of some value since 
they confirm the results of Klar and Krauss and, in addi- 
tion, show the magnitude of the change to be expected in 
the oxygen from an ordinary commercial cylinder. 

Epcar R. SMITH 

Bureau of Standards, 

Washington, D. C., 
March 27, 1934. 


1 Publication approved by the Director of the Bureau of 
Standards, U. S. Department of Commerce. 

2R. Klar and A. Krauss, Naturwissenschaften 22, 119 
(1934). 


The Magnetic Susceptibility of Ammonium Hexabromo-Hypoantimonate 


At the suggestion of Professor Linus Pauling, I have de- 
termined the magnetic susceptibility of ammonium hexa- 
bromo-hypoantimonate, (NH,)2SbBre, in order to obtain 
some information about the valence of antimony in this 
compound. 

The substance was prepared by the method of Ephraim 
and Weinberg.! It was dried by centrifuging. The product 





was a deep black powder, apparently similar to Ephraim’s 
compound. 

The determination of the specific magnetic susceptibility 
by the Gouy method led to a value, x = —0.036 X 10~ cg, 
the substance thus being shown to be diamagnetic. 


1 Ephraim and Weinberg, Berichte 42, 4450 (1909). 
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Ammonium hexabromo-hypoantimonate, therefore, does 
not contain unpolymerized complexes involving a quadriva- 
lent antimony atom, for such complexes would con- 
tain an odd number of electrons and hence give rise to 
strong paramagnetism. It is probable that the sub- 
stance is a double salt containing one trivalent and 
one quinquivalent antimony atom corresponding to the 
formula (NH,4),Sb™'Sb’ Bris. The observed diamagnetism 
could be made compatible with the presence of quadriva- 
lent antimony by assuming that the quadrivalent antimony 
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atoms occur in groups of two joined by shared electron- 
pair bonds, Sb : Sb, similar to the group Hg : Hg** in 
mercurous salts. This, however, seems improbable because 
the stoichiometric formula, (NH4)2SbBrs, suggests that 
each antimony atom is surrounded by six bromine atoms 
which would prevent close approach of the antimony 
atoms. 
NORMAN ELLIoTT 
Pasadena, California, 
April 11, 1934. 


The Infrared Spectrum of Methyl Deuteride 


As a consequence of its high degree of geometric sym- 
metry the methane molecule, with nine degrees of freedom 
of vibration, has only two active fundamental absorption 
bands, each triply degenerate. While the character of these 
bands justifies the tetrahedral model, they yield very 
little information regarding the dimensions and elastic 
constants of the system. The replacement of one of the 
hydrogens by an atom of deuterium introduces a preferred 
axis, with respect to which the molecule is symmetrical, 
and increases the number of active vibrations from two 
to six. Of these, three are single, introducing electric 
moments parallel to the axis, and three are double, giving 
rise to perpendicular bands. 

We have observed all of these fundamental bands, and 
have determined the rotational structure of one of them, 
the lowest frequency parallel vibration. This consists of a 
strong and narrow zero branch at approximately 1090 
cm~!, and the usual positive and negative branches, the 
intervals between rotation lines being about 4.5 cm™. 
From this interval the maximum moment of inertia may be 
computed at once, J~ 12.28 X 10~*°, and the H-H distances 
are about 2.51X10-§ cm. More bands must be analyzed 
and the molecular constants evaluated before these dimen- 
sions can be given with precision. Assuming the atomic 
distances to be the same for ordinary methane, its moment 
of inertia comes out approximately 10.45<10-*. As 


would be expected, this value is about twice as large as 
that determined from the rotational Raman effect by 
Dickinson, Dillon and Rasetti,! since they assumed that 
the successive lines have J values differing by two. They 
found J=5.17X10~. 

Methyl deuteride is of especial interest in the study of 
rotational and vibrational interaction in perpendicular 
bands, since it is the first molecule available in which both 
moments of inertia may be determined precisely. Because 
the bands lie rather close together and overlap con- 
siderably, their analysis must await further study, and 
will be reported later. 

A very interesting test of the assumption that the 
atomic distances are the same for both isotopes of hydrogen 
will be provided by the study of the methane molecule 
with three hydrogen atoms replaced by deuterium. These 
measurements we hope to publish within a short time, as 
well as the fundamental frequencies of tetra-deutero 
methane. 

E. F. BARKER 
NATHAN GINSBURG 
University of Michigan, 
April 13, 1934. 


1 Dickinson, Dillon and Rasetti, Phys. Rev. 34, 582 
(1929). 


Resolution of Crossed Bond Eigenfunctions in Terms of the Uncrossed Set 


The method of repeatedly applying the theorem for 
uncrossing any pair of bonds expressed by the relation 


x=||4+= 


is adequate for resolving every bond eigenfunction in 
terms of the linearly independent set. However, the labor 
of doing this rises rapidly as the number of bonds in- 
creases. Since in obtaining bases for irreducible representa- 
tions! the resolution of crossed bond eigenfunctions is 
practically unavoidable, some simplification seems neces- 
sary. Very simple rules have been found for the cases 
usually encountered, i.e., those involving 5 bonds or 
fewer, 

We assume we have a crossed bond eigenfunction A with 
2norbitals arranged symmetrically on a circle connected in 
pairs by m bonds. 


Rule I 


The number of times this process of uncrossing will 
yield a particular bond eigenfunction B with no crossed 
bonds equals the number of different ways a set of m con- 
tinuous, entirely noncoincident (except perhaps at points), 
noncrossing paths may be drawn along the lines in A, 
such that they connect each pair of orbitals bound to- 
gether in B. Ordinarily this rule enables one to determine 
coefficients 0 or 1 by inspection, but in any case the 
determination is very simple. A subsidiary rule is necessary 
in cases yielding larger coefficients. Such cases arise when 
there is a point common to more than three paths. 


1 Eyring, Frost and Turkevich, J. Chem. Phys. 1, 777 
(1933). 
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Rule II 


Draw the shortest set of » paths of the kind described in 
rule I. If there is a point common to p paths (only points 
for which is 3 or greater are to be considered) the number 
of times B will occur in uncrossing A is (p—2). In de- 
termining » we note that frequently for an eigenfunction 
A parts of the lines in A are not traversed by the set of 
paths. All points connected by untraversed lines are for 
this purpose to be counted as one point. Having done this 
one further procedure may be necessary in applying rule IT. 
Any path which passes through a point and which with the 
perimeter circle separates other paths common to the 
point into two groups is omitted in determining p for the 
point. This discarding is continued until there are no 
more such paths. 
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The coefficients for all uncrossed bonds are positive if 
the convention is adopted in every B and A of arranging 
the orbitals in order around a circle and always directing 
the arrow representing a bond toward the orbitals with the 
higher ordinal number. 

These rules have been checked for all possible patterns 
of 5 bonds or less. A more complicated rule is required for 
direct application to more than 5 bonds. However, eigen- 
functions involving 6 or more bonds are readily reduced by 
comparatively few uncrossings to cases to which our 
simple rules apply. 

HENRY EyRinG 
CHENG E. Sun 
Frick Chemical Laboratory, 
Princeton University, 
April 13, 1934. 


On Diffusion Equilibria of Electrolytes 


A study of the fundamental elements contained in the 
theory of electrolytic conductivity which will be reported 
later, has led among others, to some results about the 
limiting case, i.e., diffusion equilibrium under the influence 
of external forces. 

Consider for the purpose of illustration the distribution 
of charged colloidal particles and ions in a gravitational 
field. If there are only very few ions present which carry 
charges of the same sign as those present on the colloidal 
particles, the problem can be solved with sufficient accuracy 
for dilute solutions. It is found that there occur large 
deviations from the Boltzmann distribution which tend to 
make the latter more uniform. At the same time the rate 
of establishment of equilibrium is greatly reduced. An 
increase in concentration tends to diminish the deviation 


from the Boltzmann distribution. The presence of jons 
carrying charges of the same sign as the colloidal particles 
acts in the same direction. 

Similarly it can be shown that the diffusion potential in 
electrolytes does not assume its classical value. 

These effects can perhaps be described by the statement 
that the activity coefficients of the charged particles lead 
to an apparent reduction in the value of their diffusion 
coefficient. 

The writer hopes to be able to present shortly a more 
detailed account of these investigations. 

Otto HALPERN 

New York University, 

University Heights, 
April 14, 1934. 


Observations on the Synthesis of Tetra-Deutero Methane 


In the course of the synthesis of tetra-deutero methane 
with the object of studying its physical properties, we 
have made some observations on the relative velocities of 
reactions between the proto water and aluminum carbide 
and the deutero water and aluminum carbide. 

Carefully prepared aluminum carbide in fairly finely 
divided form was placed in a sealed container fitted with 
a small reflux condenser, and the proto water added at 
low temperatures. The whole system was evacuated and 
the mixture allowed to come to room temperature. Gas 
bubbles formed at ordinary temperatures and in the 
course of 10 or 15 minutes raised the pressure in the 
apparatus to atmospheric. The mixture was then raised to 
about 80°C and the methane generated collected over 
mercury. It required on the average 2 minutes to collect 
100 cc of methane. 

With the same apparatus and aluminum carbide from 
the same batch and of the same fineness of division, the 
same method was attempted using the deutero water. 
At room temperatures, no reaction took place in the 
course of two hours. We then filled the apparatus with 
CO, at atmospheric pressure and raised the temperature 


to 80°. It now required 46.5 minutes to collect the first 
100 cc of methane gas. This gives a ratio of the velocities 
for the two reactions of approximately 23. This seems to 
be the highest ratio so far observed for reactions involving 
compounds of the two hydrogens. Not until the reaction 
vessel was heated to about 65° were we able to see any 
bubbles of methane form in the reaction vessel. The next 
100 cc were collected at 90°. At this temperature, the rate 
of evolution of the gas was about three times that at 80°. 

The temperatures were controlled by maintaining the 
temperature of a water bath within a degree or two of the 
temperatures given. There is a possibility, of course, of 
slight error in this way, but that there is an enormous 
difference in the velocities of the reactions cannot be 
doubted. 

Harotp C. UREY 
DoNALD PRICE 
Department of Chemistry, 
Columbia University, 
New York, New York, 
April 16, 1934. 
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LETTERS TO THE EDITOR 


The Diffusion of Hydrogen Through Metals: Fractionating the Hydrogen Isotopes 


Dr. L. B. Tuckerman of the Bureau of Standards, 
referring to our note appearing on page 105, February, 
1934, Journal of Chemical Physics, has called our attention 
to the fact that diffusion of nascent hydrogen through 
metals is nothing new. We certainly did not want to give 
that impression. We only observed what to us was a 
striking example of a well-known phenomenon. Particularly 
noteworthy are the researches of Sainte-Claire Deville 
(Comptes Rendus 59, 102 (1864)), Cailletet (Comptes 
Rendus 66, 847 (1868)), and Graham (Proc. Roy. Soc. 
Al6, 422 (1868)), and the later work by Bellati and 
Lusana (Atti. Instit. Veneto (7) 1, 1173 (1890); 2, 987 
(1891)). The work of Cailletet is especially interesting. 
He allowed electrolytic hydrogen to pass through an 
iron tube and developed pressures of the order of 20 
atmospheres. 


Our note was specifically directed to the possibility of 
separating the hydrogen isotopes by means of selective 
diffusion of these through metals at room temperature. 
“We decided to determine whether or not the hydrogen 
thus collected within the metal pillow was substantially 
the lighter of the two hydrogen isotopes.”’ If we have 
given our readers the impression that we thought that the 
diffusion of electrolytic hydrogen through metals was new, 
we are glad to take this opportunity to correct any such 
impression. 

CoLtn G. FINK 
HAROLD C. UREY 
D. B. LAKE 
Columbia University, 
New York, New York, 
April 10, 1934. 





